Three-dimensional imaging of vortex structure in a ferroelectric nanoparticle driven by an electric field.

Topological defects of spontaneous polarization are extensively studied as templates for unique physical phenomena and in the design of reconfigurable electronic devices. Experimental investigations of the complex topologies of polarization have been limited to surface phenomena, which has restricted the probing of the dynamic volumetric domain morphology in operando. Here, we utilize Bragg coherent diffractive imaging of a single BaTiO3 nanoparticle in a composite polymer/ferroelectric capacitor to study the behavior of a three-dimensional vortex formed due to competing interactions involving ferroelectric domains. Our investigation of the structural phase transitions under the influence of an external electric field shows a mobile vortex core exhibiting a reversible hysteretic transformation path. We also study the toroidal moment of the vortex under the action of the field. Our results open avenues for the study of the structure and evolution of polar vortices and other topological structures in operando in functional materials under cross field configurations.Imaging of topological states of matter such as vortex configurations has generally been limited to 2D surface effects. Here Karpov et al. study the volumetric structure and dynamics of a vortex core mediated by electric-field induced structural phase transition in a ferroelectric BaTiO3 nanoparticle.

Morphology dictated heterogeneous dynamics in two-dimensional aggregates.

Particulate aggregates occur in a variety of non-equilibrium steady-state morphologies ranging from finite-size compact crystalline structures to non-compact string-like conformations. This diversity is due to the competition between pair-wise short range attraction and long range repulsion between particles. We identify different microscopic mechanisms in action by following the simulated particle trajectories for different morphologies in two dimensions at a fixed density and temperature. In particular, we show that the compact clusters are governed by symmetric caging of particles by their nearest neighbors while sidewise asymmetric binding of particles leads to non-compact aggregates. The measured timescales for these two mechanisms are found to be distinctly different providing phenomenological evidence of a relation between microstructure and dynamics of particulate aggregates. Supporting these findings, the time dependent diffusivity is observed to differ across the morphological hierarchy, while the average long-time dynamics is, in general, sub-diffusive at 'low' temperatures. Finally, one generic relation between diffusivity and structural randomness, applicable to simple equilibrium systems, is validated for complex aggregate forming systems through further analysis of the same system at different temperatures.

Machine learning bandgaps of double perovskites.

The ability to make rapid and accurate predictions on bandgaps of double perovskites is of much practical interest for a range of applications. While quantum mechanical computations for high-fidelity bandgaps are enormously computation-time intensive and thus impractical in high throughput studies, informatics-based statistical learning approaches can be a promising alternative. Here we demonstrate a systematic feature-engineering approach and a robust learning framework for efficient and accurate predictions of electronic bandgaps of double perovskites. After evaluating a set of more than 1.2 million features, we identify lowest occupied Kohn-Sham levels and elemental electronegativities of the constituent atomic species as the most crucial and relevant predictors. The developed models are validated and tested using the best practices of data science and further analyzed to rationalize their prediction performance.

Classification of octet AB-type binary compounds using dynamical charges: A materials informatics perspective.

The role of dynamical (or Born effective) charges in classification of octet AB-type binary compounds between four-fold (zincblende/wurtzite crystal structures) and six-fold (rocksalt crystal structure) coordinated systems is discussed. We show that the difference in the dynamical charges of the fourfold and sixfold coordinated structures, in combination with Harrison's polarity, serves as an excellent feature to classify the coordination of 82 sp-bonded binary octet compounds. We use a support vector machine classifier to estimate the average classification accuracy and the associated variance in our model where a decision boundary is learned in a supervised manner. Finally, we compare the out-of-sample classification accuracy achieved by our feature pair with those reported previously.

Classification of ABO3 perovskite solids: a machine learning study.

We explored the use of machine learning methods for classifying whether a particular ABO3 chemistry forms a perovskite or non-perovskite structured solid. Starting with three sets of feature pairs (the tolerance and octahedral factors, the A and B ionic radii relative to the radius of O, and the bond valence distances between the A and B ions from the O atoms), we used machine learning to create a hyper-dimensional partial dependency structure plot using all three feature pairs or any two of them. Doing so increased the accuracy of our predictions by 2-3 percentage points over using any one pair. We also included the Mendeleev numbers of the A and B atoms to this set of feature pairs. Doing this and using the capabilities of our machine learning algorithm, the gradient tree boosting classifier, enabled us to generate a new type of structure plot that has the simplicity of one based on using just the Mendeleev numbers, but with the added advantages of having a higher accuracy and providing a measure of likelihood of the predicted structure.

A minimal description of morphological hierarchy in two-dimensional aggregates.

A dimensionless parameter Λ is proposed to describe a hierarchy of morphologies in two-dimensional (2D) aggregates formed due to varying competition between short-range attraction and long-range repulsion. Structural transitions from finite non-compact to compact to percolated structures are observed in the configurations simulated by molecular dynamics at a constant temperature and density. Configurational randomness across the transition, measured by the two-body excess entropy S2, exhibits data collapse with the average potential energy [small epsilon, Greek, macron] of the systems. Independent master curves are presented among S2, the reduced second virial coefficient B2* and Λ, justifying this minimal description. This work lays out a coherent basis for the study of 2D aggregate morphologies relevant to diverse nano- and bio-processes.

Mechanical loss in multiferroic materials at high frequencies: friction and the evolution of ferroelastic microstructures.

Energy absorption in multiferroic materials stems typically from strain relaxation which can be strong even when no extrinsic defects exist in the material. Computer simulations of a simple two-dimensional model on a generic, proper ferroelastic material identify the dissipative mechanisms associated with the dynamical motion as: a) advance and retraction of needle-shaped twin domains and, b) movement of kinks inside twin boundaries. Both movements involve friction losses.

Proposal for manipulating functional interface properties of composite organic semiconductors with addition of designed macromolecules.

The arrangement of the electronic levels in an interface between organic semiconductors is crucial for the operation of devices such as solar cells and light emitting diodes. With the addition of designed macromolecules, we show that it is possible to control the relative position of the highest occupied molecular orbital and lowest unoccupied molecular orbital levels, and consequently improve the performance. The designed macromolecules consist of two end segments, each compatible with one of the interface components, and a central segment which adds functionality to the interface. The tails control the position and the orientation of the functional units. When the central functional unit is an electric dipole, an electrostatic field is created due to the orientation of the dipoles, which shifts the electronic levels in a controlled way. We develop a theoretical framework, based on self-consistent field theory, to study the concentration and the orientation of the central functional units. We find that the levels can shift by as much as several tenths of an eV.

High junction and twin boundary densities in driven dynamical systems.

A novel mechanism for the generation of device materials with very high domain boundary densities is described: we shear the sample in a computer experiment and achieve higher twin densities than in rapid quench. These domain patterns are very stable. Elastically soft materials (image with 6.4$ \times $10(5) atoms) has greater twin densities than hard materials, even for nano-crystals.

Magnetic nonuniformity and thermal hysteresis of magnetism in a manganite thin film.

We measured the chemical and magnetic depth profiles of a single crystalline (La(1-x)Pr(x))(1-y)Ca(y)MnO(3-δ) (x=0.52±0.05, y=0.23±0.04, δ=0.14±0.10) film grown on a NdGaO(3) substrate using x-ray reflectometry, electron microscopy, electron energy-loss spectroscopy, and polarized neutron reflectometry. Our data indicate that the film exhibits coexistence of different magnetic phases as a function of depth. The magnetic depth profile is correlated with a variation of chemical composition with depth. The thermal hysteresis of ferromagnetic order in the film suggests a first-order ferromagnetic transition at low temperatures.

First-principles study of stability of the bcc and ω phases of a low Al concentration Nb1-xAlx alloy.

The phase stability and site occupancy of bcc (body centered cubic) Nb(5)Al and slightly rearranged atomic structures have been examined by means of first-principles calculations. In order to use first-principles methods, a periodic cell is required and we used ordered Nb(5)Al compounds as a tractable example of a low Al concentration Nb(1 - x)Al(x) alloy (in this case, for about 17 at.% Al). The instability against an ω-structure atomic displacement was also studied, since this structure is detrimental to ductility. Mulliken population analysis was used to provide an understanding of the hybridization between the atoms and the electronic origin of the site occupancy and instability of the underlying bcc structures. By making calculations for several different configurations of the Nb-Al system we estimated the strengths of the Nb-Nb and Nb-Al bonds. It is shown that the stability of the underlying bcc phases is directly related to Nb-Nb and Nb-Al first-nearest-neighbor interactions. The first-principles calculations were extended to finite temperature by including various contributions to the free energy. In particular, the vibrational free energy was calculated within the quasiharmonic approximation, and it is shown that the contribution of the low energy modes to the lattice entropy helps to stabilize ordered bcc phases against ω-type phase transformations. Semi-quasi-random structures were employed to study the stability of the ordered and disordered bcc phases. Our study showed, in agreement with experiment, that the ω, ordered, and disordered phases can coexist in a nonequilibrium state at finite temperature.

Phonon mechanisms and transformation paths in Pu.

A long-standing problem in Pu science is the crystallographic mechanism for the delta-->alpha' (fcc-->monoclinic) transformation. Orientation relations between the two structures impose severe restrictions on the possible mechanisms and require the transition to be reconstructive, which we describe as a sequence of three displacive transitions: fcc-->trigonal-->hexagonal-->monoclinic. We predict instabilities along the Lambda and Sigma branches in the phonon dispersion of the delta phase and formulate a free energy to describe the displacement of atoms across the transition. We suggest that the delta-->alpha' transition in Pu lies at the threshold of a change in character of the orientation relationship from lighter to heavier actinides, correlating with changes in electron itinerancy, magnetism, and volume.

Stress distributions in diblock copolymers.

We demonstrate how a generalized self-consistent field theory for polymer melts that includes elastic stress and strain fields can be applied to the study of AB diblock copolymers melts. By obtaining the stress distributions for volume conserving strain loadings where lamellar and hexagonal morphologies are stable, we show that the local stress is reduced at the domain interface but slightly enhanced in the immediate vicinity of the interface. The overall stress profile is the result of the combined effects of chain connectivity across the interface, which yields a positive contribution, and the immiscible nature of the monomers, which leads to a stress reduction because of interfacial tension.

Pinning frequencies of the collective modes in alpha-uranium.

Uranium is the only known element that features a charge-density wave (CDW) and superconductivity. We report a comparison of the specific heat of single-crystal and polycrystalline alpha-uranium. In the single crystal we find excess contributions to the heat capacity at 41 K, 38 K, and 23 K, with a Debye temperature ThetaD = 265 K. In the polycrystalline sample the heat capacity curve is thermally broadened (ThetaD = 184 K), but no excess heat capacity was observed. The excess heat capacity Cphi (taken as the difference between the single-crystal and polycrystal heat capacities) is well described in terms of collective-mode excitations above their respective pinning frequencies. This attribution is represented by a modified Debye spectrum with two cutoff frequencies, a pinning frequency V0 for the pinned CDW (due to grain boundaries in the polycrystal), and a normal Debye acoustic frequency occurring in the single crystal.

Simulations of spinodal nucleation in systems with elastic interactions.

Systems with long-range interactions quenched into a metastable state near the pseudospinodal exhibit nucleation that is qualitatively different from classical nucleation near the coexistence curve. We observe nucleation droplets in Langevin simulations of a two-dimensional model of martensitic transformations and determine that the structure of the nucleating droplet differs from the stable martensite structure. Our results, together with experimental measurements of the phonon dispersion curve, allow us to predict the nature of the droplet. The results have implications for nucleation in many solid-solid transitions and the structure of the final state.

Adaptability and "intermediate phase" in randomly connected networks.

We present a simple model that enables us to analytically characterize a floppy to rigid transition and an associated self-adaptive intermediate phase in a random bond network. In this intermediate phase, the network adapts itself to lower the stress due to constraints. Our simulations verify this picture. We use these insights to identify applications of these ideas in computational problems such as vertex cover and K-satisfiability.

Improved convergence in block copolymer self-consistent field theory by Anderson mixing.

A modification to real space polymeric self-consistent field theory algorithms that greatly improves the convergence properties is presented. The method is based on Anderson mixing [D. G. Anderson, J. Assoc. Comput. Mach. 12, 547 (1965)], and each iteration computed takes negligibly longer to perform than with other methods, but the number of iterations required to reach a high accuracy solution is greatly reduced. No a priori knowledge of possible phases is required to apply this method. We apply our approach to a standard diblock copolymer melt, and demonstrate iteration reductions of more than a factor of 5 in some cases.

Elastic moduli of multiblock copolymers in the lamellar phase.

We study the linear elastic response of multiblock copolymer melts in the lamellar phase, where the molecules are composed of tethered symmetric AB diblock copolymers. We use a self-consistent field theory method, and introduce a real space approach to calculate the tensile and shear moduli as a function of block number. The former is found to be in qualitative agreement with experiment. We find that the increase in bridging fraction with block number, that follows the increase in modulus, is not responsible for the increase in modulus. It is demonstrated that the change in modulus is due to an increase in mixing of repulsive A and B monomers. Under extension, this increase originates from a widening of the interface, and more molecules pulled free of the interface. Under compression, only the second of these two processes acts to increase the modulus.

Ordering mechanisms in triblock copolymers.

The ordering mechanisms for an ABC triblock copolymer system are studied using self-consistent field theory. We find a two-phase mechanism, similar to what has been suggested experimentally (two-step mechanism). Analysis of free energy components shows that the two-phase process comes about through a competition between stretching energy and interfacial energy. The mechanism is found to be sufficiently robust so as to make it potentially useful for device applications.

Strain-induced metal-insulator phase coexistence in perovskite manganites.

The coexistence of distinct metallic and insulating electronic phases within the same sample of a perovskite manganite, such as La(1-x-y)Pr(y)Ca(x)MnO3, presents researchers with a tool for tuning the electronic properties in materials. In particular, colossal magnetoresistance in these materials--the dramatic reduction of resistivity in a magnetic field--is closely related to the observed texture owing to nanometre- and micrometre-scale inhomogeneities. Despite accumulated data from various high-resolution probes, a theoretical understanding for the existence of such inhomogeneities has been lacking. Mechanisms invoked so far, usually based on electronic mechanisms and chemical disorder, have been inadequate to describe the multiscale, multiphase coexistence within a unified picture. Moreover, lattice distortions and long-range strains are known to be important in the manganites. Here we show that the texturing can be due to the intrinsic complexity of a system with strong coupling between the electronic and elastic degrees of freedom. This leads to local energetically favourable configurations and provides a natural mechanism for the self-organized inhomogeneities over both nanometre and micrometre scales. The framework provides a physical understanding of various experimental results and a basis for engineering nanoscale patterns of metallic and insulating phases.