Depolarizing-Field Effects in Epitaxial Capacitor Heterostructures.

We identify a transient enhancement of the depolarizing field, leading to an unexpected quench of net polarization, during the growth of a prototypical metal-ferroelectric-metal epitaxial system made of BaTiO_{3} and SrRuO_{3}. Reduced conductivity and, hence, charge screening efficiency in the early growth stage of the SrRuO_{3} top electrode promotes a breakdown of ferroelectric BaTiO_{3} into domains. We demonstrate how a thermal annealing procedure can recover the single-domain state. By tracking the polarization state in situ, using optical second harmonic generation, we bring new understanding to interface-related electrostatic effects in ferroelectric capacitors.

On the effect of confined fluid molecular structure on nonequilibrium phase behaviour and friction.

A detailed understanding of the behaviour of confined fluids is critical to a range of industrial applications, for example to control friction in engineering components. In this study, a combination of tribological experiments and confined nonequilibrium molecular dynamics simulations has been used to investigate the effect of base fluid molecular structure on nonequilibrium phase behaviour and friction. An extensive parameter study, including several lubricant and traction fluid molecules subjected to pressures (0.5-2.0 GPa) and strain rates (104-1010 s-1) typical of the elastohydrodynamic lubrication regime, reveals clear relationships between the friction and flow behaviour. Lubricants, which are flexible, broadly linear molecules, give low friction coefficients that increase with strain rate and pressure in both the experiments and the simulations. Conversely, traction fluids, which are based on inflexible cycloaliphatic groups, give high friction coefficients that only weakly depend on strain rate and pressure. The observed differences in friction behaviour can be rationalised through the stronger shear localisation which is observed for the traction fluids in the simulations. Higher pressures lead to more pronounced shear localisation, whilst increased strain rates lead to a widening of the sheared region. The methods utilised in this study have clarified the physical mechanisms of important confined fluid behaviour and show significant potential in both improving the prediction of elastohydrodynamic friction and developing new molecules to control it.

Boundary-controlled barostats for slab geometries in molecular dynamics simulations.

Molecular dynamics simulation barostat schemes are derived for achieving a given normal pressure for a thin liquid or solid layer confined between two parallel walls. This work builds on the boundary-controlled barostat scheme of Lupkowski and van Swol [J. Chem. Phys. 93, 737 (1990)]. Two classes of barostat are explored, one in which the external load is applied to a virtual regular lattice to which the wall atoms are bound using a tethering potential. The other type of barostat applies the external force directly to the wall atoms, which are not tethered. The extent to which the wall separation distribution is Gaussian is shown to be an effective measure of the quality of the barostat. The first class of barostat can suffer from anomalous dynamical signatures, even resonances, which are sensitive to the effective mass of the virtual lattice, whose value lacks any rigorous definition. The second type of barostat performs much better under equilibrium and wall-sliding nonequilibrium conditions and in not being so prone to resonance instabilities in the wall separation and does not require so many largely arbitrary parameters. The results of exploratory simulations which characterize the dynamical response of the model systems for both dry and wet or lubricated systems using the different barostats are presented. The barostats which have an inherent damping mechanism, such as the ones analogous to a damped harmonic oscillator, reduce the occurrence of large fluctuations and resonances in the separation between the two walls, and they also achieve a new target pressure more quickly. Near a nonequilibrium phase boundary the attributes of the barostat can have a marked influence on the observed behavior.

Molecular dynamics studies of the dissociated screw dislocation in silicon.

Characterizing the motion of dislocations through covalent, high Peierls barrier materials is a key problem in materials science, while despite the progress in experimental studies the actual observation of the atomistic behaviour involved in core migration remains limited. We have applied a hybrid embedding scheme to investigate the dissociated screw dislocation in silicon, consisting of two 30° partials separated by a stacking fault ribbon, under the influence of a constant external strain. Our 'learn on the fly' hybrid technique allows us to calculate the forces on atoms in the vicinity of the core region using the tight binding Kwon potential, whilst the remainder of the bulk matrix is treated within a classical approximation. Applying a 5% strain to the dissociated screw dislocation, for a simulation time of 100 ps at a temperature of 600 K, we observe movement of the partials through two different mechanisms: double kink formation and square ring diffusion at the core. Our results suggest that in these conditions, the role of solitons or anti-phase defects in seeding kink formation and subsequent migration is an important one, which should be taken into account in future studies.