Search

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.

A pair of simulated left and right circularly polarized ultra-fast laser pulses of duration 20 femtoseconds that induce a mixture of excited states are applied to ethane. The response of the electron dynamics is investigated within the next generation quantum theory of atoms in molecules (NG-QTAIM) using third-generation eigenvector-trajectories which are introduced in this work. This enables an analysis of the mechanical and chiral properties of the electron dynamics of ethane without needing to subject the C-C bond to external torsions as was the case for second-generation eigenvector-trajectories. The mechanical properties, in particular, the bond-flexing and bond-torsion were found to increase depending on the plane of the applied laser pulses. The bond-flexing and bond-torsion, depending on the plane of polarization, increases or decreases after the laser pulses are switched off. This is explainable in terms of directionally-dependent effects of the long-lasting superpositions of excited states. The chiral properties correspond to the ethane molecule being classified as formally achiral consistent with previous NG-QTAIM investigations. Future planned investigations using ultra-fast circularly polarized lasers are briefly discussed., (© 2023 Wiley Periodicals LLC.)

A nonionizing ultrafast laser pulse of 20-fs duration with a peak amplitude electric-field ±E = 200 × 10 -4 a.u. was simulated. It was applied to the ethene molecule to consider its effect on the electron dynamics, both during the application of the laser pulse and for up to 100 fs after the pulse was switched off. Four laser pulse frequencies ω = 0.2692, 0.2808, 0.2830, and 0.2900 a.u. were chosen to correspond to excitation energies mid-way between the (S 1 ,S 2 ), (S 2 ,S 3 ), (S 3 ,S 4 ) and (S 4 ,S 5 ) electronic states, respectively. Scalar quantum theory of atoms in molecules (QTAIM) was used to quantify the shifts of the C1C2 bond critical points (BCPs). Depending on the frequencies ω selected, the C1C2 BCP shifts were up to 5.8 times higher after the pulse was switched off compared with a static E-field with the same magnitude. Next generation QTAIM (NG-QTAIM) was used to visualize and quantify the directional chemical character. In particular, polarization effects and bond strengths, in the form of bond-rigidity vs. bond-flexibility, were found, for some laser pulse frequencies, to increase after the laser pulse was switched off. Our analysis demonstrates that NG-QTAIM, in partnership with ultrafast laser irradiation, is useful as a tool in the emerging field of ultrafast electron dynamics, which will be essential for the design, and control of molecular electronic devices., (© 2023 Wiley Periodicals LLC.)

Active bi-metallic nanoparticles are of key importance in catalysis and renewable energy. Here, the in situ formation of bi-metallic nanoparticles is investigated by exsolution on 200 nm diameter perovskite fibers. The B-site co-doped perovskite fibers display a high degree of exsolution, decorated with NiCo or Ni 3 Fe bi-metallic nanoparticles with average diameter about 29 and 35 nm, respectively. The perovskite fibers are utilized as cathode materials in pure CO 2 electrolysis cells due to their redox stability in the CO/CO 2 atmosphere. After in situ electrochemical switching, the nanoparticles exsolved from the perovskite fiber demonstrate an enhanced performance in pure CO 2 electrolysis. At 900 °C, the current density of solid oxide electrolysis cell (SOEC) with 200 µm YSZ electrolyte supported NiFe doped perovskite fiber anode reaches 0.75 Acm -2 at 1.6 V superior to the NiCo doped perovskite fiber anode (about 1.5 times) in pure CO 2 . According to DFT calculations (PBE-D3 level) the superior CO 2 conversion on NiFe compared to NiCo bi-metallic species is related to an enhanced driving force for C-O cleavage under formation of CO chemisorbed on the nanoparticle and a reduced binding energy of CO required to release this product., (© 2022 The Authors. Small published by Wiley-VCH GmbH.)

We used next-generation QTAIM (NG-QTAIM) to explain the cis -effect for two families of molecules: C 2 X 2 (X = H, F, Cl) and N 2 X 2 (X = H, F, Cl). We explained why the cis -effect is the exception rather than the rule. This was undertaken by tracking the motion of the bond critical point ( BCP ) of the stress tensor trajectories T σ ( s ) used to sample the U σ -space cis- and trans-characteristics. The T σ (s) were constructed by subjecting the C1-C2 BCP and N1-N2 BCP to torsions ± θ and summing all possible T σ (s) from the bonding environment. During this process, care was taken to fully account for multi-reference effects. We associated bond-bending and bond-twisting components of the T σ (s) with cis- and trans- characteristics, respectively, based on the relative ease of motion of the electronic charge density ρ ( r b ). Qualitative agreement is found with existing experimental data and predictions are made where experimental data is not available.

The effect of a static electric (E)-field and an unchirped and chirped laser pulse field on the cycl[3.3.3]azine molecule was investigated using next-generation quantum theory of atoms in molecules (NG-QTAIM). Despite the magnitude of the E-field of the laser pulses being an order of magnitude lower than for the static E-field, the variation of the energy gap between the lowest lying singlet (S 1 ) and triplet (T 1 ) excited states was orders of magnitude greater for the laser pulse than for the static E-field. Insights into the response of the electronic structure were captured by NG-QTAIM, where differences in the inverted singlet-triplet gap due to the laser pulses were significant larger compared to those induced by the static E-field. The response of the S 1 and T 1 excited states, as determined by NG-QTAIM, switched discontinuously between weak and strong chemical character for the static E-field. In contrast, the response to the laser pulses, determined by NG-QTAIM, is to induce a continuous range of chemical character, indicating the unique ability of the laser pulses to induce polarization effects in the form of "mixed" bond types. Our analysis demonstrates that NG-QTAIM is a useful tool for understanding the response to laser irradiation of the lowest-lying singlet S 1 and triplet T 1 excited states of emitters exhibiting thermally activated delayed fluorescence. The chirped laser pulse led to more frequent instances of the desired outcome of an inverted singlet-triplet gap than the unchirped pulse, indicating its usefulness as a tool to design more efficient organic light-emitting diode devices., (© 2021 Wiley Periodicals LLC.)

The association between molecular chirality and helical characteristics known as the chirality-helicity equivalence is determined for the first time by quantifying a chirality-helicity measure consistent with photoexcitation circular dichroism experiments. This is demonstrated using a formally achiral S N 2 reaction and a chiral S N 2 reaction. Both the achiral and chiral S N 2 reactions possess significant values of the chirality-helicity measure and display a Walden inversion, i. e. an inversion of the chirality between the reactant and product. We also track the chirality-helicity measure along the reaction path and discover the presence of chirality at the transition state and two intermediate structures for both reactions. We demonstrate the need for the chirality-helicity measure to differentiate between steric effects due to eclipsed conformations and chiral behaviors in formally achiral species. We explain the significance of this work for asymmetric synthetic reactions including the intermediate structures where the Cahn-Ingold-Prelog (CIP) rules cannot be used., (© 2021 Wiley-VCH GmbH.)

Using a model molecule, we show that it is possible to create molecules that show the required properties for use as elements in a molecular circuit or computer: two conformations with similar energy but different electric conductivity, and the possibility to switch between those by applying an external electric field.

The effect of a varying, directional E x , E y , and E z electric field on the ethene molecule was investigated using next-generation quantum theory of atoms in molecules (QTAIM). Despite using low E-field strengths that are within the realm of experiment and do not measurably alter the molecular geometry, significant changes to the QTAIM properties were observed. Using conventional QTAIM, the shifting of the C─C and C─H bond critical points (BCPs) demonstrates polarization through an interchange in the size of the atoms involved in a bond, since a BCP is located on the boundary between a pair of bonded atoms. Next-generation QTAIM, however, demonstrates the polarization effect more directly with a change in morphology of the 3-D envelope around the BCP. Modest increases of ≈ 2% in the ellipticity ε of the BCP were uncovered when the C─C bond was aligned parallel or anti-parallel to the applied E x -field. Significant asymmetries were found in the response of the next-generation QTAIM 3-D paths of the C─H bonds to the applied E-field. When the E-field coincided with the C─C bond, the BCP moved in response and was accompanied by the envelope constructed from 3-D next-generation paths. The response displayed a polarization effect that increased with increasing magnitude of the E x -field parallel and anti-parallel to the C─C bond. Our analysis demonstrates that next-generation QTAIM is a useful tool for understanding the response of molecules to E-fields, for example, for the screening of molecular wires for the design of molecular electronic devices. © 2019 Wiley Periodicals, Inc., (© 2019 Wiley Periodicals, Inc.)

Next-generation quantum theory of atoms in molecules (QTAIM) was used to investigate the competition between hydrogen bonding and halogen bonding for the recently proposed (Y = Br, I, At)/halogenabenzene/NH 3 complex. Differences between using the SR-ZORA Hamiltonian and effective core potentials (ECPs) to account for relativistic effects with increased atomic mass demonstrated that next-generation QTAIM is a much more responsive tool than conventional QTAIM. Subtle details of the competition between halogen bonding and hydrogen bonding were observed, indicating a mixed chemical character shown in the 3-D paths constructed from the bond-path framework set B. In addition, the use of SR-ZORA reduced or entirely removed spurious features of B on the site of the halogen atoms.

The effect of an electric field on a recently proposed molecular switch based on a quinone analogue was investigated using next-generation quantum theory of atoms in molecules (QTAIM) methodology. The reversal of a homogenous external electric field was demonstrated to improve the "OFF" functioning of the switch. This was achieved by destabilization of the H atom participating in the tautomerization process along the hydrogen bond that defines the switch. The "ON" functioning of the switch, from the position of the tautomerization barrier, is also improved by the reversal of the homogenous external electric field: this result was previously inaccessible. The "ON" and "OFF" functioning of the switch was visualized in terms of the response of the most preferred directions of motion of the electronic charge density to the applied external field. All measures from QTAIM and the stress tensor provide consistent results for the factors affecting the "ON" and "OFF" switch performance. Our analysis therefore demonstrates use for future design of molecular electronic devices. © 2019 Wiley Periodicals, Inc., (© 2019 Wiley Periodicals, Inc.)

The structure of a chloride terminated copper monolayer electrodeposited onto Au(111) from a CuSO 4 /KCl electrolyte was investigated ex situ by three complementary experimental techniques (scanning tunneling microscopy (STM), photoelectron spectroscopy (PES), X-ray standing wave (XSW) excitation) and density functional theory (DFT) calculations. STM at atomic resolution reveals a stable, highly ordered layer which exhibits a Moiré structure and is described by a (5 × 5) unit cell. The XSW/PES data yield a well-defined position of the Cu layer and the value of 2.16 Å above the topmost Au layer suggests that the atoms are adsorbed in threefold hollow sites. The chloride exhibits some distribution around a distance of 3.77 Å in agreement with the observed Moiré pattern due to a higher order commensurate lattice. This structure, a high order commensurate Cl overlayer on top of a commensurate (1 × 1) Cu layer with Cu at threefold hollow sites, is corroborated by the DFT calculations.

A series of para-oligophenylene mono- and dicarboxylic acids (R-(C6H4)nCOOH, n = 1-3, R = H,COOH) was studied. Adsorbed on Au(111)/mica modified by an underpotential deposited bilayer of Ag, the self-assembled monolayers (SAMs) were analyzed by near-edge X-ray absorption fine structure spectroscopy, X-ray photoelectron spectroscopy, and scanning tunneling microscopy. In all cases SAMs are formed with molecules adopting an upright orientation and anchored to the substrate by a carboxylate. Except benzoic acid, all SAMs could be imaged at molecular resolution, which revealed highly crystalline layers with a dense molecular packing. The structures of the SAMs are described by a rectangular (5 × √3) unit cell for the prevailing phase of the monocarboxylic acids and an oblique ([Formula: see text]) unit cell for the dicarboxylic acids, thus evidencing a pronounced influence of the second COOH moiety on the SAM structure. Density functional theory calculations suggest that hydrogen bonding between the SAM-terminating COOH moieties accounts for the difference. Contrasting other classes of SAMs, the systems studied here are determined by intermolecular interactions whereas molecule-substrate interactions play a secondary role. Thus, eliminating problems arising from the mismatch between the molecular and the substrate lattices, coordinatively bonded carboxylic acids on silver should provide considerable flexibility in the design of SAM structures.

A potential end-point in the miniaturization of electronic devices lies in the field of molecular electronics, where molecules perform the function of single components. To date, hydrogen tautomerism in unimolecular switches has been restricted to the central macrocycle of porphyrin-type molecules. The present work reveals how H-tautomerism is the mechanism for switching in substituted quinone derivatives, a novel class of molecules with a different chemical structure. We hence reveal that the previous restrictions applying to tautomeric molecular switches bound to a surface are not valid in general. The activation energy of switching in a prototypical quinone derivative is determined using inelastic electron tunneling. Through computational modeling, we show that the mechanism underlying this process is tautomerization of protons belonging to two amino groups. This switching property is retained upon functionalization by the addition of side groups, meaning that the switch can be chemically modified to fit specific applications.

Benzotriazole (BTAH) has been used as a copper corrosion inhibitor since the 1950s; however, the molecular level detail of how inhibition occurs remains a matter of debate. The onset of BTAH adsorption on a Cu(111) single crystal was investigated via scanning tunnelling microscopy (STM), vibrational spectroscopy (RAIRS) and supporting DFT modelling. BTAH adsorbs as anionic (BTA(-)), CuBTA is a minority species, while Cu(BTA)2, the majority of the adsorbed species, form chains, whose sections appear to diffuse in a concerted manner. The copper surface appears to reconstruct in a (2 × 1) fashion.

The tri-para-carboxylic polychlorotriphenylmethyl radical adsorption on a Cu/Au(111) surface has been investigated in ultra-high vacuum environment. The presence of copper favours the formation of metal-organic assemblies, which are analysed using scanning tunnelling microscopy (STM) and high resolution electron energy loss spectroscopy (HREELS). DFT methods indicate that the unpaired electron survives after adsorption, thus creating systems which could be regarded as potential candidates for spintronics applications.

(51)V NMR parameters have been calculated for VOCl(3), the reference compound in (51)V NMR spectroscopy, in order to capture environmental effects in both the neat liquid and the solid state. Using a combination of periodic geometry optimizations and Car-Parrinello molecular dynamics simulations with embedded cluster NMR calculations, we are able to test the ability of current computational approaches to reproduce (51)V NMR properties (isotropic shifts, anisotropic shifts and quadrupole coupling constants) in the gas, liquid and solid states, for direct comparison with liquid and solid-state experimental data. The results suggest that environmental effects in the condensed phases can be well captured by an embedded cluster approach and that the remaining discrepancy with experiment may be due to the approximate density functionals in current use. The predicted gas-to-liquid shift on the isotropic shielding constant is small, validating the common practice to use a single VOCl(3), molecule as reference in (51)V NMR computations.

Resonance tunneling spectroscopy and density functional theory calculations are employed to explore local variations in the electronic surface potential of a single graphene layer grown on Rh(111). A work function modulation of 220 meV is experimentally measured, indicating that the chemical bonding strength varies significantly across the supercell of the Moiré pattern formed when graphene is bonded to Rh(111). In combination with high-resolution images, which provide precise knowledge of the local atomic registry at the carbon-metal interface, we identify experimentally, and confirm theoretically, the atomic configuration of maximum chemical bonding to the substrate. Our observations are at odds with reported trends for other transition metal substrates. We explain why this is the case by considering the various factors that contribute to the bonding at the graphene/metal interface.

The surface-confined coupling reaction between melamine (1,3,5-triazine-2,4,6-triamine) and 1,4-phenylene diisocyanate has been investigated on Au(111) by scanning tunneling microscopy. Diisocyanate species are stabilized at the edges of melamine arrays and coupling reactions to form small urea oligomers may be initiated at room temperature. These oligomers are incorporated into the two-dimensional melamine array. Annealing accelerates the formation of larger oligomers with multiple urea linkages. The oligomers can themselves form ordered 2-D structures stabilized by intermolecular H-bonding. At higher annealing temperatures, oligomers containing as many as seven or eight urea linkages were identified. These oligomers were able to form 2-D porous structures via interoligomer H-bonding interactions. We discuss the composition of all of the phases observed and identify how covalent and noncovalent interactions stabilize each phase.

The aim of this publication is to present a general strategy to engineer more efficient photoluminescent dendritic molecules based on polyhedral oligomeric silsesquioxane (POSS) cores. A series of chromophores were grafted on POSS cores to form dendritic molecules for which steric hindrance was used as a trigger to tune their photophysical properties. For fluorescence in the blue/near-UV spectral ranges, 4-vinylbiphenyl molecules were chosen as model chromophores to present a general approach based on stable chemistry and bulky groups grafted to the chromophores to enhance photoluminescence efficiency of the dendritic molecules. Photoluminescence quantum yields as well as steady-state and time-resolved solution spectroscopy along with molecular dynamics investigation and electronic structure calculations on a family of new materials are reported. We highlight an apparent contrast between free chromophore and dendritic molecules photophysical properties and show that chromophores' engineering and confinement around an inorganic core allows the design of more efficient photoluminescent dendritic molecules relevant to sensors and hybrid light-emitting diodes.