Ab initio calculation of anharmonic vibrational states of polyatomic systems: Electronic structure combined with vibrational self-consistent field

Performance of Effective Harmonic Oscillator Approach for the Calculations of Vibrational Transition Energies of Large Molecules

The accuracy and performance of the effective harmonic oscillator approximation for the description of anharmonic vibrational structure calculations are tested for large molecular systems and compared with experimental values along with vibrational self-consistent field and second-order perturbation theories. The effective harmonic oscillator approach is an effective single-particle approximation where the variational parameters are the centroids and widths of the multidimensional Gaussian product functions posited as the vibrational wave functions. A comprehensive calculation for 849 transitions that include the fundamentals, two and three quanta overtone transitions, and several combination bands of three polyaromatic hydrocarbons and one DNA nucleobase with a total of 231 normal modes are assessed. A comparison of EHO results with the experimental values is done for the polyaromatic hydrocarbons, and a close agreement is found between the two results. It also offers anharmonic eigenstates and eigenfunctions that are nearly identical with vibrational self-consistent field theory. An extensive analysis on the resultant wave functions of the excited states is performed. The overall root-mean-square deviation (RMSD) between these two methods for 849 transitions understudy is only about 8.3 cm-1, suggesting the effective harmonic oscillator as a viable alternative for the reliable calculations of transition energies of large molecular systems.

Analyses of hydrogen local environments in metals and intermetallic compounds using inelastic neutron scattering calculations based on first-principles hydrogen adiabatic potentials

This study re-evaluates the theoretical approach to analyzing inelastic neutron spectra of hydrogen-containing metals and intermetallic compounds. Previously, these analyses utilized hydrogen quantum nuclear states, modeled as solutions to the Schrödinger equation. The potential surfaces in these models were approximated from the total energies derived from first-principles electronic structure calculations. The current study improves upon this method by employing more efficient and accurate treatments for sampling the potential surface. It utilizes symmetrically irreducible sampling points arranged on densely populated mesh grids for the first-principles calculations. A comparative analysis of the theoretical predictions with experimental spectra for hydrides of Ti2Sb and Ti3Sb, as well as a LaNi5hydrogen primary solid solution, demonstrates that this approach is promising for elucidating the unknown local environments of hydrogen atoms in systems where the approximate potential well describes the hydrogen quantum states.

Exploring Alternate Methods for the Calculation of High-Level Vibrational Corrections of NMR Spin-Spin Coupling Constants

Traditional nuclear magnetic resonance (NMR) calculations typically treat systems with a Born-Oppenheimer-derived electronic wave function that is solved for a fixed nuclear geometry. One can numerically account for this neglected nuclear motion by averaging over property values for all nuclear geometries with a vibrational wave function and adding this expectation value as a correction to an equilibrium geometry property value. Presented are benchmark coupled-cluster singles and doubles (CCSD) vibrational corrections to spin-spin coupling constants (SSCCs) computed at the level of vibrational second-order perturbation theory (VPT2) using the vibrational averaging driver of the CFOUR program. As CCSD calculations of vibrational corrections are very costly, cheaper electronic structure methods are explored via a newly developed Python vibrational averaging program within the Dalton Project. Namely, results obtained with the second-order polarization propagator approximation (SOPPA) and density functional theory (DFT) with the B3LYP and PBE0 exchange-correlation functionals are compared to the benchmark CCSD//CCSD(T) and experimental values. CCSD//CCSD(T) corrections are also combined with literature CC3 equilibrium geometry values to form the highest-order vibrationally corrected values available (i.e., CC3//CCSD(T) + CCSD//CCSD(T)). CCSD//CCSD(T) statistics showed favorable statistics in comparison to experimental values, albeit at an unfavorably high computational cost. A cheaper CCSD//CCSD(T) + B3LYP method showed quite similar mean absolute deviation (MAD) values as CCSD//CCSD(T), concluding that CCSD//CCSD(T) + B3LYP is optimal in terms of cost and accuracy. With reference to experimental values, a vibrational correction was not worth the cost for all of the other methods tested. Finally, deviation statistics showed that CC3//CCSD(T) + CCSD//CCSD(T) vibrational-corrected equilibrium values deteriorated in comparison to CCSD//CCSD(T) attributed to the use of a smaller basis set or lack of solvation effects for the CC3 equilibrium calculations.

A Critical Assessment on Calculating Vibrational Spectra in Nanostructured Materials

Vibrational spectroscopy is an omnipresent spectroscopic technique to characterize functional nanostructured materials such as zeolites, metal-organic frameworks (MOFs), and metal-halide perovskites (MHPs). The resulting experimental spectra are usually complex, with both low-frequency framework modes and high-frequency functional group vibrations. Therefore, theoretically calculated spectra are often an essential element to elucidate the vibrational fingerprint. In principle, there are two possible approaches to calculate vibrational spectra: (i) a static approach that approximates the potential energy surface (PES) as a set of independent harmonic oscillators and (ii) a dynamic approach that explicitly samples the PES around equilibrium by integrating Newton's equations of motions. The dynamic approach considers anharmonic and temperature effects and provides a more genuine representation of materials at true operating conditions; however, such simulations come at a substantially increased computational cost. This is certainly true when forces and energy evaluations are performed at the quantum mechanical level. Molecular dynamics (MD) techniques have become more established within the field of computational chemistry. Yet, for the prediction of infrared (IR) and Raman spectra of nanostructured materials, their usage has been less explored and remain restricted to some isolated successes. Therefore, it is currently not a priori clear which methodology should be used to accurately predict vibrational spectra for a given system. A comprehensive comparative study between various theoretical methods and experimental spectra for a broad set of nanostructured materials is so far lacking. To fill this gap, we herein present a concise overview on which methodology is suited to accurately predict vibrational spectra for a broad range of nanostructured materials and formulate a series of theoretical guidelines to this purpose. To this end, four different case studies are considered, each treating a particular material aspect, namely breathing in flexible MOFs, characterization of defects in the rigid MOF UiO-66, anharmonic vibrations in the metal-halide perovskite CsPbBr3, and guest adsorption on the pores of the zeolite H-SSZ-13. For all four materials, in their guest- and defect-free state and at sufficiently low temperatures, both the static and dynamic approach yield qualitatively similar spectra in agreement with experimental results. When the temperature is increased, the harmonic approximation starts to fail for CsPbBr3 due to the presence of anharmonic phonon modes. Also, the spectroscopic fingerprints of defects and guest species are insufficiently well predicted by a simple harmonic model. Both phenomena flatten the potential energy surface (PES), which facilitates the transitions between metastable states, necessitating dynamic sampling. On the basis of the four case studies treated in this Review, we can propose the following theoretical guidelines to simulate accurate vibrational spectra of functional solid-state materials: (i) For nanostructured crystalline framework materials at low temperature, insights into the lattice dynamics can be obtained using a static approach relying on a few points on the PES and an independent set of harmonic oscillators. (ii) When the material is evaluated at higher temperatures or when additional complexity enters the system, e.g., strong anharmonicity, defects, or guest species, the harmonic regime breaks down and dynamic sampling is required for a correct prediction of the phonon spectrum. These guidelines and their illustrations for prototype material classes can help experimental and theoretical researchers to enhance the knowledge obtained from a lattice dynamics study.

Determining internal coordinate sets for optimal representation of molecular vibration

Arising from the harmonic approximation in solving the vibrational Schrödinger equation, normal modes dissect molecular vibrations into distinct degrees of freedom. Normal modes are widely used as they give rise to descriptive vibrational notations and are convenient for expanding anharmonic potential energy surfaces as an alternative to higher-order Taylor series representations. Usually, normal modes are expressed in Cartesian coordinates, which bears drawbacks that can be overcome by switching to internal coordinates. Considering vibrational notations, normal modes with delocalized characters are difficult to denote, but internal coordinates offer a route to clearer notations. Based on the Hessian, normal mode decomposition schemes for a given set of internal coordinates can describe a normal mode by its contributions from internal coordinates. However, choosing a set of internal coordinates is not straightforward. While the Hessian provides unique sets of normal modes, various internal coordinate sets are possible for a given system. In the present work, we employ a normal mode decomposition scheme to choose an optimal set. Therefore, we screen reasonable sets based on topology and symmetry considerations and rely on a metric that minimizes coupling between internal coordinates. Ultimately, the Nomodeco toolkit presented here generates internal coordinate sets to find an optimal set for representing molecular vibrations. The resulting contribution tables can be used to clarify vibrational notations. We test our scheme on small to mid-sized molecules, showing how the space of definable internal coordinate sets can significantly be reduced.

Ab initio calculations on the X∼2B1 → X∼1A1 photoelectron spectrum of thioformaldehyde negative ion

The equilibrium structures, multidimensional potential energy surfaces, and anharmonic vibrational frequencies of thioformaldehyde (H2CS) and its 2B1 radical anion were obtained at the (U)CCSD(T)-F12a/cc-pVTZ-F12 level; and the photoelectron spectrum of H2CS-X∼2B1→ H2CS X∼1A1 was simulated by calculating the Franck-Condon factors (FCFs). The additional electrons cause a significant change in the bond length of the C = S bond of H2CS, which affects the potential energy surface and the C = S bond stretching mode, while the interaction between electrons turns the originally planar structure of H2CS into a slightly out-of-plane umbrella structure, resulting in an unusual change in the vibrational frequency of the out-of-plane bending mode. Finally, a time-independent eigenstate-free Raman wave function approach (RWF) was used to calculate the photoelectron spectrum of H2CS and to point out the changes in the band under the influence of different temperatures. These data provides a theoretical basis for interstellar observations and experimental studies.

Vibrationally correlated calculations in polyspherical coordinates: Taylor expansion-based kinetic energy operators

The efficiency of quantum chemical simulations of nuclear motion can in many cases greatly benefit from the application of curvilinear coordinate systems. This is rooted in the fact that a set of smartly selected curvilinear coordinates may represent the motion naturally well, thus decreasing the couplings between motions in these coordinates. In this study, we assess the validity of different Taylor expansion-based approximations of kinetic energy operators in a (curvilinear) polyspherical parametrization. To this end, we investigate the accuracy as well as the numerical performance of the approximations in time-independent vibrational coupled cluster and full vibrational interaction calculations for several test cases ranging from tri- to penta-atomic molecules. We find that several of the proposed schemes reproduce the vibrational ground state and excitation energies to a decent accuracy, justifying their application in future investigations. Furthermore, due to the restricted mode coupling and their inherent sum-of-products form, the new approximations open up the possibility of treating large molecular systems with efficient vibrational coupled cluster schemes in general coordinates.

Vibrational heat-bath configuration interaction with semistochastic perturbation theory using harmonic oscillator or VSCF modals

Vibrational heat-bath configuration interaction (VHCI)-a selected configuration interaction technique for vibrational structure theory-has recently been developed in two independent works [J. H. Fetherolf and T. C. Berkelbach, J. Chem. Phys. 154, 074104 (2021); A. U. Bhatty and K. R. Brorsen, Mol. Phys. 119, e1936250 (2021)], where it was shown to provide accuracy on par with the most accurate vibrational structure methods with a low computational cost. Here, we eliminate the memory bottleneck of the second-order perturbation theory correction using the same (semi)stochastic approach developed previously for electronic structure theory. This allows us to treat, in an unbiased manner, much larger perturbative spaces, which are necessary for high accuracy in large systems. Stochastic errors are easily controlled to be less than 1 cm-1. We also report two other developments: (i) we propose a new heat-bath criterion and an associated exact implicit sorting algorithm for potential energy surfaces expressible as a sum of products of one-dimensional potentials; (ii) we formulate VHCI to use a vibrational self-consistent field (VSCF) reference, as opposed to the harmonic oscillator reference configuration used in previous reports. Our tests are done with quartic and sextic force fields, for which we find that with VSCF, the minor improvements to accuracy are outweighed by the higher computational cost associated the matrix element evaluations. We expect VSCF-based VHCI to be important for more general potential representations, for which the harmonic oscillator basis function integrals are no longer analytic.

CRYSTAL23: A Program for Computational Solid State Physics and Chemistry

The Crystal program for quantum-mechanical simulations of materials has been bridging the realm of molecular quantum chemistry to the realm of solid state physics for many years, since its first public version released back in 1988. This peculiarity stems from the use of atom-centered basis functions within a linear combination of atomic orbitals (LCAO) approach and from the corresponding efficiency in the evaluation of the exact Fock exchange series. In particular, this has led to the implementation of a rich variety of hybrid density functional approximations since 1998. Nowadays, it is acknowledged by a broad community of solid state chemists and physicists that the inclusion of a fraction of Fock exchange in the exchange-correlation potential of the density functional theory is key to a better description of many properties of materials (electronic, magnetic, mechanical, spintronic, lattice-dynamical, etc.). Here, the main developments made to the program in the last five years (i.e., since the previous release, Crystal17) are presented and some of their most noteworthy applications reviewed.

Accuracy of Different Electronic Basis Set Families for Anharmonic Molecular Vibrations: A Comprehensive Benchmark Study

In this work, the accuracy and convergence of different electronic basis set families for the computation of anharmonic molecular vibrational spectroscopic calculations are benchmarked. A series of 39 different basis sets from different families following their hierarchy are assessed on VSCF and VSCF-PT2 algorithms with commonly used MP2 and DFT based B3LYP-D potentials for a set of molecular systems. Such an effort has been validated in a previous work ( J. Phys. Chem. A 2020, 124, 9203-9221) with split-valence basis sets for fundamentals and intensities. Here, fundamental transitions, vibrationally excited states, and intensities are compared with the experimental data to estimate the accuracy for a series of Jensen, Dunning, Calendar, Karlsruhe, and Sapporo basis set families. The convergence of basis sets are also compared with the large ANO basis set. Comprehensive statistical error analysis in terms of accuracy and precision was carried out to assess the performance of each basis set. It is observed that the improvement for the calculated harmonic and anharmonic values from the smaller basis sets to the medium (i.e., triple-ξ) is considerable. Beyond this, from medium to large basis sets, the convergence is slow and mostly posits nearly converged values. Basis sets with and without diffuse functions offer characteristically different accuracies and convergence patterns. Finally, recommendations are given on the choice of basis set chosen as black-box which can balance between accuracy and computational time, estimation of the errors, and their selections especially for large molecules.

Accessing Position Space Wave Functions in Band Structure Calculations of Periodic Systems─A Generalized, Adapted Numerov Implementation for One-, Two-, and Three-Dimensional Quantum Problems

In this work, a generalized, adapted Numerov implementation capable of determining band structures of periodic quantum systems is outlined. Based on the input potential, the presented approach numerically solves the Schrödinger equation in position space at each momentum space point. Thus, in addition to the band structure, the method inherently provides information about the state functions and probability densities in position space at each momentum space point considered. The generalized, adapted Numerov framework provided reliable estimates for a variety of increasingly complex test suites in one, two, and three dimensions. The accuracy of the proposed methodology was benchmarked against results obtained for the analytically solvable Kronig-Penney model. Furthermore, the presented numerical solver was applied to a model potential representing a 2D optical lattice being a challenging application relevant, for example, in the field of quantum computing.

The first HyDRA challenge for computational vibrational spectroscopy

Vibrational spectroscopy in supersonic jet expansions is a powerful tool to assess molecular aggregates in close to ideal conditions for the benchmarking of quantum chemical approaches. The low temperatures achieved as well as the absence of environment effects allow for a direct comparison between computed and experimental spectra. This provides potential benchmarking data which can be revisited to hone different computational techniques, and it allows for the critical analysis of procedures under the setting of a blind challenge. In the latter case, the final result is unknown to modellers, providing an unbiased testing opportunity for quantum chemical models. In this work, we present the spectroscopic and computational results for the first HyDRA blind challenge. The latter deals with the prediction of water donor stretching vibrations in monohydrates of organic molecules. This edition features a test set of 10 systems. Experimental water donor OH vibrational wavenumbers for the vacuum-isolated monohydrates of formaldehyde, tetrahydrofuran, pyridine, tetrahydrothiophene, trifluoroethanol, methyl lactate, dimethylimidazolidinone, cyclooctanone, trifluoroacetophenone and 1-phenylcyclohexane-cis-1,2-diol are provided. The results of the challenge show promising predictive properties in both purely quantum mechanical approaches as well as regression and other machine learning strategies.

Computational infrared and Raman spectra by hybrid QM/MM techniques: a study on molecular and catalytic material systems

Vibrational spectroscopy is one of the most well-established and important techniques for characterizing chemical systems. To aid the interpretation of experimental infrared and Raman spectra, we report on recent theoretical developments in the ChemShell computational chemistry environment for modelling vibrational signatures. The hybrid quantum mechanical and molecular mechanical approach is employed, using density functional theory for the electronic structure calculations and classical forcefields for the environment. Computational vibrational intensities at chemical active sites are reported using electrostatic and fully polarizable embedding environments to achieve more realistic vibrational signatures for materials and molecular systems, including solvated molecules, proteins, zeolites and metal oxide surfaces, providing useful insight into the effect of the chemical environment on the signatures obtained from experiment. This work has been enabled by the efficient task-farming parallelism implemented in ChemShell for high-performance computing platforms.  This article is part of a discussion meeting issue 'Supercomputing simulations of advanced materials'.

Molecular Structures and Vibrational Spectra of trans- and cis-Polyacetylene and Their Oligoenes Revisited Using Density Functional Theory Calculations

Polyacetylene, the most representative synthetic conducting polymer, has attracted much attention because it exhibits high conductivity upon doping. In this paper, molecular structures, electronic excitation energies, and Raman and infrared spectra were calculated using density functional theory for trans- and cis-oligoenes with various chain lengths up to the number of C═C bonds (n) of 100 and trans- and cis-polyacetylenes under one-dimensional periodic boundary condition. The harmonic vibrational frequencies obtained at the B3LYP/6-311G(d,p) level were scaled by the scaling factors determined with respect to the anharmonic vibrational frequencies using the B2PLYP method, in which the coefficients of the functional were optimized for trans-oligoenes. The calculated infrared and Raman frequencies reproduce reasonably well the observed frequencies for trans- and cis-polyacetylene. Based on the chain-length dependence of the calculated Raman spectra of trans-oligoenes, we proposed the possibility of longer conjugated trans-segments observed in the resonance Raman spectra of trans-polyacetylene excited at longer wavelengths of 647.1 and 1064 nm. We also elucidated the origin of the excitation-wavelength dependence of the resonance Raman spectra of trans-polyacetylene and the structure of isomerization intermediates from cis-form to trans-form. In addition, the previous assignments of Raman and infrared spectra of trans- and cis-polyacetylene were reexamined in the present study based on the chain-length dependence of the spectra.

Performance of Vibrational Self-Consistent Field Theory for Accurate Potential Energy Surfaces: Fundamentals, Excited States, and Intensities

The performance of vibrational structure calculations beyond harmonic approximation in the framework of the vibrational self-consistent field method with second-order perturbation corrections (VSCF-PT2) is investigated in conjunction with very accurate potential energy surfaces (PESs) given by various coupled-cluster electronic structure theories. The quality of anharmonic calculations depends on the accuracy of the underlying multidimensional PES obtained from its functional form, which is given by the level of electronic structure theory. Two such highest levels of typical coupled-cluster electronic structure methods, CCSD and the ″gold standard″ CCSD(T), along with their variants such as CCD, CR-CCL (completely renormalized CR-CC(2,3) approach), and CCSD(TQ) are tested for the construction of accurate anharmonic potentials without any fitting or ad hoc scaling and using cc-pVTZ basis sets. The accuracy of VSCF-PT2 theory in comparison to experimental values is tested for a series of 16 molecules with 135 fundamental bands, 64 overtones, and combination bands and also for 39 intensities. It is found that CCD and CCSD bind the potential tighter than CCSD(T) and the computed VSCF-PT2 transitions are more blue-shifted showing higher deviation from the experiment. In general, VSCF-PT2 results computed at the CCSD(T) potential offer a good cost/accuracy ratio, with the mean absolute deviation and the mean absolute percentage error with the experiment being ∼16 cm^-1 and 1.38, respectively, for fundamentals. Additionally, while the CR-CCL and CCSD(TQ) methods offer similar levels of accuracies as compared to CCSD(T), the former offers a better accuracy/cost ratio than the latter and is a suitable alternative to CCSD(T).

A high-accuracy machine-learning water model for exploring water nanocluster structures

Water, the most important molecule on the Earth, possesses many essential and unique physical properties that are far from completely understood, partly due to serious difficulties in identifying the precise microscopic structures of water. Hence, identifying the structures of water nanoclusters is a fundamental and challenging issue for studies on the relationship between the macroscopic physical properties of water and its microscopic structures. For large-scale simulations (at the level of nm and ns) of water nanoclusters, a calculation method with simultaneous accuracy at the level of quantum chemistry and efficiency at the level of an empirical potential method is in great demand. Herein, a machine-learning (ML) water model was utilized to explore the microscopic structural features at different length scales for water nanoclusters with a size up to several nm. The ML water model can be employed to efficiently predict the structures of water nanoclusters with a similar accuracy to that of density functional theory and with substantially lower computational resource demands. To validate the low-lying structure search results with experimental spectral results, an ML water model combined with velocity autocorrelation function analysis was used to simulate the vibrational spectra of water nanoclusters with up to thousands of water molecules. By comparing the simulated and experimentally recorded vibrational spectra, the atomic structures determined by a simulation based on the ML water model are all verified. To demonstrate its ability to represent water's structural evolution at large length and time scales, the ML water model was employed to model the structural evolution during the crystal-liquid transition, and the phase transition temperatures of water clusters with different sizes were precisely predicted. The ML water model provides an efficient theoretical calculation tool for exploring the structures and physical properties of water and their relationships, especially for clusters with relatively large sizes and processes with relatively long durations.

My Trajectory in Molecular Reaction Dynamics and Spectroscopy

This is the story of a career in theoretical chemistry during a time of dramatic changes in the field due to phenomenal growth in the availability of computational power. It is likewise the story of the highly gifted graduate students and postdoctoral fellows that I was fortunate to mentor throughout my career. It includes reminiscences of the great mentors that I had and of the exciting collaborations with both experimentalists and theorists on which I built much of my research. This is an account of the developments of exciting scientific disciplines in which I was involved: vibrational spectroscopy, molecular reaction mechanisms and dynamics, e.g., in atmospheric chemistry, and the prediction of new, exotic molecules, in particular noble gas molecules. From my very first project to my current work, my career in science has brought me the excitement and fascination of research. What a wonderful pursuit!

Exploring the anharmonic vibrational structure of carbon dioxide trimers

Our previously developed mbCO2 potential [O. Sode and J. N. Cherry, J. Comput. Chem. 38, 2763 (2017)] is used to describe the vibrational structure of the intermolecular motions of the CO₂ trimers: barrel-shaped and cyclic trimers. Anharmonic corrections are accounted for using the vibrational self-consistent field theory, vibrational second-order Møller-Plesset perturbation (VMP2) theory, and vibrational configuration interaction (VCI) methods and compared with experimental observations. For the cyclic structure, we revise the assignments of two previously observed experimental peaks based on our VCI and VMP2 results. We note that the experimental band observed near 13 cm^-1 is the out-of-phase out-of-plane degenerate motion with E″ symmetry, while the peak observed at 18 cm^-1 likely corresponds to the symmetric out-of-plane torsion A″ vibration. Since the VCI treatment of the vibrational motions accounts for vibrational mixing and delocalization, overtones and combination bands were also observed and quantified in the intermolecular regions of the two trimer isomers.

Vibrational heat-bath configuration interaction

We introduce vibrational heat-bath configuration interaction (VHCI) as an accurate and efficient method for calculating vibrational eigenstates of anharmonic systems. Inspired by its origin in electronic structure theory, VHCI is a selected CI approach that uses a simple criterion to identify important basis states with a pre-sorted list of anharmonic force constants. Screened second-order perturbation theory and simple extrapolation techniques provide significant improvements to variational energy estimates. We benchmark VHCI on four molecules with 12-48 degrees of freedom and use anharmonic potential energy surfaces truncated at fourth and sixth orders. When compared to other methods using the same truncated potentials, VHCI produces vibrational spectra of tens or hundreds of states with sub-wavenumber accuracy at low computational cost.

Ionic Liquid Clusters Generated from Electrospray Thrusters: Cold Ion Spectroscopic Signatures of Size-Dependent Acid-Base Interactions

We determine the intramolecular distortions at play in the 2-hydroxyethylhydrazinium nitrate (HEHN) ionic liquid (IL) propellant, which presents the interesting case that the HEH⁺ cation has multiple sites (i.e., hydroxy, primary amine, and secondary ammonium groups) available for H-bonding with the nitrate anion. These interactions are quantified by analyzing the vibrational band patterns displayed by cold cationic clusters, (HEH⁺)_n(NO₃^-)_n-1, n = 2-6, which are obtained using IR photodissociation of the cryogenically cooled, mass-selected ions. The strong interaction involving partial proton transfer of the acidic N-H proton in HEH⁺ cation to the nitrate anion is strongly enhanced in the ternary n = 2 cluster but is suppressed with increasing cluster size. The cluster spectra recover the bands displayed by the bulk liquid by n = 5, thus establishing the minimum domain required to capture this aspect of macroscopic behavior.

Multiresolution continuous wavelet transform for studying coupled solute-solvent vibrations via ab initio molecular dynamics

Vibrational analysis in solution and the theoretical determination of infrared and Raman spectra are of key importance in many fields of chemical interest. Vibrational band dynamics of molecules and their sensitivity to the environment can also be captured by these spectroscopies in their time dependent version. However, it is often difficult to provide an interpretation of the experimental data at the molecular scale, such as molecular mechanisms or the processes hidden behind them. In this work, we present a theoretical-computational protocol based on ab initio molecular dynamics simulations and a combination of normal-like (generalized) mode analysis of solute-solvent clusters with a wavelet transform, for the first time. The case study is the vibrational dynamics of N-methyl-acetamide (NMA) in water solution, a well-known model of hydration of peptides and proteins. Amide modes are typical bands of peptide and protein backbone, and their couplings with the environment are very challenging in terms of the accurate prediction of solvent induced intensity and frequency shifts. The contribution of water molecules surrounding NMA to the composition of generalized and time resolved modes is introduced in our vibrational analysis, showing unequivocally its influence on the amide mode spectra. It is also shown that such mode compositions need the inclusion of the first shell solvent molecules to be accurately described. The wavelet analysis is proven to be strongly recommended to follow the time evolution of the spectra, and to capture vibrational band couplings and frequency shifts over time, preserving at the same time a well-balanced time-frequency resolution. This peculiar feature also allows one to perform a combined structural-vibrational analysis, where the different strengths of hydrogen bond interactions can quantitatively affect the amide bands over time at finite temperature. The proposed method allows for the direct connection between vibrational modes and local structural changes, providing a link from the spectroscopic observable to the structure, in this case the peptide backbone, and its hydration layouts.

Hardware efficient quantum algorithms for vibrational structure calculations

We introduce a framework for the calculation of ground and excited state energies of bosonic systems suitable for near-term quantum devices and apply it to molecular vibrational anharmonic Hamiltonians. Our method supports generic reference modal bases and Hamiltonian representations, including the ones that are routinely used in classical vibrational structure calculations. We test different parametrizations of the vibrational wavefunction, which can be encoded in quantum hardware, based either on heuristic circuits or on the bosonic Unitary Coupled Cluster Ansatz. In particular, we define a novel compact heuristic circuit and demonstrate that it provides a good compromise in terms of circuit depth, optimization costs, and accuracy. We evaluate the requirements, number of qubits and circuit depth, for the calculation of vibrational energies on quantum hardware and compare them with state-of-the-art classical vibrational structure algorithms for molecules with up to seven atoms.

NWChem: Past, present, and future

Specialized computational chemistry packages have permanently reshaped the landscape of chemical and materials science by providing tools to support and guide experimental efforts and for the prediction of atomistic and electronic properties. In this regard, electronic structure packages have played a special role by using first-principle-driven methodologies to model complex chemical and materials processes. Over the past few decades, the rapid development of computing technologies and the tremendous increase in computational power have offered a unique chance to study complex transformations using sophisticated and predictive many-body techniques that describe correlated behavior of electrons in molecular and condensed phase systems at different levels of theory. In enabling these simulations, novel parallel algorithms have been able to take advantage of computational resources to address the polynomial scaling of electronic structure methods. In this paper, we briefly review the NWChem computational chemistry suite, including its history, design principles, parallel tools, current capabilities, outreach, and outlook.

Vibrationally resolved emission spectra of luminescent conjugated oligothiophenes from anharmonic calculations

We report on accurate and efficient calculations of vibrationally resolved emission spectra for oligothiophenes from anharmonic vibrational configuration interaction wave-function calculations in reduced vibrational spaces. These reduced spaces are chosen based on the independent mode displaced harmonic oscillator model. Good agreement with experiment is obtained for all-trans oligothiophenes with two to five rings also when employing only a few active modes. Vibrational modes incorporating inter-ring carbon-carbon stretches and a ring breathing mode are found to be the main players in the vibrational progression for the emission from the first excited electronic state for all investigated oligothiophene derivatives. The presented framework is here illustrated for oligothiophenes, but we have made no underlying system-dependent assumptions and believe it to become a valuable tool for the rational design of fluorescence biomarkers.

Understanding the anharmonic vibrational structure of the carbon dioxide dimer

Understanding the vibrational structure of the CO₂ system is important to confirm the potential energy surface and interactions in such van der Waals complexes. In this work, we use our previously developed mbCO2 potential function to explore the vibrational structure of the CO₂ monomer and dimer. The potential function has been trained to reproduce the potential energies at the CCSD(T)-F12b/aug-cc-pVTZ level of electronic structure theory. The harmonic approximation, as well as anharmonic corrections using vibrational structure theories such as vibrational self-consistent field, vibrational second-order Møller-Plesset perturbation, and vibrational configuration interaction (VCI), is applied to address the vibrational motions. We compare the vibrational results using the mbCO2 potential function with traditional electronic structure theory results and to experimental frequencies. The anharmonic results for the monomer most closely match the experimental data to within 3 cm^-1, including the Fermi dyad frequencies. The intermolecular and intramolecular dimer frequencies were treated separately and show good agreement with the most recent theoretical and experimental results from the literature. The VCI treatment of the dimer vibrational motions accounts for vibrational mixing and delocalization, such that we observe the dimer Fermi resonance phenomena, both in the intramolecular and intermolecular regions.

Structural sensitivity of CH vibrational band in methyl benzoate

The CH vibrational bands of methyl benzoate are studied to understand its coupling pattern with other vibrational bands of the biological molecule. This will facilitate to understand the biological structure and dynamics in spectroscopic as well as in microscopic study. Due to the congested spectroscopic pattern, near degeneracy, and strong anharmonicity of the CH stretch vibrations, assignment of the CH vibrational frequencies are often misleading. Anharmonic vibrational frequency calculation with multidimensional potential energy surface interprets the CH vibrational spectra more accurately. In this article we have presented the importance of multidimensional potential energy surface in anharmonic vibrational frequency calculation and discuss the unexpected red shift of asymmetric CH stretch vibration of methyl group. The CD stretch vibrational band which is splitted to double peaks due to the Fermi resonance is also discussed here.

Towards the SMART workflow system for computational spectroscopy

Is it possible to convert highly specialized research in the field of computational spectroscopy into robust and user-friendly aids to experiments and industrial applications?

Anharmonic vibrational spectra from double incremental potential energy and dipole surfaces

Using incremental approaches, size limitations for property surface generations are pushed significantly, enabling accurate large molecule anharmonic vibrational spectra calculations.