Second generation Car-Parrinello molecular dynamics

Born-Oppenheimer Molecular Dynamics with a Linear Combination of Atomic Orbitals and Hybrid Functionals for Condensed Matter Simulations Made Possible. Theory and Performance for the Microcanonical and Canonical Ensembles

The implementation of an original Born-Oppenheimer molecular dynamics module is presented, which is able to perform simulations of large and complex condensed phase systems for sufficiently long time scales at the level of density functional theory with hybrid functionals, in the microcanonical (NVE) and canonical (NVT) ensembles. The algorithm is fully integrated in the Crystal code, a program for quantum mechanical simulations of materials, whose peculiarity stems from the use of atom-centered basis functions within a linear combination of atomic orbitals to describe the wave function. The corresponding efficiency in the evaluation of the exact Fock exchange series has led to the implementation of a rich variety of hybrid density functionals at a low computational cost. In addition, the molecular dynamics implementation benefits also from the effective MPI parallelization of the code, suited to exploit high-performance computing resources available on current generation supercomputer architectures. Furthermore, the information contained in the trajectory of the dynamics is extracted through a series of postprocessing algorithms that provide the radial distribution function, the diffusion coefficient and the vibrational density of states. In this work, we present a detailed description of the theoretical framework and the algorithmic implementation, followed by a critical evaluation of the accuracy and parallel performance (e.g., strong and weak scaling) of this approach, when ice and liquid water simulations are performed in the microcanonical and canonical ensembles.

Water in PHI Nanopores: Modeling Adsorption, Solvent Structure, and Thermodynamics

We modeled the uptake of water molecules into the nanopores of potassium-polyheptazineimide (K-PHI), a 2D covalent material that is one of the best water-splitting photocatalysts to date possessing experimentally reported strong water binding. In the current models, we find that first water molecules are bound with -94.5 kJ/mol, i.e., twice the cohesion energy of water and one of the highest adsorption enthalpies reported so far. This strong binding proceeds unexpectedly on a similar enthalpy level until the pore is filled, while the binding strength is passed through a conjugated water network. The tight binding is also expressed in calculated, strongly shortened O-O distances, which are on average about 5% shorter than in bulk water, which corresponds to a much higher water density, for a 2D structure above 1.1 g/ cm³. The H-bridges are strongly aligned in the direction perpendicular to the covalent planes, which could give reasons for the experimentally observed ultrahigh ion fluxes and conductivity of K-PHI membranes. Decomposition of the adsorption energy into components reveals an unexpectedly high charge transfer contribution, where the partly naked K⁺ ions play a key role. The latter fact not only offers a new structural lead motif for the design of more strongly, but reversibly binding adsorption materials involving metal ions on their surface but also puts cations as known cofactors in enzymes into a new light.

Reactive Martini: Chemical Reactions in Coarse-Grained Molecular Dynamics Simulations

Chemical reactions are ubiquitous in both materials and the biophysical sciences. While coarse-grained (CG) molecular dynamics simulations are often needed to study the spatiotemporal scales present in these fields, chemical reactivity has not been explored thoroughly in CG models. In this work, a new approach to model chemical reactivity is presented for the widely used Martini CG Martini model. Employing tabulated potentials with a single extra particle for the angle dependence, the model provides a generic framework for capturing bonded topology changes using nonbonded interactions. As a first example application, the reactive model is used to study the macrocycle formation of benzene-1,3-dithiol molecules through the formation of disulfide bonds. We show that starting from monomers, macrocycles with sizes in agreement with experimental results are obtained using reactive Martini. Overall, our reactive Martini framework is general and can be easily extended to other systems. All of the required scripts and tutorials to explain its use are provided online.

A Dynamic Proton Bond: MH+·H2O ⇌ M·H3O+ Interconversion in Loosely Coordinated Environments

The interaction of organic molecules with oxonium cations within their solvation shell may lead to the emergence of dynamic supramolecular structures with recurrently changing host-guest chemical identity. We illustrate this phenomenon in benchmark proton-bonded complexes of water with polyether macrocyles. Despite the smaller proton affinity of water versus the ether group, water in fact retains the proton in the form of H₃O⁺, with increasing stability as the coordination number increases. Hindrance in many-fold coordination induces dynamic reversible (ether)·H₃O⁺ ⇌ (etherH⁺)·H₂O interconversion. We perform infrared action ion spectroscopy over a broad spectral range to expose the vibrational signatures of the loose proton bonding in these systems. Remarkably, characteristic bands for the two limiting proton bonding configurations are observed in the experimental vibrational spectra, superimposed onto diffuse bands associated with proton delocalization. These features cannot be described by static equilibrium structures but are accurately modeled within the framework of ab initio molecular dynamics.

A variational formulation of the Harris functional as a correction to approximate Kohn-Sham density functional theory

Accurate descriptions of intermolecular interactions are of great importance in simulations of molecular liquids. We present an electronic structure method that combines the accuracy of the Harris functional approach with the computational efficiency of approximately linear-scaling density functional theory (DFT). The non-variational nature of the Harris functional has been addressed by constructing a Lagrangian energy functional, which restores the variational condition by imposing stationarity with respect to the reference density. The associated linear response equations may be treated with linear-scaling efficiency in an atomic orbital based scheme. Key ingredients to describe the structural and dynamical properties of molecular systems are the forces acting on the atoms and the stress tensor. These first-order derivatives of the Harris Lagrangian have been derived and implemented in consistence with the energy correction. The proposed method allows for simulations with accuracies close to the Kohn-Sham DFT reference. Embedded in the CP2K program package, the method is designed to enable ab initio molecular dynamics simulations of molecular solutions for system sizes of several thousand atoms. Available subsystem DFT methods may be used to provide the reference density required for the energy correction at near linear-scaling efficiency. As an example of production applications, we applied the method to molecular dynamics simulations of the binary mixtures cyclohexane-methanol and toluene-methanol, performed within the isobaric-isothermal ensemble, to investigate the hydrogen bonding network in these non-ideal mixtures.

Unravelling the Complex LiOH-Based Cathode Chemistry in Lithium-Oxygen Batteries

The LiOH-based cathode chemistry has demonstrated potential for high-energy Li-O₂ batteries. However, the understanding of such complex chemistry remains incomplete. Herein, we use the combined experimental methods with ab initio calculations to study LiOH chemistry. We provide a unified reaction mechanism for LiOH formation during discharge via net 4 e^- oxygen reduction, in which Li₂ O₂ acts as intermediate in low water-content electrolyte but LiHO₂ as intermediate in high water-content electrolyte. Besides, LiOH decomposes via 1 e^- oxidation during charge, generating surface-reactive hydroxyl species that degrade organic electrolytes and generate protons. These protons lead to early removal of LiOH, followed by a new high-potential charge plateau (1 e^- water oxidation). At following cycles, these accumulated protons lead to a new high-potential discharge plateau, corresponding to water formation. Our findings shed light on understanding of 4 e^- cathode chemistries in metal-air batteries.

Molecular Dynamics Simulations of Nitric Oxide Scattering Off Graphene

We performed classical molecular dynamics simulations to model the scattering process of nitric oxide, NO, off graphene supported on gold. This is motivated by our desire to probe the energy transfer in collisions with graphene. Since many of these collision systems comprising of graphene and small molecules have been shown to scatter non-reactively, classical molecular dynamics appear to describe such systems sufficiently. We directed thousands of trajectories of NO molecules onto graphene along the surface normal, while varying impact position, but also speed, orientation, and rotational excitation of the nitric oxide, and compare the results with experimental data. While experiment and theory do not match quantitatively, we observe agreement that the relative amount of kinetic energy lost during the collision increases with increasing initial kinetic energy of the NO. Furthermore, while at higher collision energies, all NO molecules lose some energy, and the vast majority of NO is scattered back, in contrast at low impact energies, the fraction of those nitric oxide molecules that are trapped at the surface increases, and some NO molecules even gain some kinetic energy during the collision process. The collision energy seems to preferentially go into the collective motion of the carbon atoms in the graphene sheet.

"Red Carbon": A Rediscovered Covalent Crystalline Semiconductor

Carbon suboxide (C₃ O₂ ) is a unique molecule able to polymerize spontaneously into highly conjugated light-absorbing structures at temperatures as low as 0 °C. Despite obvious advantages, little is known about the nature and the functional properties of this carbonaceous material. In this work, the aim is to bring "red carbon," a forgotten polymeric semiconductor, back to the community's attention. A solution polymerization process is adapted to simplify the synthesis and control the structure. This allows one to obtain this crystalline covalent material at low temperatures. Both spectroscopic and elemental analyses support the chemical structure represented as conjugated ladder polypyrone ribbons. Density functional theory calculations suggest a crystalline structure of AB stacks of polypyrone ribbons and identify the material as a direct bandgap semiconductor with a medium bandgap that is further confirmed by optical analysis. The material shows promising photocatalytic performance using blue light. Moreover, the simple condensation-aromatization route described here allows the straightforward fabrication of conjugated ladder polymers and can be inspiring for the synthesis of carbonaceous materials at low temperatures in general.

Accurate diffusion coefficients of the excess proton and hydroxide in water via extensive ab initio simulations with different schemes

Despite its simple molecular formula, obtaining an accurate in silico description of water is far from straightforward. Many of its very peculiar properties are quite elusive, and in particular, obtaining good estimations of the diffusion coefficients of the solvated proton and hydroxide at a reasonable computational cost has been an unsolved challenge until now. Here, I present extensive results of several unusually long ab initio molecular dynamics (MD) simulations employing different combinations of the Born-Oppenheimer and second-generation Car-Parrinello MD propagation methods with different ensembles (NVE and NVT) and thermostats, which show that these methods together with the RPBE-D3 functional provide a very accurate estimation of the diffusion coefficients of the solvated H₃O⁺ and OH^- ions, together with an extremely accurate description of several properties of neutral water (such as the structure of the liquid and its diffusion and shear viscosity coefficients). In addition, I show that the estimations of D_H3O⁺ and D_OH^- depend dramatically on the simulation length, being necessary to reach timescales in the order of hundreds of picoseconds to obtain reliable results.

Structure and interactions at the Mg(0001)/water interface: An ab initio study

A molecular level understanding of metal/bulk water interface structure is key for a wide range of processes, including aqueous corrosion, which is our focus, but their buried nature makes experimental investigation difficult and we must mainly rely on simulations. We investigate the Mg(0001)/water interface using second generation Car–Parrinello molecular dynamics (MD) to gain structural information, combined with static density functional theory calculations to probe the atomic interactions and electronic structure (e.g., calculating the potential of zero charge). By performing detailed structural analyses of both metal–surface atoms and the near-surface water, we find that, among other insights: (i) water adsorption causes significant surface roughening (the planar distribution for top-layer Mg has two peaks separated by [Formula: see text]), (ii) strongly adsorbed water covers only [Formula: see text] of available surface sites, and (iii) adsorbed water avoids clustering on the surface. Static calculations are used to gain a deeper understanding of the structuring observed in MD. For example, we use an energy decomposition analysis combined with calculated atomic charges to show that adsorbate clustering is unfavorable due to Coulombic repulsion between adsorption site surface atoms. Results are discussed in the context of previous simulations carried out on other metal/water interfaces. The largest differences for the Mg(0001)/water system appear to be the high degree of surface distortion and the minimal difference between the metal work function and metal/water potential of zero charge (at least compared to other interfaces with similar metal–water interaction strengths). The structural information, in this paper, is important for understanding aqueous Mg corrosion, as the Mg(0001)/water interface is the starting point for key reactions. Furthermore, our focus on understanding the driving forces behind this structuring leads to important insights for general metal/water interfaces.

Biomimetic Nanochannels: From Fabrication Principles to Theoretical Insights

Biological nanochannels which can regulate ionic transport across cell membranes intelligently play a significant role in physiological functions. Inspired by these nanochannels, numerous artificial nanochannels have been developed during recent years. The exploration of smart solid-state nanochannels can lay a solid foundation, not only for fundamental studies of biological systems but also practical applications in various fields. The basic fabrication principles, functional materials, and diverse applications based on artificial nanochannels are summarized in this review. In addition, theoretical insights into transport mechanisms and structure-function relationships are discussed. Meanwhile, it is believed that improvements will be made via computer-guided strategy in designing more efficient devices with upgrading accuracy. Finally, some remaining challenges and perspectives for developments in both novel conceptions and technology of this inspiring research field are stated.

Electronic and Vibrational Properties of Allene Carotenoids

Carotenoids are conjugated linear molecules built from the repetition of terpene units, which display a large structural diversity in nature. They may, in particular, contain several types of side or end groups, which tune their functional properties, such as absorption position and photochemistry. We report here a detailed experimental study of the absorption and vibrational properties of allene-containing carotenoids, together with an extensive modeling of these experimental data. Our calculations can satisfactorily explain the electronic properties of vaucheriaxanthin, where the allene group introduces the equivalent of one C=C double bond into the conjugated C=C chain. The position of the electronic absorption of fucoxanthin and butanoyloxyfucoxanthin requires long-range corrections to be found correctly on the red side of that of vaucheriaxanthin; however, these corrections tend to overestimate the effect of the conjugated and nonconjugated C=O groups in these molecules. We show that the resonance Raman spectra of these carotenoids are largely perturbed by the presence of the allene group, with the two major Raman contributions split into two components. These perturbations are satisfactorily explained by modeling, through a gain in the Raman intensity of the C=C antisymmetric stretching mode, induced by the presence of the allene group in the carotenoid C=C chain.

Full Assignment of Ab-Initio Raman Spectra at Finite Temperatures Using Wannier Polarizabilities: Application to Cyclohexane Molecule in Gas Phase

We demonstrate how to fully ascribe Raman peaks simulated using ab initio molecular dynamics to specific vibrations in the structure at finite temperatures by means of Wannier functions. Here, we adopt our newly introduced method for the simulation of the Raman spectra in which the total polarizability of the system is expressed as a sum over Wannier polarizabilities. The assignment is then based on the calculation of partial Raman activities arising from self- and/or cross-correlations between different types of Wannier functions in the system. Different types of Wannier functions can be distinguished based on their spatial spread. To demonstrate the predictive power of this approach, we applied it to the case of a cyclohexane molecule in the gas phase and were able to fully assign the simulated Raman peaks.

Adiabatic Molecular Orbital Tracking in Ab Initio Molecular Dynamics

The ab initio molecular dynamics (AIMD) method provides a computational route for the real-time simulation of reactive chemistry. An often-overlooked capability of this approach is the opportunity to examine the electronic evolution of a chemical system. For AIMD trajectories based on Hartree-Fock or density functional theory methods, the real-time evolution of single-particle molecular orbitals (MOs) can provide detailed insights into the time-dependent electronic structure of molecules. The evolving electronic Hamiltonians at each MD step pose problems for tracking and visualizing a given MO's character, ordering, and associated phase throughout an MD trajectory, however. This report presents and assesses a simple algorithm for correcting these deficiencies by exploiting similarity projections of the electronic structure between neighboring MD steps. Two aspects bring this analysis beyond a simple step-to-step projection scheme. First, the challenging case of coincidental orbital degeneracies is resolved via a quadrupole-field perturbation that nonetheless rigorously preserves energy conservation. Second, the resulting orbitals are shown to evolve adiabatically, in spite of the "preservation of character" concept that undergirds a projection of neighboring steps' MOs. The method is tested on water clusters, which exhibit considerable dynamic degeneracies, as well as a classic organic nucleophilic substitution reaction, in which the adiabatic evolution of the bonding orbitals clarifies textbook interpretations of the electronic structure during this reactive collision.

Cooperativity and Anticooperativity in Ion-Water Interactions: Implications for the Aqueous Solvation of Ions

Non-additive effects in hydrogen bonds (HB) take place as a consequence of electronic charge transfers. Therefore, it is natural to expect cooperativity and anticooperativity in ion-water interactions. Nevertheless, investigations on this matter are scarce. This paper addresses the interactions of (i) the cations Li⁺ , Na⁺ , K⁺ , Be²⁺ , Mg²⁺ , and Ca²⁺ together with (ii) the anions F^- , Cl^- , Br^- , NO₃ ^- and SO₄ ^2- with water clusters (H₂ O)_n , n=1-8, and the effects of these ions on the HBs within the complete molecular adducts. We used quantum chemical topology tools, specifically the quantum theory of atoms in molecules and the interacting quantum atoms energy partition to investigate non-additive effects among the interactions studied herein. Our results show a decrease on the interaction energy between ions and the first neighbouring water molecules with an increment of the coordination number. We also found strong cooperative effects in the interplay between HBs and ion-dipole interactions within the studied systems. Such cooperativity affects considerably the interactions among ions with their first and second solvation shells in aqueous environments. Overall, we believe this article provides valuable information about how ion-dipole contacts interact with each other and how they relate to other interactions, such as HBs, in the framework of non-additive effects in aqueous media.

Modeling LiF and FLiBe Molten Salts with Robust Neural Network Interatomic Potential

Lithium-based molten salts have attracted significant attention due to their applications in energy storage, advanced fission reactors, and fusion devices. Lithium fluorides and particularly 66.6%LiF-33.3%BeF₂ (Flibe) are of considerable interest in nuclear systems, as they show an excellent combination of favorable heat transfer, neutron moderation, and transmutation characteristics. For nuclear salts, the range of possible local structures, compositions, and thermodynamic conditions presents significant challenges in atomistic modeling. In this work, we demonstrate that atom-centered neural network interatomic potentials (NNIPs) provide a fast method for performing molecular dynamics of molten salts that is as accurate as ab initio molecular dynamics. For LiF, these potentials are able to accurately reproduce ab initio interactions of dimers, crystalline solids under deformation, crystalline LiF near the melting point, and liquid LiF at high temperatures. For Flibe, NNIPs accurately predict the structures and dynamics at normal operating conditions, high-temperature-pressure conditions, and in the crystalline solid phase. Furthermore, we show that NNIP-based molecular dynamics of molten salts are scalable to reach long time scales (e.g., nanosecond) and large system sizes (e.g., 10⁵ atoms) while maintaining ab initio density functional theory accuracy and providing more than 3 orders of magnitude of computational speedup for calculating structure and transport properties.

Progress and challenges in ab initio simulations of quantum nuclei in weakly bonded systems

Atomistic simulations based on the first-principles of quantum mechanics are reaching unprecedented length scales. This progress is due to the growth in computational power allied with the development of new methodologies that allow the treatment of electrons and nuclei as quantum particles. In the realm of materials science, where the quest for desirable emergent properties relies increasingly on soft weakly bonded materials, such methods have become indispensable. In this Perspective, an overview of simulation methods that are applicable for large system sizes and that can capture the quantum nature of electrons and nuclei in the adiabatic approximation is given. In addition, the remaining challenges are discussed, especially regarding the inclusion of nuclear quantum effects (NQEs) beyond a harmonic or perturbative treatment, the impact of NQEs on electronic properties of weakly bonded systems, and how different first-principles potential energy surfaces can change the impact of NQEs on the atomic structure and dynamics of weakly bonded systems.

A Review on Combination of Ab Initio Molecular Dynamics and NMR Parameters Calculations

This review focuses on a combination of ab initio molecular dynamics (aiMD) and NMR parameters calculations using quantum mechanical methods. The advantages of such an approach in comparison to the commonly applied computations for the structures optimized at 0 K are presented. This article was designed as a convenient overview of the applied parameters such as the aiMD type, DFT functional, time step, or total simulation time, as well as examples of previously studied systems. From the analysis of the published works describing the applications of such combinations, it was concluded that including fast, small-amplitude motions through aiMD has a noticeable effect on the accuracy of NMR parameters calculations.

Quantum Chemical Microsolvation by Automated Water Placement

We developed a quantitative approach to quantum chemical microsolvation. Key in our methodology is the automatic placement of individual solvent molecules based on the free energy solvation thermodynamics derived from molecular dynamics (MD) simulations and grid inhomogeneous solvation theory (GIST). This protocol enabled us to rigorously define the number, position, and orientation of individual solvent molecules and to determine their interaction with the solute based on physical quantities. The generated solute-solvent clusters served as an input for subsequent quantum chemical investigations. We showcased the applicability, scope, and limitations of this computational approach for a number of small molecules, including urea, 2-aminobenzothiazole, (+)-syn-benzotriborneol, benzoic acid, and helicene. Our results show excellent agreement with the available ab initio molecular dynamics data and experimental results.

The Many-Body Expansion for Aqueous Systems Revisited: II. Alkali Metal and Halide Ion-Water Interactions

We present a detailed study of the many-body expansion (MBE) for alkali metal and halide ion-water interactions and quantify the effect of these ions on the strength of the surrounding aqueous hydrogen bonding environment. Building on our previous work on neutral water clusters [J. P. Heindel and S. S. Xantheas, J. Chem. Theor. Comput. 16 (11), 6843-6855 (2020)], we carry out the MBE for the alkali metal and halide ion-water clusters, Z^+/-(H₂O)₉, where Z = Li⁺, Na⁺, K⁺, Rb⁺, Cs⁺, F^-, Cl^-, Br^-, and I^- and compare them with the results for a pure water cluster with the same number of "bodies", viz., (H₂O)₁₀. The 2-B ion-water (I-W) interaction accounts for a larger percentage of the total cluster binding energy compared to a pure water cluster of the same size, with the total 3-B term being smaller and of opposite sign (repulsive), whereas higher order terms are essentially negligible. The same oscillating behavior around zero for the MBE terms higher than the 5-B with a basis set that was reported for water clusters is also observed for the ion-water clusters considered here, with the basis set superposition error (BSSE) corrections amending this as in the water cluster case. A remarkable, linear anticorrelation between the total 2-B (I-W), the total 2-B (W-W), and also the 3-B (W-W-W) interactions is found, quantifying the effect of the different ions in disrupting and altering (weakening) the neighboring hydrogen bonded water network: stronger (I-W) interactions result in weaker (W-W) interactions. Additional linear correlations across the two series of alkali metals and halide ions were found between the 3-B (I-W-W) and the 2-B (I-W) as well as between the 3-B (I-W-W) and the 3-B (W-W-W) interactions, suggesting the existence of previously unrealized underlying physics governing these 2-B intermolecular and 3-B collective interactions. Our results further suggest a universal behavior of the two different families of ions (alkali metals and halides) for both the correlations of the various components of the total binding energies and the estimate of the 2-B BSSE correction, which is reported to follow a common profile for ion-water and water-water interactions when cast in terms of reduced distances and energies of the respective dimers. We expect the current results that quantify the interplay between ion-water and water-water interactions in aqueous clusters to impact the development of classical, ab initio-based force fields for monatomic ion solvation, whereas the insights into the nature of the BSSE to be critical in future ab initio-based, many-body molecular dynamics studies.

Neutral gold clusters studied by the isothermal Brownian-type molecular dynamics and metadynamics molecular dynamics simulations

The DFTB theory was combined with the isothermal Brownian-type molecular dynamics (MD) and metadynamics molecular dynamics (MMD) algorithms to perform simulation studies for Au clusters. Two representative DFTB parametrizations were investigated. In one parametrization, the DFTB-A, the Slater-Koster parameters in the DFTB energy function were determined focusing on the ionic repulsive energy part, E_rep and the other, the DFTB-B, due attention was paid to the electronic band-structure energy part, E_band . Minimized structures of these two parametrizations were separately applied in MD and MMD simulations to generate unbiased and biased trajectories in collective variable (CV) space, respectively. Here, we found the MD simulations monitored at 300 K manifest fluxional characteristics in planar cluster Au₉ /DFTB-A, but give no discernible tracts of fluxionality for planar Au₈ /DFTB-A and Au₈ /DFTB-B, for nonplanar Au₁₀ /DFTB-A and, to some extent, for nonplanar Au₉ /DFTB-B; they are plausibly being hindered by higher-than k_B T energy barriers. Very recent FIR-MPD spectroscopy measurements, however, were reported to have detected at 300 K both the planar and nonplanar neutral Au_n clusters in the size range 5 ≤ n ≤ 13. The failure of MD simulations has prompted us to apply the MMD simulation and construct the free energy landscape (FEL) in CV space. Through scrutinizing the FELs of these clusters and their associated structures, we examine the relative importance of E_rep /DFTB-A and E_band /DFTB-B in unraveling the covalent-like behavior of valence electrons in Au_n . Most important of all, we shall evaluate the DFTB parametrization in MMD strategy through comparing extensively the simulation data recorded with the gas-phase experimental data.

Hydrogen bond dynamics of interfacial water molecules revealed from two-dimensional vibrational sum-frequency generation spectroscopy

Vibrational sum-frequency generation (vSFG) spectroscopy allows the study of the structure and dynamics of interfacial systems. In the present work, we provide a simple recipe, based on a narrowband IR pump and broadband vSFG probe technique, to computationally obtain the two-dimensional vSFG spectrum of water molecules at the air-water interface. Using this technique, to study the time-dependent spectral evolution of hydrogen-bonded and free water molecules, we demonstrate that at the interface, the vibrational spectral dynamics of the free OH bond is faster than that of the bonded OH mode.

"On-The-Fly" Calculation of the Vibrational Sum-Frequency Generation Spectrum at the Air-Water Interface

In the present work, we provide an electronic structure based method for the "on-the-fly" determination of vibrational sum frequency generation (v-SFG) spectra. The predictive power of this scheme is demonstrated at the air-water interface. While the instantaneous fluctuations in dipole moment are obtained using the maximally localized Wannier functions, the fluctuations in polarizability are approximated to be proportional to the second moment of Wannier functions. The spectrum henceforth obtained captures the signatures of hydrogen bond stretching, bending, as well as low-frequency librational modes.

CP2K: An electronic structure and molecular dynamics software package - Quickstep: Efficient and accurate electronic structure calculations

CP2K is an open source electronic structure and molecular dynamics software package to perform atomistic simulations of solid-state, liquid, molecular, and biological systems. It is especially aimed at massively parallel and linear-scaling electronic structure methods and state-of-the-art ab initio molecular dynamics simulations. Excellent performance for electronic structure calculations is achieved using novel algorithms implemented for modern high-performance computing systems. This review revisits the main capabilities of CP2K to perform efficient and accurate electronic structure simulations. The emphasis is put on density functional theory and multiple post-Hartree-Fock methods using the Gaussian and plane wave approach and its augmented all-electron extension.

Water structure near the surface of Weyl semimetals as catalysts in photocatalytic proton reduction

In this work, second-generation Car-Parrinello-based mixed quantum-classical mechanics molecular dynamics simulations of small nanoparticles of NbP, NbAs, TaAs, and 1T-TaS₂ in water are presented. The first three materials are topological Weyl semimetals, which were recently discovered to be active catalysts in photocatalytic water splitting. The aim of this research was to correlate potential differences in the water structure in the vicinity of the nanoparticle surface with the photocatalytic activity of these materials in light induced proton reduction. The results presented herein allow explaining the catalytic activity of these Weyl semimetals: the most active material, NbP, exhibits a particularly low water coordination near the surface of the nanoparticle, whereas for 1T-TaS₂, with the lowest catalytic activity, the water structure at the surface is most ordered. In addition, the photocatalytic activity of several organic and metalorganic photosensitizers in the hydrogen evolution reaction was experimentally investigated with NbP as the proton reduction catalyst. Unexpectedly, the charge of the photosensitizer plays a decisive role for the photocatalytic performance.

Molecular Structure and Modeling of Water-Air and Ice-Air Interfaces Monitored by Sum-Frequency Generation

From a glass of water to glaciers in Antarctica, water-air and ice-air interfaces are abundant on Earth. Molecular-level structure and dynamics at these interfaces are key for understanding many chemical/physical/atmospheric processes including the slipperiness of ice surfaces, the surface tension of water, and evaporation/sublimation of water. Sum-frequency generation (SFG) spectroscopy is a powerful tool to probe the molecular-level structure of these interfaces because SFG can specifically probe the topmost interfacial water molecules separately from the bulk and is sensitive to molecular conformation. Nevertheless, experimental SFG has several limitations. For example, SFG cannot provide information on the depth of the interface and how the orientation of the molecules varies with distance from the surface. By combining the SFG spectroscopy with simulation techniques, one can directly compare the experimental data with the simulated SFG spectra, allowing us to unveil the molecular-level structure of water-air and ice-air interfaces. Here, we present an overview of the different simulation protocols available for SFG spectra calculations. We systematically compare the SFG spectra computed with different approaches, revealing the advantages and disadvantages of the different methods. Furthermore, we account for the findings through combined SFG experiments and simulations and provide future challenges for SFG experiments and simulations at different aqueous interfaces.

Modeling Dynamic Conformations of Organic Molecules: Alkyne Carotenoids in Solution

Calculating the spectroscopic properties of complex conjugated organic molecules in their relaxed state is far from simple. An additional complexity arises for flexible molecules in solution, where the rotational energy barriers are low enough so that nonminimum conformations may become dynamically populated. These metastable conformations quickly relax during the minimization procedures preliminary to density functional theory calculations, and so accounting for their contribution to the experimentally observed properties is problematic. We describe a strategy for stabilizing these nonminimum conformations in silico, allowing their properties to be calculated. Diadinoxanthin and alloxanthin present atypical vibrational properties in solution, indicating the presence of several conformations. Performing energy calculations in vacuo and polarizable continuum model calculations in different solvents, we found three different conformations with values for the δ dihedral angle of the end ring ca. 0, 180, and 90° with respect to the plane of the conjugated chain. The latter conformation, a nonglobal minimum, is not stable during the minimization necessary for modeling its spectroscopic properties. To circumvent this classical problem, we used a Car–Parinello MD supermolecular approach, in which diadinoxanthin was solvated by water molecules so that metastable conformations were stabilized by hydrogen-bonding interactions. We progressively removed the number of solvating waters to find the minimum required for this stabilization. This strategy represents the first modeling of a carotenoid in a distorted conformation and provides an accurate interpretation of the experimental data.

On the Possibility of Helium Adsorption in Nitrogen Doped Graphitic Materials

The potassium salt of polyheptazine imide (K-PHI) is a promising photocatalyst for various chemical reactions. From powder X-ray diffraction data an idealized structural model of K-PHI has been derived. Using atomic coordinates of this model we defined an energetically optimized K-PHI structure, in which the K ions are present in the pore and between the PHI-planes. The distance between the anion framework and K+ resembles a frustrated Lewis pair-like structure, which we denote as frustrated Coulomb pair that results in an interesting adsorption environment for otherwise non-adsorbing, non-polar gas molecules. We demonstrate that even helium (He) gas molecules, which are known to have the lowest boiling point and the lowest intermolecular interactions, can be adsorbed in this polarized environment with an adsorption energy of  - 4.6 kJ mol-1 per He atom. The interaction between He atoms and K-PHI is partially originating from charge transfer, as disclosed by our energy decomposition analysis based on absolutely localized molecular orbitals. Due to very small charge transfer interactions, He gas adsorption saturates at 8 at%, which however can be subject to further improvement by cation variation.

Energy transfer within the hydrogen bonding network of water following resonant terahertz excitation

Energy dissipation in water is very fast and more efficient than in many other liquids. This behavior is commonly attributed to the intermolecular interactions associated with hydrogen bonding. Here, we investigate the dynamic energy flow in the hydrogen bond network of liquid water by a pump-probe experiment. We resonantly excite intermolecular degrees of freedom with ultrashort single-cycle terahertz pulses and monitor its Raman response. By using ultrathin sample cell windows, a background-free bipolar signal whose tail relaxes monoexponentially is obtained. The relaxation is attributed to the molecular translational motions, using complementary experiments, force field, and ab initio molecular dynamics simulations. They reveal an initial coupling of the terahertz electric field to the molecular rotational degrees of freedom whose energy is rapidly transferred, within the excitation pulse duration, to the restricted translational motion of neighboring molecules. This rapid energy transfer may be rationalized by the strong anharmonicity of the intermolecular interactions.

Synergistic Approach of Ultrafast Spectroscopy and Molecular Simulations in the Characterization of Intramolecular Charge Transfer in Push-Pull Molecules

The comprehensive characterization of Intramolecular Charge Transfer (ICT) stemming in push-pull molecules with a delocalized π-system of electrons is noteworthy for a bespoke design of organic materials, spanning widespread applications from photovoltaics to nanomedicine imaging devices. Photo-induced ICT is characterized by structural reorganizations, which allows the molecule to adapt to the new electronic density distribution. Herein, we discuss recent photophysical advances combined with recent progresses in the computational chemistry of photoactive molecular ensembles. We focus the discussion on femtosecond Transient Absorption Spectroscopy (TAS) enabling us to follow the transition from a Locally Excited (LE) state to the ICT and to understand how the environment polarity influences radiative and non-radiative decay mechanisms. In many cases, the charge transfer transition is accompanied by structural rearrangements, such as the twisting or molecule planarization. The possibility of an accurate prediction of the charge-transfer occurring in complex molecules and molecular materials represents an enormous advantage in guiding new molecular and materials design. We briefly report on recent advances in ultrafast multidimensional spectroscopy, in particular, Two-Dimensional Electronic Spectroscopy (2DES), in unraveling the ICT nature of push-pull molecular systems. A theoretical description at the atomistic level of photo-induced molecular transitions can predict with reasonable accuracy the properties of photoactive molecules. In this framework, the review includes a discussion on the advances from simulation and modeling, which have provided, over the years, significant information on photoexcitation, emission, charge-transport, and decay pathways. Density Functional Theory (DFT) coupled with the Time-Dependent (TD) framework can describe electronic properties and dynamics for a limited system size. More recently, Machine Learning (ML) or deep learning approaches, as well as free-energy simulations containing excited state potentials, can speed up the calculations with transferable accuracy to more complex molecules with extended system size. A perspective on combining ultrafast spectroscopy with molecular simulations is foreseen for optimizing the design of photoactive compounds with tunable properties.

Accessing the Accuracy of Density Functional Theory through Structure and Dynamics of the Water-Air Interface

Density functional theory-based molecular dynamics simulations are increasingly being used for simulating aqueous interfaces. Nonetheless, the choice of the appropriate density functional, critically affecting the outcome of the simulation, has remained arbitrary. Here, we assess the performance of various exchange-correlation (XC) functionals, based on the metrics relevant to sum-frequency generation spectroscopy. The structure and dynamics of water at the water-air interface are governed by heterogeneous intermolecular interactions, thereby providing a critical benchmark for XC functionals. We find that the XC functionals constrained by exact functional conditions (revPBE and revPBE0) with the dispersion correction show excellent performance. The poor performance of the empirically optimized density functional (M06-L) indicates the importance of satisfying the exact functional condition. Understanding the performance of different XC functionals can aid in resolving the controversial interpretation of the interfacial water structure and direct the design of novel, improved XC functionals better suited to describing the heterogeneous interactions in condensed phases.

Car-Parrinello Monitor for More Robust Born-Oppenheimer Molecular Dynamics

Born-Oppenheimer molecular dynamics (BOMD) is a promising simulation method for exploring the possible reaction pathways of a chemical system, but one significant challenge is the increased difficulty of converging the self-consistent field (SCF) calculation that often accompanies the breaking and forming of chemical bonds. To address this challenge, we developed an enhancement to the BOMD simulation method called the Car-Parrinello monitor (CPMonitor) that uses Car-Parrinello molecular dynamics (CPMD) to recover from SCF convergence failures. CPMonitor works by detecting SCF convergence failures in BOMD and switching to a CPMD Hamiltonian to propagate through the region of configuration space where the SCF calculation is unable to converge, then switching back to BOMD when good convergence behavior is re-established. We present a series of simulation studies that use CPMonitor, including detailed studies of the thermodynamic and dynamical properties of simple systems, as well as ab initio nanoreactor simulations containing transition metal atoms that were previously not possible to simulate using standard BOMD methods. Our studies show that CPMonitor can make BOMD simulations robust to SCF convergence difficulties and improve simulation performance and stability in reaction discovery applications.

Enhancement of the local asymmetry in the hydrogen bond network of liquid water by an ultrafast electric field pulse

Condensed phase electron decomposition analysis based on density functional theory has recently revealed an asymmetry in the hydrogen-bond network in liquid water, in the sense that a significant population of water molecules are simultaneously donating and accepting one strong hydrogen-bond and another substantially weaker one. Here we investigate this asymmetry, as well as broader structural and energetic features of water's hydrogen-bond network, following the application of an intense electric field square pulse that invokes the ultrafast reorientation of water molecules. We find that the necessary field-strength required to invoke an ultrafast alignment in a picosecond time window is on the order of 108 Vm-1. The resulting orientational anisotropy imposes an experimentally measurable signature on the structure and dynamics of the hydrogen-bond network, including its asymmetry, which is strongly enhanced. The dependence of the molecular reorientation dynamics on the field-strength can be understood by relating the magnitude of the water dipole-field interaction to the rotational kinetic energy, as well as the hydrogen-bond energy.

Atmospheric Spectroscopy and Photochemistry at Environmental Water Interfaces

The air-water interface is ubiquitous in nature, as manifested in the form of the surfaces of oceans, lakes, and atmospheric aerosols. The aerosol interface, in particular, can play a crucial role in atmospheric chemistry. The adsorption of atmospheric species onto and into aerosols modifies their concentrations and chemistries. Moreover, the aerosol phase allows otherwise unlikely solution-phase chemistry to occur in the atmosphere. The effect of the air-water interface on these processes is not entirely known. This review summarizes recent theoretical investigations of the interactions of atmosphere species with the air-water interface, including reactant adsorption, photochemistry, and the spectroscopy of reactants at the water surface, with an emphasis on understanding differences between interfacial chemistries and the chemistries in both bulk solution and the gas phase. The results discussed here enable an understanding of fundamental concepts that lead to potential air-water interface effects, providing a framework to understand the effects of water surfaces on our atmosphere.

Computational Approach to Molecular Catalysis by 3d Transition Metals: Challenges and Opportunities

Computational chemistry provides a versatile toolbox for studying mechanistic details of catalytic reactions and holds promise to deliver practical strategies to enable the rational in silico catalyst design. The versatile reactivity and nontrivial electronic structure effects, common for systems based on 3d transition metals, introduce additional complexity that may represent a particular challenge to the standard computational strategies. In this review, we discuss the challenges and capabilities of modern electronic structure methods for studying the reaction mechanisms promoted by 3d transition metal molecular catalysts. Particular focus will be placed on the ways of addressing the multiconfigurational problem in electronic structure calculations and the role of expert bias in the practical utilization of the available methods. The development of density functionals designed to address transition metals is also discussed. Special emphasis is placed on the methods that account for solvation effects and the multicomponent nature of practical catalytic systems. This is followed by an overview of recent computational studies addressing the mechanistic complexity of catalytic processes by molecular catalysts based on 3d metals. Cases that involve noninnocent ligands, multicomponent reaction systems, metal-ligand and metal-metal cooperativity, as well as modeling complex catalytic systems such as metal-organic frameworks are presented. Conventionally, computational studies on catalytic mechanisms are heavily dependent on the chemical intuition and expert input of the researcher. Recent developments in advanced automated methods for reaction path analysis hold promise for eliminating such human-bias from computational catalysis studies. A brief overview of these approaches is presented in the final section of the review. The paper is closed with general concluding remarks.

Recent Advances in Metallic Glass Nanostructures: Synthesis Strategies and Electrocatalytic Applications

Recent advances in metallic glass nanostructures (MGNs) are reported, covering a wide array of synthesis strategies, computational discovery, and design solutions that provide insight into distinct electrocatalytic applications. A brief introduction to the development and unique features of MGNs with an overview of top-down and bottom-up synthesis strategies is presented. Specifically, the morphology and structural analysis of several examples applying MGNs as electrodes are highlighted. Subsequently, a comprehensive discussion of commonly employed kinetic parameters and their connection with the unique material structures of MGNs on individual electrocatalytic reactions is made, including the hydrogen evolution reaction, oxygen reduction reaction, and alcohol (methanol or ethanol) oxidation reaction. Finally, a summary of the challenges and perspective on the future research and development relevant to MGNs as electrocatalysts is provided.

On the Hydrogen Bond Strength and Vibrational Spectroscopy of Liquid Water

In the present work, we introduce two new metrics i.e. hydrogen-bond strength and charge-transfer between the donor/acceptor water molecules as a measure of hydrogen-bond rearrangement dynamics. Further, we also employ a simple model based on energy flux through the donor-acceptor water pairs to quantify the extent of the local hydrogen-bond network reorganization. Most importantly, we report a linear relationship between the OH stretch frequency and the charge and energy transfer through donor-acceptor water pairs. We demonstrate that the vibrational frequency fluctuations, which are used to determine third-order non-linear spectroscopic observables like the short-time slope of three pulse photon echo, can be used as an analog of the fluctuations in the hydrogen-bond strength and charge-transfer. The timescales obtained from our hydrogen-bond strength correlation and charge-transfer correlation decay are in excellent agreement with the computed frequency-time correlation function, as well as with recent vibrational echo experiments.

Impact of finite-temperature and condensed-phase effects on theoretical X-ray absorption spectra of transition metal complexes

The impact of condensed-phase and finite-temperature effects on the theoretical X-ray absorption spectra of transition metal complexes is assessed. The former are included in terms of the all-electron Gaussian and augmented plane-wave approach, whereas the latter are taken into account by extensive ensemble averaging along second-generation Car-Parrinello ab initio molecular dynamics trajectories. We find that employing the periodic boundary conditions and including finite-temperature effects systematically improves the agreement between our simulated X-ray absorption spectra and experimental measurements. © 2018 Wiley Periodicals, Inc.