The role of interlayer gases and surface asperities in compression-induced intermetallic formation in Ni/Al nanocomposites.

Reactive composites comprising alternating nano- or microscale layers of Ni and Al are known to undergo self-sustaining alloying reactions under compression loading, however the effect of infiltrated gas within the microstructure of such reactive nanolaminates-as well as the presence of asperities on the free surfaces of such composites-is not well understood. This work presents atomistic molecular dynamics simulation and analysis of the mechanical dynamics and thermal evolution of planar Ni/Al nanolaminates under a variety of scenarios of layer dimensions, surface asperity shape and orientation, and interlayer gas identity and concentration. These simulations indicate that the rate of the alloying reaction is inversely correlated with the layer width of the nanolaminate, recapitulating experimental results. The presence of surface asperities of comparable scale to the nanolaminate layer thickness enhances short-term intermetallic mixing but has a marginal accelerant effect on the reaction. Interlayer argon gas acts as a mechanical interferent to reaction, whilst interlayer nitrogen gas-modelled here with a novel interatomic potential-is shown to enhance heat production. These calculations also characterise compression wave dynamics in compression-loaded Ni/Al nanolaminates in greater detail than prior studies and illustrate significant qualitative and quantitative differences between extant embedded atom model (EAM) parameterisations of Ni/Al. Speeds of sound for each metal and EAM are also reported. These differences have implications for the interpretation and comparability of EAM-based modelling of Ni/Al reactions moving forward, and may also have wider implications for EAM modelling of intermetallic systems in general.

The Hitchhiker's Guide to the Wave Function∥.

The electronic wave function of molecules is 3N-dimensional and inseparable in the coordinates of the N electrons. Whereas molecular orbitals are often invoked to visualize the electronic structure, they are nonunique, with the same 3N-dimensional wave function being represented by an infinite number of 3-D, one-electron functions (orbitals). Furthermore, multireference wave functions cannot be described by an antisymmetrized product of a single set of occupied orbitals. What is required is a way to visualize the full dimensionality of the wave function, including the effects of correlation, as a 3N-dimensional being would be able to do. In the past 5 years, we have been developing a way to analyze and visualize highly dimensional wave functions by focusing on the structure of the repeating unit demanded by fermionic behavior. This 3N-dimensional repeating unit, the wave function "tile", can be projected onto the three dimensions of each electron, in turn, to reveal the complete electronic structure. It is found that the tile reproduces canonical chemical motifs such as core-electrons, single bonds and lone pairs. Multiple bonds emerge as the "banana" bonds favored by Pauling. As a function of the reaction coordinate, electron motions are visualized that correspond to the curly arrow notation of organic chemists. Excited states can also be inspected. Analyzing a wave function in terms of fermionic tiling allows for insight not facilitated by the inspection of orbitals or configuration interaction vectors: The wave function tiles of resonance structures reveal that electron correlation in benzene pushes opposing spin electrons to occupy alternate Kekulé structures, and in C2, the emerging structure supports the notion of a triply bonded structure with a weak, fourth bonding contribution.

Prominent out-of-plane diffraction in helium scattering from a methyl-terminated Si(111) surface.

Due to their electrochemical and oxidative stability, organic-terminated semiconductor surfaces are well suited to applications in, for example, photoelectrodes and electrochemical cells, which explains the lively interest in their detailed characterization. Helium atom scattering (HAS) is a useful tool to carry out such characterization. Here, we have simulated HAS in He/CH3-Si(111) based on density functional theory (DFT) potential energy surfaces (PESs) and multi-configuration time-dependent Hartree (MCTDH) dynamics. Our analysis of HAS shows that most diffraction taking place in this system corresponds to high-order out-of-plane peaks. This is a general trend that does not depend on the specific features of the simulations, such as the inclusion or not of the van der Waals long-range effects. This is the first and only He-surface system for which such huge out-of-plane diffraction has been described. This striking theoretical finding should encourage new experimental developments to confirm this previously unreported effect.

Non-adiabatic quantum molecular dynamics by the basis expansion leaping multi-configuration Gaussian (BEL MCG) method: Multi-set and single-set formalisms.

Non-adiabatic transitions are quite often of critical importance in chemical reactions. We have recently developed the basis expansion leaping multi-configuration Gaussian (BEL MCG) method to obtain time-propagated wave packets describing multidimensional reactive molecular systems such as quantum tunneling [T. Murakami and T. J. Frankcombe, J. Chem. Phys. 149, 134113 (2018)]. In this work, we develop BEL MCG for multiple electronic state problems. We present two formalisms for the BEL MCG description of multi-state wave packets, namely, "multi-set" and "single-set." We pay particular attention to investigate what is required to yield accurate dynamics. When there is low population on an electronic state, it is important in the "multi-set" case that the reexpression on that electronic state is applied rigorously. The sharing of basis functions in the single-set approach leads to needing a lower number of basis functions than in the multi-set approach, making it preferable for direct dynamics.

Photoactivity and Stability Co-Enhancement: When Localized Plasmons Meet Oxygen Vacancies in MgO.

Durability is still one of the key obstacles for the further development of photocatalytic energy-conversion systems, especially low-dimensional ones. Encouragingly, recent studies show that nanoinsulators such as SiO2 and MgO exhibit substantially enhanced photocatalytic durability than the typical semiconductor p25 TiO2 . Extending this knowledge, MgO-Au plasmonic defect nanosystems are developed that combine the stable photoactivity from MgO surface defects with energy-focusing plasmonics from Au nanoparticles (NPs), where Au NPs are anchored onto monodispersed MgO nanotemplates. Theoretical calculations reveal that the midgap defect (MGD) states in MgO are generated by oxygen vacancies, which provide the main avenues for upward electron transitions under photoexcitation. These electrons drive stable proton photoreduction to H2 gas via water splitting. A synergistic interaction between Au's localized plasmons and MgO's oxygen vacancies is observed here, which enhances MgO's photoactivity and stability simultaneously. Such co-enhancement is attributed to the stable longitudinal-plasmon-free Au NPs, which provide robust hot electrons capable of overcoming the interband transition barrier (≈1.8 eV) to reach proton reduction potential for H2 generation. The demonstrated plasmonic defect nanosystems are expected to open a new avenue for developing highly endurable photoredox systems for the integration of multifunctionalities in energy conversion, environmental decontamination, and climate change mitigation.

Accurate quantum molecular dynamics for multidimensional systems by the basis expansion leaping multi-configuration Gaussian (BEL MCG) method.

Quantum phenomena are quite often of critical importance in chemical reactions. Thus the development of quantum molecular dynamics approaches is required to study the role of quantum effects such as tunnelling in chemical processes. The basis expansion leaping multi-configuration Gaussian (BEL MCG) method has been developed to obtain time-propagated wave packets describing reactive molecular systems. Here we examine the applicability of BEL MCG to double well problems in several dimensions. We pay particular attention to what is required to yield highly accurate dynamics with respect to several key features of the BEL MCG propagation. The importance of using basis functions of a width appropriate to the nature of the potential energy surface in the region of configuration space where each basis function is located is highlighted, which has implications for virtually all quantum molecular dynamics methods utilising Gaussian basis functions.

Challenges facing an understanding of the nature of low-energy excited states in photosynthesis.

While the majority of the photochemical states and pathways related to the biological capture of solar energy are now well understood and provide paradigms for artificial device design, additional low-energy states have been discovered in many systems with obscure origins and significance. However, as low-energy states are naively expected to be critical to function, these observations pose important challenges. A review of known properties of low energy states covering eight photochemical systems, and options for their interpretation, are presented. A concerted experimental and theoretical research strategy is suggested and outlined, this being aimed at providing a fully comprehensive understanding.

Chemical bonding motifs from a tiling of the many-electron wavefunction.

A method is presented to partition the 3N-dimensional space of a many-electron wavefunction into hyper-regions related by permutation symmetry. These hyper-regions represent unit cells, or "tiles" of the wavefunction from which the wavefunction may be regenerated in its entirety upon application of the set of permutations of like-spin electrons. The method, wherein a Voronoi diagram is constructed from the (even permutations of the) average position of a swarm of Monte Carlo walkers sampling |Ψ|(2), determines a self-consistent partitioning of the wavefunction. When one of the identical 3N-dimensional Voronoi sites is projected onto the coordinates of each electron, chemical motifs naturally appear, such as core electrons, lone-pairs, single-bonds and banana-bonds. The structures determined for N2, O2, F2, and other molecules correspond to the double-quartet theory of Linnett. When the procedure is applied to C2, we arrive at an interpretation of its bonding in terms of a near triple bond with singlet-coupled outer electrons.

Explicit calculation of the excited electronic states of the photosystem II reaction centre.

The excited states of sets of the cofactors found in the photosystem II reaction centre have been calculated directly as a multi-monomer supermolecule for the first time. Time-dependent density functional theory was used with the CAM-B3LYP functional. Multiple excited states for each cofactor were found at lower energies than the lowest energy state corresponding to charge transfer states (in which an electron is shifted from one cofactor to another). The electrostatic environment was found to have a dramatic impact on the excited state energies, with the effect of a surrounding dielectric medium being less significant.

H2 Adsorption in a Porous Crystal: Accurate First-Principles Quantum Simulation.

A general method is presented for constructing, from ab initio quantum chemistry calculations, the potential energy surface (PES) for H2 absorbed in a porous crystalline material. The method is illustrated for the metal-organic framework material MOF-5. Rigid body quantum diffusion Monte Carlo simulations are used in the construction of the PES and to evaluate the quantum ground state of H2 in MOF-5, the zero-point energy, and the enthalpy of adsorption at 0 K.

Electron-pinned defect-dipoles for high-performance colossal permittivity materials.

The immense potential of colossal permittivity (CP) materials for use in modern microelectronics as well as for high-energy-density storage applications has propelled much recent research and development. Despite the discovery of several new classes of CP materials, the development of such materials with the required high performance is still a highly challenging task. Here, we propose a new electron-pinned, defect-dipole route to ideal CP behaviour, where hopping electrons are localized by designated lattice defect states to generate giant defect-dipoles and result in high-performance CP materials. We present a concrete example, (Nb+In) co-doped TiO2 rutile, that exhibits a largely temperature- and frequency-independent colossal permittivity (> 10(4)) as well as a low dielectric loss (mostly < 0.05) over a very broad temperature range from 80 to 450 K. A systematic defect analysis coupled with density functional theory modelling suggests that 'triangular' In2(3+)Vo(••)Ti(3+) and 'diamond' shaped Nb2(5+)Ti(3+)A(Ti) (A = Ti(3+)/In(3+)/Ti(4+)) defect complexes are strongly correlated, giving rise to large defect-dipole clusters containing highly localized electrons that are together responsible for the excellent CP properties observed in co-doped TiO2. This combined experimental and theoretical work opens up a promising feasible route to the systematic development of new high-performance CP materials via defect engineering.

Using Hessian update formulae to construct modified Shepard interpolated potential energy surfaces: application to vibrating surface atoms.

Modified Shepard interpolation based on second order Taylor series expansions has proven to be a flexible tool for constructing potential energy surfaces in a range of situations. Extending this to gas-surface dynamics where surface atoms are allowed to move represents a substantial increase in the dimensionality of the problem, reflected in a dramatic increase in the computational cost of the required Hessian (matrix of second derivatives) evaluations. This work demonstrates that using approximate Hessians derived from well known Hessian update formulae and a single accurate Hessian can provide an effective way to avoid this expensive accurate Hessian determination.

Basis expansion leaping: a new method to solve the time-dependent Schrödinger equation for molecular quantum dynamics.

A wide variety of molecular systems that have recently come into the reach of experimental and theoretical investigation is dominated by quantum phenomena. However, even state of the art quantum propagation techniques are either unsuitable for general application to molecular systems with strong interference and tunneling characteristics or are computationally prohibitive for systems with more than a few degrees of freedom. In this Letter, we introduce a novel quantum propagation technique with wide applicability, controllable accuracy, and efficient utilization of computational resources. Its performance is validated for tunneling and dissociating systems with 1, 2, and 3 degrees of freedom, and the scaling behavior with respect to system dimensionality and requested accuracy is discussed.

Comparative first-principles structural and vibrational properties of rutile and anatase TiO2.

The structural and vibrational properties of two polymorphs of TiO2, rutile and anatase, have been investigated by first-principles methods at different levels of exchange-correlational (XC) energy functionals in density functional theory (DFT). Reports in the literature to date are contradictory regarding the stability of the rutile phase using DFT XC-functionals more sophisticated than simple local-density approximation. Here the PBEsol generalized gradient approximation (GGA), TPSS meta-GGA, and HSE06 hybrid functionals have been employed to demonstrate the XC-functional effects on the calculated structural, phonon and thermodynamic properties of rutile and anatase TiO2. Lattice and elastic parameters correctly calculated with these XC-functionals show good agreement with the experimental values. Calculated phonon frequencies generated stable phonon dispersion relations for both rutile and anatase TiO2when correctly converged, in agreement with the experimental observations. The phonon frequencies along high symmetry Brillouin zone paths and their corresponding phonon density of states showed sensitivity to different levels of XC-functional employed in phonon dispersion prediction. Nevertheless, the thermodynamic properties of rutile and anatase TiO2estimated by harmonic approximations are in excellent experimental agreement and are effectively invariant to the level of theory employed in the DFT XC-functional.

Reactive Bimetallic Nanostructures Based on Triply Periodic Minimal Surfaces: A Molecular Dynamics Study toward the Limits of Performance.

A variety of intermetallic compounds possesses high enthalpies of formation. These compounds may be formed from reactive compacts or nanostructures comprised of unreacted precursor metals. These precursor structures support self-propagating high temperature synthesis (SHS) reactions which afford very high specific energy densities and rates, with excellent spatial control and a variety of useful applications. The present work compares the reactivity of notional bimetallic nanostructures based on well-known triply periodic minimal surfaces (TPMSes) with the popular reactive nanolaminate (RNL) modality for the Ni/Al system, using a molecular dynamics approach. TPMS-derived nanostructures were found to have lower ignition energies and faster reaction rates than RNLs of comparable periodicity, while the maximum achievable temperature of ignitions was found to be modulated by a complex interplay of factors including reaction rate and specific metal/metal interface density. Nanostructure reactivity and thermochemistry is also affected by effective diffusion dimensionality and recalescent precipitation of intermetallic crystallites. The TPMS-derived reactive nanostructures presented herein anticipate plausible advances in nanofabrication technology.

An adiabatic capture theory and quasiclassical trajectory study of C + NO and O + CN on the 2A', 2A", and 4A" potential energy surfaces.

The adiabatic capture centrifugal sudden approximation (ACCSA) has been applied to the C + NO and O + CN reactions, along with quasiclassical trajectory simulations. Existing global analytic fits to the potential energy surfaces of the CNO system in the (2)A', (2)A", and (4)A" electronic states have been used. Thermal rate constants for reaction in each of the electronic states have been calculated. In all cases a strong temperature dependence is evident in the calculated rate constants. The agreement between the calculated adiabatic capture and quasiclassical trajectory rate constants is excellent in some cases, but these rate constants differ considerably in other cases. This behavior is analyzed in terms of the anisotropy of the potential energy surfaces. On the basis of this analysis, we propose a new diagnostic for the reliability of ACCSA capture calculations.

Modified Shepard interpolation of gas-surface potential energy surfaces with strict plane group symmetry and translational periodicity.

A new formulation of modified Shepard interpolation of potential energy surface data for gas-surface reactions has been developed. The approach has been formulated for monoatomic or polyatomic adsorbates interacting with crystalline solid surfaces of any plane group symmetry. The interpolation obeys the two dimensional translational periodicity and plane group symmetry of the solid surface by construction. The interpolation remains continuous and smooth everywhere. The interpolation developed here is suitable for constructing potential energy surfaces by sampling classical trajectories using the Grow procedure. A model function has been used to demonstrate the method, showing the convergence of the classical gas-surface reaction probability.

A new method for screening potential sII and sH hydrogen clathrate hydrate promoters with model potentials.

A new predictive computational method for classifying clathrate hydrate promoter molecules is presented, based on the interaction energies between potential promoters and the water networks of sII and sH clathrates. The motivation for this work is identifying promoters for storing hydrogen compactly in clathrate hydrates. As a first step towards achieving this goal, we have developed a general method aimed at distinguishing between molecules that form sII clathrate hydrates and molecules that can-together with a weakly interacting help gas-form sH clathrate hydrates. The new computational method calculates differences in estimated formation energies of the sII and the sH clathrate hydrate. Model interaction potentials have been used, including the electrostatic interactions with newly calculated partial charges for all the considered potential promoter molecules. The methodology can discriminate between the clathrate structure types (sII or sH) formed by each potential promoter with good selectivity, i.e., better than achieved with a simple van der Waals diameter criterion.

Potential energy surfaces for gas-surface reactions.

A method for constructing the potential energy surface for reactions of a molecule with the surface of cleaved non-conducting crystals is reported. The method uses systematic fragmentation to express the total potential in terms of potential energy surfaces which describe reactions of relatively small molecules in the gas phase. The approach is illustrated by an application to the reaction of hydrogen atoms with a hydrogen-terminated silicon(111) surface.

Hole-Pinned Defect Clusters for a Large Dielectric Constant up to GHz in Zinc and Niobium Codoped Rutile SnO2.

High permittivity materials for a gigahertz (GHz) communication technology have been actively sought for some time. Unfortunately, in most materials, the dielectric constant starts to drop as frequencies increase through the megahertz (MHz) range. In this work, we report a large dielectric constant of ∼800 observed in defect-mediated rutile SnO2 ceramics, which is nearly frequency and temperature independent over the frequency range of 1 mHz to 35 GHz and temperature range of 50-450 K. Experimental and theoretical investigations demonstrate that the origin of the high dielectric constant can be attributed to the formation of locally well-defined Zn2+-Nb4+ defect clusters, which create hole-pinned defect dipoles. We believe that this work provides a promising strategy to advance dipole polarization theory and opens up a direction for the design and development of high frequency, broadband dielectric materials for use in future communication technology.