334 Search Results

Large atomic-orbital (AO) basis sets of at least triple and preferably quadruple-ζ (QZ) size are required to adequately converge Kohn–Sham density functional theory (DFT) calculations toward the complete basis set limit. However, incrementing the cardinal number by one nearly doubles the AO basis dimension, and the computational cost scales as the cube of the AO dimension, so this is very computationally demanding. Here, in this work, we develop and test a threshold-based natural atomic orbital (NAO) scheme in which ϵ-NAOs are obtained as eigenfunctions of atomic blocks of the density matrix in a one-center orthogonalized representation. This enables compression ofmore » the AO basis that is optimal for a given threshold, 10^–ϵ, by discarding NAOs with occupation numbers below that threshold. Extensive pilot test calculations using the Hartree–Fock functional and taking the converged density matrix as input suggest that a threshold of 10^–5 can yield a compression factor (ratio of AO to compressed ϵ-NAO dimension) between 2.5 and 4.5 for the QZ pc-3 basis. The errors in relative energies are typically less than 0.1 kcal/mol when the compressed basis is used instead of the uncompressed basis. Between 10 and 100 times smaller errors (i.e., usually less than 0.01 kcal/mol) can be obtained with a threshold 10^–7, while the compression factor is typically between 2 and 2.5.« less

Carbon–carbon (C–C) bond formation is necessary for hydrocarbon (and oxygenate) synthesis beyond methane (and formate/formic acid) during electrochemical CO and CO₂ reduction (ECOR and ECO₂R). Cu has notable ability to form hydrocarbons compared to other pure metals. In particular, the (100) facet of face-centered cubic Cu forms ethylene competitively with H₂ and methane during both ECOR and ECO₂R. Past simulations based on density functional theory (DFT) with standard exchange-correlation functional approximations predict fast nonelectroactive C–C bond formation channels involving adsorbed (*) CO together with another *CO, formyl (*CHO), or hydroxymethylidyne (*COH), forming OC*–*CO, OC*–CHO*, and OC*–*COH, respectively. Such simulations supportmore » the prevailing hypothesis that emergence of C₂ products is kinetically determined at the early stages of the reduction chemistry. Here we show, via simulations with more accurate many-body, i.e., “correlated”, wavefunction theory (enabled by an embedding scheme), that the coupling of *CO with a *CO or a *COH (previously predicted at the same level of theory to kinetically dominate over *CHO as the one-electron reduction product of *CO) is highly activated (kinetically impeded), with free energy barriers >1 eV, in contradiction to previous DFT-based simulations. Intriguingly, we find that the coupling of two adjacent *COHs incurs only a small barrier (<0.3 eV) and is exoergic (< –1 eV); however, given the predicted low surface mobility of *COH, the emergence of HOC*–*COH is also improbable, at least at low *COH coverages. We therefore conclude that it is highly unlikely for *CO to participate in nonelectroactive C–C bond formation on pristine Cu(100), contrary to conventional wisdom, and that the energetically favorable *COH dimerization may occur only after substantial buildup of *COH on the surface.« less

Amines play several key roles in chemistry and biology and are involved in numerous industrial processes, often with significant economic impacts. Recently, amines are also garnering interest as catalysts for polymer synthesis and for CO₂ fixation, incentivizing the need to rapidly design and screen new amino compounds. Hence, developing reliable methods to predict their physicochemical properties, e.g., the base dissociation constant (pKb), is pivotal. Here, a density functional theory (DFT)-based approach was employed to compute the pKb of substituted amines, exploring the impact of several key parameters, including (i) the number of explicit water molecules at the reaction center, (ii)more » the van der Waals (vdW) surface, and (iii) solvent polarizability. In previous work, it was determined that including two explicit water molecules at the reaction center resulted in highly accurate pKb estimates for primary amines. Here, we find that including a third water molecule at the reaction center is essential for accurate pKb for secondary and tertiary amines. The revised methodology was then applied to a wider selection of amines, obtaining a minimum average error (MAE) < 0.4. In conclusion, this result represents an extension of our “easy-to-use method,” a simple and direct DFT approach exploiting CAM-B3LYP/SMD/6-311G+(d,p) to compute pKb without post facto modifications.« less

Metal–organic frameworks (MOFs), with their distinctive porous structures and tunable chemical properties, have shown immense promise in the separation and storage of gases. Currently, the accurate simulation of their adsorptive properties remains challenging, especially for systems where the molecules fit very tightly into the pores. Traditional simulation methods often approximate the frameworks as rigid and do not account for the framework flexibility seen in materials such as NbOFFIVE-1-Ni. First-principles molecular dynamics (FPMD) simulations offer the desired accuracy in modeling this flexibility but are limited by their extensive computational demands, rendering them impractical for long simulations. Conversely, classical force field-based simulationsmore » offer computational efficiency but lack the necessary accuracy. Here, to break this accuracy-efficiency trade-off, we have developed machine learning interatomic potentials trained on energies and forces from FPMD to model the framework flexibility of NbOFFIVE-1-Ni in the presence of water over nanosecond time scales. Furthermore, by integrating MLIP-driven molecular dynamics (MLIP-MD) with grand canonical Monte Carlo (GCMC) simulations, we further incorporated framework flexibility into adsorption predictions, yielding water adsorption isotherms that better align with experimental data compared to those of conventional GCMC simulations. These advances offer new opportunities for the design and optimization of MOFs in gas storage and separation applications.« less

Vapor-phase infiltration (VPI) of inorganic materials in polymers is increasingly becoming popular for synthesizing various functional hybrid materials. While AlO_x infiltration using trimethylaluminum (TMA) has been extensively studied, the mechanism of diethylzinc (DEZ)-based ZnO_x infiltration, especially one that is initiated by AlO_x priming, has not received much attention because highly reactive hydroxyl groups generated by AlO_x-priming are expected to dominate the initial binding of DEZ, thus enabling the overall ZnO_x VPI. Here, we interrogate the ZnO_x infiltration mechanism in AlO_x-primed poly(methyl methacrylate) (PMMA) in comparison to the control AlO_x-only infiltration by utilizing a suite of complementary characterizations, including quartz crystalmore » microbalance mass gain measurement, transmission electron microscopy, infrared reflection–absorption spectroscopy (IRRAS), and synchrotron X-ray absorption spectroscopy (XAS). The multivalent TMA precursor and associated hyperbranched AlO_x network can quickly saturate the AlO_x infiltration by clogging the polymer-free volume near the top. On the contrary, the ZnO_x infiltration using divalent DEZ precursor, once activated via AlO_x-priming, can lead to accelerated ZnOx infiltration. With the help of IRRAS, XAS, and density functional theory (DFT) simulations, we uncover that the AlO_x-priming enhances the reactivity of neighboring carbonyl groups toward DEZ and opens up simultaneous reaction pathways, leading to accelerated high-fidelity infiltration of ZnO_x.« less

Key electronic processes related to molecular excitonic states of finite stacks of indolonaphthyridine molecules are analyzed via the non-orthogonal configuration interaction with fragments (NOCI-F) method. Indolonaphthyridine is an organic chromophore that can undergo several electronic photoexcitation-related intermolecular processes, such as exciton and electron transfer. The structures studied here are noncrystalline arrangements built as either ordered stacks of indolonaphthyridine or stacks extracted from molecular dynamics simulations including thermal disorder. Taking dimers or trimers from either model, we performed CASSCF and NOCI-F calculations to quantify the intermolecular electronic couplings governing singlet fission, excited singlet and triplet diffusion, and hole and electron diffusionmore » processes. Also, comparing the results for the different models, we studied the effect of structural disorder and distortion on these couplings. Finally, we present a newly developed, advanced postanalysis tool. It takes the NOCI-F data as input to carry out a multifragment full Hamiltonian procedure that involves the complete stack, providing physical information not available from the dimer/trimer models, hence giving access to additional insight into the material’s properties.« less

Covalent organic frameworks (COFs), with their modular architectures and tunable functionalities, provide a versatile platform to design sorbents for the direct capture of CO₂ from air. Here, for this work, we combined density functional theory, molecular dynamics, and grand canonical Monte Carlo simulations with experiment to understand structural factors for furthering COF-999-NH₂’s performance as the precursor to COF-999 for direct air CO₂ capture. Small energy differences among laterally shifted stackings suggest intrinsic stacking heterogeneity. The simulations show pronounced layer buckling coupled to extensive amine–nitrile hydrogen bonding and persistent pore water, which initiates undesired polymerization and undermines uptake. The predicted presencemore » of water is confirmed by subsequent experiments. These insights point to a single, actionable design rule: exclude retained water by introducing hydrophobic pore environments to maximize the CO₂ capture efficiency.« less

Coupled cluster singles and doubles with perturbative triples [CCSD(T)] and single determinant fixed-node diffusion Monte Carlo (SD-DMC) have emerged as two of the most useful methods for providing benchmark reaction and interaction energies of chemical systems without strong static correlation. The errors in DMC energies are dominated by an inexact description of the nodal surfaces for electron exchange. One of the main approaches to addressing the fixed-node error is to use multideterminant (MD) trial wave functions. We consider here the energy differences between pairs of related molecules with aromatic and quinoidal structures as well as between quinoidal isomers. Quinoidal systemsmore » tend to have some diradical character, leading one to anticipate that SD-DMC calculations may face challenges in accurately describing their energetics. The MD trial wave functions were generated from the complete active space calculations. A comparison is made with the predictions of well-converged CCSD(T) calculations.« less

We have examined the electronic structure of C₅H₅MgMgC₅H₅, or dimagnesocene, using high-level coupled-cluster techniques. This research is suitable in light of the remarkable synthesis of the valence-isoelectronic diberyllocene by Boronski, Crumpton, Wales, and Aldridge. The Mg–Mg bond distance is predicted to be 2.758 Å, and the Mg–Mg bond dissociation energy is predicted to be 51.8 kcal/mol. Unique aspects to the present research is our characterization of the ionization energy and the electron affinity of this molecule, the energy of dissociation into two neutral cyclopentadienyl magnesium radicals, the determination of the neutral structure at the CCSD(T)/cc-pVTZ level, and the computation ofmore » the Raman intensities at the MP2/cc-pVDZ level. Apart from mass spectroscopy, the simplest means of experimental detection is gas-phase or matrix-isolated infrared spectroscopy, in which the A₂″ peak at 801 cm^–1 should be the most prominent, with an intensity of 581 km/mol.« less

Understanding whether a polaritonic phenomenon is fundamentally quantum or classical is essential for building accurate theoretical models and guiding experimental design. Here, in this work, we address this question in the context of polaritonic spectra and report an intriguing new feature: the twin polariton, an additional splitting beyond the primary resonant polariton splitting originating from vacuum field fluctuations. We show that the twin polariton persists in the many-molecule limit under permutationally symmetrical initial-state constraint and that it follows the same linear dependence on coupling strength as the primary polariton splitting. This establishes a novel mechanism by which a quantum featuremore » (the twin polariton) can be tuned through a classical one (the primary polariton), offering new opportunities to probe and control the fundamental nature of polaritonic systems.« less