First-principles energetics of water: a many-body analysis

Standard forms of density-functional theory (DFT) have good predictive power for many materials, but are not yet fully satisfactory for solid, liquid and cluster forms of water. We use a many-body separation of the total energy into its 1-body, 2-body (2B) and beyond-2-body (B2B) components to analyze the deficiencies of two popular DFT approximations. We show how machine-learning methods make this analysis possible for ice structures as well as for water clusters. We find that the crucial energy balance between compact and extended geometries can be distorted by 2B and B2B errors, and that both types of first-principles error are important.

First-principles energetics of water clusters and ice: a many-body analysis.

Standard forms of density-functional theory (DFT) have good predictive power for many materials, but are not yet fully satisfactory for cluster, solid, and liquid forms of water. Recent work has stressed the importance of DFT errors in describing dispersion, but we note that errors in other parts of the energy may also contribute. We obtain information about the nature of DFT errors by using a many-body separation of the total energy into its 1-body, 2-body, and beyond-2-body components to analyze the deficiencies of the popular PBE and BLYP approximations for the energetics of water clusters and ice structures. The errors of these approximations are computed by using accurate benchmark energies from the coupled-cluster technique of molecular quantum chemistry and from quantum Monte Carlo calculations. The systems studied are isomers of the water hexamer cluster, the crystal structures Ih, II, XV, and VIII of ice, and two clusters extracted from ice VIII. For the binding energies of these systems, we use the machine-learning technique of Gaussian Approximation Potentials to correct successively for 1-body and 2-body errors of the DFT approximations. We find that even after correction for these errors, substantial beyond-2-body errors remain. The characteristics of the 2-body and beyond-2-body errors of PBE are completely different from those of BLYP, but the errors of both approximations disfavor the close approach of non-hydrogen-bonded monomers. We note the possible relevance of our findings to the understanding of liquid water.

Graphite intercalation compounds under pressure

Motivated by recent experimental work, we use first-principles density functional theory methods to conduct an extensive search for low enthalpy structures of C$_6$Ca under pressure. As well as a range of buckled structures, which are energetically competitive over an intermediate range of pressures, we show that the high pressure system ($\gtrsim 18$ GPa) is unstable towards the formation of a novel class of layered structures, with the most stable compound involving carbon sheets containing five- and eight-membered rings. As well as discussing the energetics of the different classes of low enthalpy structures, we comment on the electronic structure of the high pressure compound and its implications for superconductivity.

The wave transmission coefficients and coupling loss factors of point connected structures

This analysis is concerned with the calculation of the elastic wave transmission coefficients and coupling loss factors between an arbitrary number of structural components that are coupled at a point. A general approach to the problem is presented and it is demonstrated that the resulting coupling loss factors satisfy reciprocity. A key aspect of the method is the consideration of cylindrical waves in two-dimensional components, and this builds upon recent results regarding the energetics of diffuse wavefields when expressed in cylindrical coordinates. Specific details of the method are given for beam and thin plate components, and a number of examples are presented. © 2002 Acoustical Society of America.

The limits of flexoelectricity in liquid crystals

The flexoelectric conversion of mechanical to electrical energy in nematic liquid crystals is investigated using continuum theory. Since the electrical energy produced cannot exceed the mechanical energy supplied, and vice-versa, upper bounds are imposed on the magnitudes of the flexoelectric coefficients in terms of the elastic and dielectric coefficients. For conventional values of the elastic and dielectric coefficients, it is shown that the flexoelectric coefficients may not be larger than a few tens of pC/m. This has important consequences for the future use of such flexoelectric materials in devices and the related energetics of distorted equilibrium structures. © 2011 Author(s).

High-density remote plasma sputtering of high-dielectric-constant amorphous hafnium oxide films

Hafnium oxide (HfOx) is a high dielectric constant (k) oxide which has been identified as being suitable for use as the gate dielectric in thin film transistors (TFTs). Amorphous materials are preferred for a gate dielectric, but it has been an ongoing challenge to produce amorphous HfOx while maintaining a high dielectric constant. A technique called high target utilization sputtering (HiTUS) is demonstrated to be capable of depositing high-k amorphous HfOx thin films at room temperature. The plasma is generated in a remote chamber, allowing higher rate deposition of films with minimal ion damage. Compared to a conventional sputtering system, the HiTUS technique allows finer control of the thin film microstructure. Using a conventional reactive rf magnetron sputtering technique, monoclinic nanocrystalline HfOx thin films have been deposited at a rate of ∼1.6nmmin-1 at room temperature, with a resistivity of 1013Ωcm, a breakdown strength of 3.5MVcm-1 and a dielectric constant of ∼18.2. By comparison, using the HiTUS process, amorphous HfOx (x=2.1) thin films which appear to have a cubic-like short-range order have been deposited at a high deposition rate of ∼25nmmin-1 with a high resistivity of 1014Ωcm, a breakdown strength of 3MVcm-1 and a high dielectric constant of ∼30. Two key conditions must be satisfied in the HiTUS system for high-k HfOx to be produced. Firstly, the correct oxygen flow rate is required for a given sputtering rate from the metallic target. Secondly, there must be an absence of energetic oxygen ion bombardment to maintain an amorphous microstructure and a high flux of medium energy species emitted from the metallic sputtering target to induce a cubic-like short range order. This HfOx is very attractive as a dielectric material for large-area electronic applications on flexible substrates. A remote plasma sputtering process (high target utilization sputtering, HiTUS) has been used to deposit amorphous hafnium oxide with a very high dielectric constant (∼30). X-ray diffraction shows that this material has a microstructure in which the atoms have a cubic-like short-range order, whereas radio frequency (rf) magnetron sputtering produced a monoclinic polycrystalline microstructure. This is correlated to the difference in the energetics of remote plasma and rf magnetron sputtering processes. © 2013 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.

Computational analysis of liquid crystalline elastomer membranes: Changing Gaussian curvature without stretch energy

Liquid crystalline elastomers (LCEs) can undergo extremely large reversible shape changes when exposed to external stimuli, such as mechanical deformations, heating or illumination. The deformation of LCEs result from a combination of directional reorientation of the nematic director and entropic elasticity. In this paper, we study the energetics of initially flat, thin LCE membranes by stress driven reorientation of the nematic director. The energy functional used in the variational formulation includes contributions depending on the deformation gradient and the second gradient of the deformation. The deformation gradient models the in-plane stretching of the membrane. The second gradient regularises the non-convex membrane energy functional so that infinitely fine in-plane microstructures and infinitely fine out-of-plane membrane wrinkling are penalised. For a specific example, our computational results show that a non-developable surface can be generated from an initially flat sheet at cost of only energy terms resulting from the second gradients. That is, Gaussian curvature can be generated in LCE membranes without the cost of stretch energy in contrast to conventional materials. © 2013 Elsevier Ltd. All rights reserved.

The evolution of a stratified turbulent cloud

Localized regions of turbulence, or turbulent clouds, in a stratified fluid are the subject of this study, which focuses on the edge dynamics occurring between the turbulence and the surrounding quiescent region. Through laboratory experiments and numerical simulations of stratified turbulent clouds, we confirm that the edge dynamics can be subdivided into materially driven intrusions and horizontally travelling internal wave-packets. Three-dimensional visualizations show that the internal gravity wave-packets are in fact large-scale pancake structures that grow out of the turbulent cloud into the adjacent quiescent region. The wave-packets were tracked in time, and it is found that their speed obeys the group speed relation for linear internal gravity waves. The energetics of the propagating waves, which include waveforms that are inclined with respect to the horizontal, are also considered and it is found that, after a period of two eddy turnover times, the internal gravity waves carry up to 16 % of the cloud kinetic energy into the initially quiescent region. Turbulent events in nature are often in the form of decaying turbulent clouds, and it is therefore suggested that internal gravity waves radiated from an initial cloud could play a significant role in the reorganization of energy and momentum in the atmosphere and oceans.©2013 Cambridge University Press.

Bipedal walking and running with spring-like biarticular muscles.

Compliant elements in the leg musculoskeletal system appear to be important not only for running but also for walking in human locomotion as shown in the energetics and kinematics studies of spring-mass model. While the spring-mass model assumes a whole leg as a linear spring, it is still not clear how the compliant elements of muscle-tendon systems behave in a human-like segmented leg structure. This study presents a minimalistic model of compliant leg structure that exploits dynamics of biarticular tension springs. In the proposed bipedal model, each leg consists of three leg segments with passive knee and ankle joints that are constrained by four linear tension springs. We found that biarticular arrangements of the springs that correspond to rectus femoris, biceps femoris and gastrocnemius in human legs provide self-stabilizing characteristics for both walking and running gaits. Through the experiments in simulation and a real-world robotic platform, we show how behavioral characteristics of the proposed model agree with basic patterns of human locomotion including joint kinematics and ground reaction force, which could not be explained in the previous models.