Unraveling the synergism mechanistic insight of O-vacancy and interfacial charge transfer in WO3-x decorated on Ag2CO3/BiOBr for photocatalysis of water pollutants: Based on experimental and density functional theory (DFT) studies

Photocatalysis has been widely used as one of the most promising approaches to remove various pollutants in liquid or gas phases during the last decade. The main emphasis of the study is on the synergy of vacancy engineering and heterojunction formation, two widely used modifying approaches, to significantly alter photocatalytic performance. The vacancy-induced Ag2CO3/BiOBr/WO3-x heterojunction system has been fabricated using a co-precipitation technique to efficiently abate methylene blue (MB) dye and doxycycline (DC) antibiotic. The as-fabricated Ag2CO3/BiOBr/WO3-x heterojunction system displayed improved optoelectronic characteristic features because of the rational combination of dual charge transferal route and defect modulation. The Ag2CO3/BiOBr/WO3-x system possessed 97% and 74% photodegradation efficacy for MB and DC, respectively, with better charge isolation and migration efficacy. The ternary photocatalyst possessed a multi-fold increase in the reaction rate for both MB and DC, i.e., 0.021 and 0.0078 min-1, respectively, compared to pristine counterparts. Additionally, more insightful deductions about the photodegradation routes were made possible by the structural investigations of MB and DC using density functional theory (DFT) simulations. This study advances the understanding of the mechanisms forming visible light active dual Z-scheme heterojunction for effective environmental remediation.

The variational quantum eigensolver self-consistent field method within a polarizable embedded framework

We formulate and implement the Variational Quantum Eigensolver Self Consistent Field (VQE-SCF) algorithm in combination with polarizable embedding (PE), thereby extending PE to the regime of quantum computing. We test the resulting algorithm, PE-VQE-SCF, on quantum simulators and demonstrate that the computational stress on the quantum device is only slightly increased in terms of gate counts compared to regular VQE-SCF. On the other hand, no increase in shot noise was observed. We illustrate how PE-VQE-SCF may lead to the modeling of real chemical systems using a simulation of the reaction barrier of the Diels-Alder reaction between furan and ethene as an example.

Adsorption of HCN on cosmic silicates: a periodic quantum mechanical study

Hydrogen cyanide (HCN) represents a small but widely distributed fraction of the interstellar molecules, and it has been observed in all the environments characterizing the formation of a new planetary system. HCN can polymerize to form biomolecules, including adenine (H5C5N5), and it has drawn attention as a possible precursor of several building blocks of life due to the presence of its polymerization products in meteorites, comets and other asteroidal bodies. To elucidate the potential catalytic role that cosmic silicates have played in these processes, we have investigated, at DFT-PBE level inclusive of a posteriori dispersion correction, the energetic and spectroscopic features of the adsorption of HCN molecules on the most relevant crystalline surfaces of the mineral forsterite (Mg2SiO4), a common silicate constituent of the interstellar core grains and planetary rocky bodies. The results reveal that HCN adsorbs both in molecular and dissociative ways, within a wide range of adsorption energies (-29.4 to -466.4 kJ mol-1). Thermodynamic and kinetic results show that dissociative adsorption is dominant already at low temperatures, a fact particularly relevant at the protoplanetary conditions (i.e., the latest stages in the star system formation process). The simulated spectroscopic features of the studied adducts show a wide range of different degrees of perturbation of C-H and CN bonds. This finding agrees with previous experimental works, and our results confirm that a complex chemistry is observed when this astrochemically-relevant molecule interacts with Mg2SiO4, which may be associated with a considerable potential reactivity towards the formation of relevant prebiotic compounds.

Energy and Electronic Properties of Nanostructures Based on the CL-20 Framework with the Replacement of the Carbon Atoms by Silicon and Germanium: A Density Functional Theory Study

We consider Si_nCL-20 and Ge_nCL-20 systems with carbon atoms replaced by silicon/germanium atoms and their dimers. The physicochemical properties of the silicon/germanium analogs of the high-energy molecule CL-20 and its dimers were determined and studied using density functional theory with the B3LYP/6-311G(d,p) level of theory. It was found that the structure and geometry of Si_nCL-20/Ge_nCL-20 molecules change dramatically with the appearance of Si-/Ge-atoms. The main difference between silicon- or germanium-substituted Si_nCL-20/Ge_nCL-20 molecules and the pure CL-20 molecule is that the NO₂ functional groups make a significant rotation relative to the starting position in the classical molecule, and the effective diameter of the frame of the systems increases with the addition of Si-/Ge-atoms. Thus, the effective framework diameter of a pure CL-20 molecule is 3.208 Å, while the effective diameter of a fully silicon-substituted Si₆CL-20 molecule is 4.125 Å, and this parameter for a fully germanium-substituted Ge₆CL-20 molecule is 4.357 Å. The addition of silicon/germanium atoms to the system leads to a decrease in the binding energy. In detail, the binding energies for CL-20/Si₆CL-20/Ge₆CL-20 molecules are 4.026, 3.699, 3.426 eV/atom, respectively. However, it has been established that the framework maintains stability, with an increase in the number of substituting silicon or germanium atoms. In addition, we designed homodesmotic reactions for the CL-20 molecule and its substituted derivatives Si₆CL-20/Ge₆CL-20, and then determined the strain energy to find out in which case more energy would be released when the framework breaks. Further, we also studied the electronic properties of systems based on CL-20 molecules. It was found that the addition of germanium or silicon atoms instead of carbon leads to a decrease in the size of the HOMO-LUMO gap. Thus, the HOMO-LUMO gaps of the CL-20/Si₆CL-20/Ge₆CL-20 molecules are 5.693, 5.339, and 5.427 eV, respectively. A similar dependence is also observed for CL-20 dimers. So, in this work, we have described in detail the dependence of the physicochemical parameters of CL-20 molecules and their dimers on the types of atoms upon substitution.

Water formation on interstellar silicates: the role of Fe2+/H2 interactions in the O + H2 → H2O reaction

Water formation by reaction of H2 and O on silicate surfaces as a first step towards the generation of interstellar ice mantles is possible thanks to the activation of H2 inferred by Fe2+ ions and quantum tunnelling effects.

The E2 state of FeMoco: Hydride Formation versus Fe Reduction and a Mechanism for H2 Evolution

The iron‐molybdenum cofactor (FeMoco) is responsible for dinitrogen reduction in Mo nitrogenase. Unlike the resting state, E₀, reduced states of FeMoco are much less well characterized. The E₂ state has been proposed to contain a hydride but direct spectroscopic evidence is still lacking. The E₂ state can, however, relax back the E₀ state via a H₂ side‐reaction, implying a hydride intermediate prior to H₂ formation. This E₂→E₀ pathway is one of the primary mechanisms for H₂ formation under low‐electron flux conditions. In this study we present an exploration of the energy surface of the E₂ state. Utilizing both cluster‐continuum and QM/MM calculations, we explore various classes of E₂ models: including terminal hydrides, bridging hydrides with a closed or open sulfide‐bridge, as well as models without. Importantly, we find the hemilability of a protonated belt‐sulfide to strongly influence the stability of hydrides. Surprisingly, non‐hydride models are found to be almost equally favorable as hydride models. While the cluster‐continuum calculations suggest multiple possibilities, QM/MM suggests only two models as contenders for the E₂ state. These models feature either i) a bridging hydride between Fe₂ and Fe₆ and an open sulfide‐bridge with terminal SH on Fe₆ (E₂‐hyd) or ii) a double belt‐sulfide protonated, reduced cofactor without a hydride (E₂‐nonhyd). We suggest both models as contenders for the E₂ redox state and further calculate a mechanism for H₂ evolution. The changes in electronic structure of FeMoco during the proposed redox‐state cycle, E₀→E₁→E₂→E₀, are discussed.

Nudged Elastic Band Method for Molecular Reactions Using Energy-Weighted Springs Combined with Eigenvector Following

The climbing image nudged elastic band method (CI-NEB) is used to identify reaction coordinates and to find saddle points representing transition states of reactions. It can make efficient use of parallel computing as the calculations of the discretization points, the so-called images, can be carried out simultaneously. In typical implementations, the images are distributed evenly along the path by connecting adjacent images with equally stiff springs. However, for systems with a high degree of flexibility, this can lead to poor resolution near the saddle point. By making the spring constants increase with energy, the resolution near the saddle point is improved. To assess the performance of this energy-weighted CI-NEB method, calculations are carried out for a benchmark set of 121 molecular reactions. The performance of the method is analyzed with respect to the input parameters. Energy-weighted springs are found to greatly improve performance and result in successful location of the saddle points in less than a thousand energy and force evaluations on average (about a hundred per image) using the same set of parameter values for all of the reactions. Even better performance is obtained by stopping the calculation before full convergence and complete the saddle point search using an eigenvector following method starting from the location of the climbing image. This combination of methods, referred to as NEB-TS, turns out to be robust and highly efficient as it reduces the average number of energy and force evaluations down to a third, to 305. An efficient and flexible implementation of these methods has been made available in the ORCA software.

All-Nitrogen Cages and Molecular Crystals: Topological Rules, Stability, and Pyrolysis Paths

We combined ab initio molecular dynamics with the intrinsic reaction coordinate in order to investigate the mechanisms of stability and pyrolysis of N4 ÷ N120 fullerene-like nitrogen cages. The stability of the cages was evaluated in terms of the activation barriers and the activation Gibbs energies of their thermal-induced breaking. We found that binding energies, bond lengths, and quantum-mechanical descriptors failed to predict the stability of the cages. However, we derived a simple topological rule that adjacent hexagons on the cage surface resulted in its instability. For this reason, the number of stable nitrogen cages is significantly restricted in comparison with their carbon counterparts. As a rule, smaller clusters are more stable, whereas the earlier proposed large cages collapse at room temperature. The most stable all-nitrogen cages are the N4 and N6 clusters, which can form the van der Waals crystals with densities of 1.23 and 1.36 g/cm3, respectively. The examination of their band structures and densities of electronic states shows that they are both insulators. Their power and sensitivity are not inferior to the modern advanced high-energy nanosystems.

The reaction mechanism of Zika virus NS2B/NS3 serine protease inhibition by dipeptidyl aldehyde: a QM/MM study

Zika virus (ZIKV) infection has become a global public health problem, associated with microcephaly in newborns and Guillain-Barré syndrome in adults. Currently, there are no commercially available anti-ZIKV drugs. The viral protease NS2B/NS3, which is involved in viral replication and maturation, is a potential drug target. Peptidomimetic aldehyde inhibitors bind covalently to the catalytic S135 of the NS3 protease. Here, we apply hybrid quantum mechanics/molecular mechanics (QM/MM) free-energy simulations at the PDDG-PM3/ff14SB level to investigate the inhibition mechanism of the ZIKV protease by a dipeptidyl aldehyde inhibitor (acyl-KR-aldehyde). The results show that proton transfer from the catalytic S135 to H51 occurs in concert with nucleophilic addition on the aldehyde warhead by S135. The anionic covalent complex between the dipeptidyl aldehyde and the ZIKV protease is analogous to the tetrahedral intermediate for substrate hydrolysis. Spontaneous protonation by H51 forms the hemiacetal. In addition, we use correlated ab initio QM/MM potential energy path calculations at levels up to LCCSD(T)/(aug)-cc-pVTZ to obtain accurate potential energy profiles of the reaction, which also support a concerted mechanism. These results provide detailed insight into the mechanism of ZIKV protease inhibition by a peptidyl aldehyde inhibitor, which will guide in the design of inhibitors.

Nitro Derivatives of Silaprismanes as High-Energy Compounds: Theoretical Study

We present ab initio study of structures and properties of silaprismanes Si[Formula: see text]H[Formula: see text] and their nitro derivatives Si[Formula: see text]H[Formula: see text]NO2 ([Formula: see text]–10). We found that silaprismane Si[Formula: see text]H9NO2 possesses the highest stability among all studied cages. Attached NO2 group results in weak decreasing of the HOMO–LUMO gap. The smaller prismanes bind with NO2 groups more strongly. The comparison between silaprismanes and carbon prismanes is also discussed.

Silicate-mediated interstellar water formation: A theoretical study

Water is one of the most abundant molecules in the form of solid ice phase in the different regions of the interstellar medium (ISM). This large abundance cannot be properly explained by using only traditional low temperature gas-phase reactions. Thus, surface chemical reactions are believed to be major synthetic channels for the formation of interstellar water ice. Among the different proposals, hydrogenation of atomic O (i.e., 2H + O → H2O) is a chemically "simple" and plausible reaction toward water formation occurring on the surfaces of interstellar grains. Here, novel theoretical results concerning the formation of water adopting this mechanism on the crystalline (010) Mg2SiO4 surface (a unequivocally identified interstellar silicate) are presented. The investigated reaction aims to simulate the formation of the first water ice layer covering the silicate core of dust grains. Adsorption of the atomic O as a first step of the reaction has been computed, results indicating that a peroxo (O 2 2 - ) group is formed. The following steps involve the adsorption, diffusion and reaction of two successive H atoms with the adsorbed O atom. Results indicate that H diffusion on the surface has barriers of 4-6 kcal mol-1, while actual formation of OH and H2O present energy barriers of 22-23 kcal mol-1. Kinetic study results show that tunneling is crucial for the occurrence of the reactions and that formation of OH and H2O are the bottlenecks of the overall process. Several astrophysical implications derived from the theoretical results are provided as concluding remarks.

Electronic and reactivity characteristics of CL-20 covalent chains and networks: a density functional theory study

The low-dimensional nanostructures built with CL-20 units were studied in order to investigate the possibility of creating covalent CL-20 crystals.

Enantioselective Anion Recognition by Chiral Halogen-Bonding [2]Rotaxanes

The application of chiral interlocked host molecules for discrimination of guest enantiomers has been largely overlooked, which is surprising given their unique three-dimensional binding cavities capable of guest encapsulation. Herein, we combined the stringent linear geometric interaction constraints of halogen bonding (XB), the noncovalent interaction between an electrophilic halogen atom and a Lewis base, with highly preorganized and conformationally restricted chiral cavities of [2]rotaxanes to achieve enantioselective anion recognition. Representing the first detailed investigation of the use of chiral XB rotaxanes for this purpose, extensive ¹H NMR binding studies and molecular dynamics (MD) simulation experiments revealed that the chiral rotaxane cavity significantly enhances enantiodiscrimination compared to the non-interlocked free axle and macrocycle components. Furthermore, by examining the enantioselectivities of a family of structurally similar XB [2]rotaxanes containing different combinations of chiral and achiral macrocycle and axle components, the dominant influence of the chiral macrocycle in our rotaxane design for determining the effectiveness of chiral discrimination is demonstrated. MD simulations reveal the crucial geometric roles played by the XB interactions in orientating the bound enantiomeric anion guests for chiral selectivity, as well as the critical importance of the anions' hydration shells in governing binding affinity and enantiodiscrimination.

Proton Shuttling and Reaction Paths in Dissociative Electron Attachment to o- and p-Tetrafluorohydroquinone, an Experimental and Theoretical Study

Here we present a combined experimental and theoretical study on the fragmentation of o- and p-tetrafluorohydroquinone upon low energy electron attachment. Despite an identical ring-skeleton and identical functional groups in these constitutional isomers, they show distinctly different fragmentation patterns, a phenomenon that cannot be explained by distinct resonances or different thermochemistry. Using high-level quantum chemical calculations with the computationally affordable domain localized pair natural orbital approach, DLPNO-CCSD(T), we are able to provide a complete and accurate description of the respective reaction dynamics, revealing proton shuttling and transition states for competing channels as the explanation for the different behavior of these isomers. The results represent a "schoolbook example" of how the combination of experiment and modern high-level theory may today provide a thorough understanding of complex reaction dynamics by computationally affordable means.

Comparison of classical reaction paths and tunneling paths studied with the semiclassical instanton theory

Comparison of classical reaction paths and semiclassical instanton paths for a proton transfer reaction mechanism.

Trends in the adsorption and reactivity of hydrogen on magnesium silicate nanoclusters

We investigate the potential role of ultrasmall silicate grains in interstellar hydrogen-based chemistry by modelling H adsorption and H₂ formation/dissociation on nanosilicates.

Reaction mechanism of the bicopper enzyme peptidylglycine α-hydroxylating monooxygenase

Peptidylglycine α-hydroxylating monooxygenase is a noninteracting bicopper enzyme that stereospecifically hydroxylates the terminal glycine of small peptides for its later amidation. Neuroendocrine messengers, such as oxytocin, rely on the biological activity of this enzyme. Each catalytic turnover requires one oxygen molecule, two protons from the solvent, and two electrons. Despite this enzyme having been widely studied, a consensus on the reaction mechanism has not yet been found. Experiments and theoretical studies favor a pro-S abstraction of a hydrogen atom followed by the rebinding of an OH group. However, several hydrogen-abstracting species have been postulated; because two protons are consumed during the reaction, several protonation states are available. An electron transfer between the copper atoms could play a crucial role for the catalysis as well. This leads to six possible abstracting species. In this study, we compare them on equal footing. We perform quantum mechanics/molecular mechanics calculations, considering the glycine hydrogen abstraction. Our results suggest that the most likely mechanism is a protonation of the abstracting species before the hydrogen abstraction and another protonation as well as a reduction before OH rebinding.

B3LYP periodic study of the physicochemical properties of the nonpolar (010) Mg-pure and fe-containing olivine surfaces

B3LYP periodic simulations have been carried out to study some physicochemical properties of the bulk structures and the corresponding nonpolar (010) surfaces of Mg-pure and Fe-containing olivine systems; i.e., Mg2SiO4 (Fo) and Mg1.5Fe0.5SiO4 (Fo75). A detailed structural analysis of the (010) Fo and Fo75 surface models shows the presence of coordinatively unsaturated metal cations (Mg(2+) and Fe(2+), respectively) with shorter metal-O distances compared to the bulk ones. Energetic analysis devoted to the Fe(2+) electronic spin configuration and to the ion position in the surfaces reveals that Fe(2+) in its quintet state and placed at the outermost positions of the slab constitutes the most stable Fe-containing surface, which is related to the higher stability of high spin states when Fe(2+) is coordinatively unsaturated. Comparison of the simulated IR and the corresponding reflectance spectra indicates that Fe(2+) substitution induces an overall bathochromic shift of the spectra due to the larger mass of Fe compared to Mg cation. In contrast, the IR spectra of the surfaces are shifted to upper values and exhibit more bands compared to the corresponding bulk systems due to the shorter metal-O distances given in the coordinatively unsaturated metals and to symmetry reduction which brings nonequivalent motions between the outermost and the internal modes, respectively.

QM/MM investigations of organic chemistry oriented questions

About 35 years after its first suggestion, QM/MM became the standard theoretical approach to investigate enzymatic structures and processes. The success is due to the ability of QM/MM to provide an accurate atomistic picture of enzymes and related processes. This picture can even be turned into a movie if nuclei-dynamics is taken into account to describe enzymatic processes. In the field of organic chemistry, QM/MM methods are used to a much lesser extent although almost all relevant processes happen in condensed matter or are influenced by complicated interactions between substrate and catalyst. There is less importance for theoretical organic chemistry since the influence of nonpolar solvents is rather weak and the effect of polar solvents can often be accurately described by continuum approaches. Catalytic processes (homogeneous and heterogeneous) can often be reduced to truncated model systems, which are so small that pure quantum-mechanical approaches can be employed. However, since QM/MM becomes more and more efficient due to the success in software and hardware developments, it is more and more used in theoretical organic chemistry to study effects which result from the molecular nature of the environment. It is shown by many examples discussed in this review that the influence can be tremendous, even for nonpolar reactions. The importance of environmental effects in theoretical spectroscopy was already known. Due to its benefits, QM/MM can be expected to experience ongoing growth for the next decade.In the present chapter we give an overview of QM/MM developments and their importance in theoretical organic chemistry, and review applications which give impressions of the possibilities and the importance of the relevant effects. Since there is already a bunch of excellent reviews dealing with QM/MM, we will discuss fundamental ingredients and developments of QM/MM very briefly with a focus on very recent progress. For the applications we follow a similar strategy.

Diffusion of atomic oxygen relevant to water formation in amorphous interstellar ices

Molecular dynamics (MD) simulations together with accurate physics-based force fields are employed to determine the mobility of atomic oxygen in amorphous ice at low temperatures, characteristic for conditions in interstellar ices. From the simulations it is found that the mobility of atomic oxygen ranges from 60 to 480 Ų ns⁻¹ in amorphous ice at temperatures between 50 and 200 K. Hence, the simulations establish that atomic oxygen is mobile to a certain degree and a chemical mechanism for water formation involving oxygen mobility is a realistic scenario. This is also confirmed by the computed migration barriers for oxygen diffusion by multiple umbrella sampling simulations, which yield barriers for diffusion in the range of 0.7–1.9 kcal mol⁻¹. The physics-based force field – based on a multipolar expansion of the electrostatic interactions – yields more pronounced energetics for oxygen migration pathways compared to the conventional point-charge models employed in typical simulations. Once formed, the computed solvation free energy suggests that atomic oxygen thermodynamically prefers to be localized inside amorphous ice and is available for chemical reaction, which may be relevant to water formation in and on grains.

Accretion disc origin of the Earth's water

Earth's water is conventionally believed to be delivered by comets or wet asteroids after the Earth formed. However, their elemental and isotopic properties are inconsistent with those of the Earth. It was thus proposed that water was introduced by adsorption onto grains in the accretion disc prior to planetary growth, with bonding energies so high as to be stable under high-temperature conditions. Here, we show both by laboratory experiments and numerical simulations that water adsorbs dissociatively on the olivine {100} surface at the temperature (approx. 500-1500 K) and water pressure (approx. 10⁻⁸ bar) expected for the accretion disc, leaving an OH adlayer that is stable at least up to 900 K. This may result in the formation of many Earth oceans, provided that a viable mechanism to produce water from hydroxyl exists. This adsorption process must occur in all disc environments around young stars. The inevitable conclusion is that water should be prevalent on terrestrial planets in the habitable zone around other stars.

A first-principles characterization of water adsorption on forsterite grains

Numerical simulations examining chemical interactions of water molecules with forsterite grains have demonstrated the efficacy of nebular gas adsorption as a viable mechanism for water delivery to the terrestrial planets. Nevertheless, a comprehensive picture detailing the water-adsorption mechanisms on forsterite is not yet available. Towards this end, using accurate first-principles density functional theory, we examine the adsorption mechanisms of water on the (001), (100), (010) and (110) surfaces of forsterite. While dissociative adsorption is found to be the most energetically favourable process, two stable associative adsorption configurations are also identified. In dual-site adsorption, the water molecule interacts strongly with surface magnesium and oxygen atoms, whereas single-site adsorption occurs only through the interaction with a surface Mg atom. This results in dual-site adsorption being more stable than single-site adsorption.

The interaction of hydrogen with the {010} surfaces of Mg and Fe olivine as models for interstellar dust grains: a density functional theory study

There is no consensus as yet to account for the significant presence of water on the terrestrial planets, but suggested sources include direct hydrogen adsorption from the parent molecular cloud after the planets' formation, and delivery of hydrous material via comets or asteroids external to the zone of the terrestrial planets. Alternatively, a more recent idea is that water may have directly adsorbed onto the interstellar dust grains involved in planetary formation. In this work, we use electronic structure calculations based on the density functional theory to investigate and compare the bulk and {010} surface structures of the magnesium and iron end-members of the silicate mineral olivine, namely forsterite and fayalite, respectively. We also report our results on the adsorption of atomic hydrogen at the mineral surfaces, where our calculations show that there is no activation barrier to the adsorption of atomic hydrogen at these surfaces. Furthermore, different surface sites activate the atom to form either adsorbed hydride or proton species in the form of hydroxy groups on the same surface, which indicates that these mineral surfaces may have acted as catalytic sites in the immobilization and reaction of hydrogen atoms to form dihydrogen gas or water molecules.

Characterization of hydrogen dissociation over aluminium-doped zinc oxide using an efficient massively parallel framework for QM/MM calculations

A task-farm parallelization framework has been implemented in the ChemShell computational chemistry environment to provide a facility for parallelizing common chemical calculations, including finite-difference Hessian evaluation, the nudged elastic band method for reaction path optimization, and population-based methods for global optimization. The optimization methods are provided by a parallel interface to the DL-FIND optimization library. As ChemShell can already exploit parallel external programs for energy and gradient evaluations, the new methods result in a two-level approach to parallelization that gives significantly improved performance for massively parallel calculations. For typical systems, speed-up factors of five to eight times have been observed compared with non-task-farmed calculations. The task-farming version of ChemShell has been used to study the heterolytic dissociation of a hydrogen molecule over a polar oxygen-terminated surface of aluminium-doped zinc oxide using an embedded cluster hybrid QM/MM approach. We calculate a 42 kcal mol −1 heat of reaction and a 30 kcal mol −1 activation energy, which is equivalent to a high backward reaction barrier of 72 kcal mol −1 per H 2 molecule, in close agreement with temperature programmed desorption experiments. The dissociation path includes a stable intermediate comprising a hydride ion in an oxygen vacancy and physisorbed atomic hydrogen.

DL-FIND: an open-source geometry optimizer for atomistic simulations

Geometry optimization, including searching for transition states, accounts for most of the CPU time spent in quantum chemistry, computational surface science, and solid-state physics, and also plays an important role in simulations employing classical force fields. We have implemented a geometry optimizer, called DL-FIND, to be included in atomistic simulation codes. It can optimize structures in Cartesian coordinates, redundant internal coordinates, hybrid-delocalized internal coordinates, and also functions of more variables independent of atomic structures. The implementation of the optimization algorithms is independent of the coordinate transformation used. Steepest descent, conjugate gradient, quasi-Newton, and L-BFGS algorithms as well as damped molecular dynamics are available as minimization methods. The partitioned rational function optimization algorithm, a modified version of the dimer method and the nudged elastic band approach provide capabilities for transition-state search. Penalty function, gradient projection, and Lagrange-Newton methods are implemented for conical intersection optimizations. Various stochastic search methods, including a genetic algorithm, are available for global or local minimization and can be run as parallel algorithms. The code is released under the open-source GNU LGPL license. Some selected applications of DL-FIND are surveyed.