Engineering Papers⌕ Search

Engineering topics

Schlegel, H. Bernhard

Publications and source records attributed to Schlegel, H. Bernhard.

Roadmap on methods and software for electronic structure based simulations in chemistry and materials

This Roadmap article provides a succinct, comprehensive overview of the state of electronic structure methods and software for molecular and materials simulations. Seventeen distinct sections collect insights by 51 leading scientists in the field. Each contribution addresses the status of a particular area, as well as current challenges and anticipated future advances, with a particular eye towards software related aspects and providing key references for further reading. Foundational sections cover density functional theory and its implementation in real-world simulation frameworks, Green's function based many-body perturbation theory, wave-function based and stochastic electronic structure approaches, relativistic effects and semiempirical electronic structure theory approaches. Subsequent sections cover nuclear quantum effects, real-time propagation of the electronic structure, challenges for computational spectroscopy simulations, and exploration of complex potential energy surfaces. The final sections summarize practical aspects, including computational workflows for complex simulation tasks, the impact of current and future high-performance computing architectures, software engineering practices, education and training to maintain and broaden the community, as well as the status of and needs for electronic structure based modeling from the vantage point of industry environments. Overall, the field of electronic structure software and method development continues to unlock immense opportunities for future scientific discovery, based on the growing ability of computations to reveal complex phenomena, processes and properties that are determined by the make-up of matter at the atomic scale, with high precision.

36 MATERIALS SCIENCE↗

A bioinspired cobalt catalyst based on a tripodal imidazole/pyridine platform capable of water reduction and oxidation

The prototypical drug carrier [Co II (L 1 )Cl]PF 6 (1), where L 1 is a tripodal amine bound to pyridine and methyl-imidazoles, had its electrocatalytic water splitting activity studied under different pH conditions. This species contains a high-spin 3d 7 Co II metal center, and is capable of generating both H 2 from water reduction and O 2 from water oxidation. Turnover numbers reach 390 after 3 h for water reduction. Initial water oxidation activity is molecular, with TONs of 71 at pH 7 and 103 at pH 11.5. The results reveal that species 1 can undergo several redox transformations, including reduction to the 3d 8 Co I species that precedes a LS 3d 6 hydride for water reduction, as well as nominal Co IV $=$O and Co III -OOH species required for water oxidation. Finally, post-catalytic analyses confirm the molecular nature of reduction and support initial molecular activity for oxidation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Ionization of HCCI Neutral and Cations by Strong Laser Fields Simulated With Time Dependent Configuration Interaction

Strong field ionization of neutral iodoacetylene (HCCI) can produce a coherent superposition of the X and A cations. This superposition results in charge migration between the CC π orbital and the iodine π-type lone pair which can be monitored by strong field ionization with short, intense probe pulses. Strong field ionization of the X and A states of HCCI cation was simulated with time-dependent configuration interaction using singly ionized configurations and singly excited, singly ionized configurations (TD-CISD-IP) and an absorbing boundary. Studies with static fields were used to obtain the 3-dimensional angular dependence of instantaneous ionization rates by strong fields and the orbitals involved in producing the cations and dications. The frequency of charge oscillation is determined by the energy separation of the X and A states; this separation can change depending on the direction and strength of the field. Furthermore, fields along the molecular axis can cause extensive mixing between the field-free X and A configurations. For coherent superpositions of the X and A states, the charge oscillations are characterized by two frequencies–the driving frequency of the laser field of the probe pulse and the intrinsic frequency due to the energy separation between the X and A states. For linear and circularly polarized pulses, the ionization rates show marked differences that depend on the polarization direction of the pulse, the carrier envelope phase and initial phase of the superposition. Varying the initial phase of the superposition at the beginning of the probe pulse is analogous to changing the delay between the pump and probe pulses. The charge oscillation in the coherent superposition of the X and A states results in maxima and minima in the ionization yield as a function of the superposition phase

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Distinct Bimetallic Cooperativity Among Water Reduction Catalysts Containing [Co III Co III ], [Ni II Ni II ], and [Zn II Zn II ] Cores

We report three binuclear species [LCo III 2 (μ-Pz) 2 ](ClO 4 ) 3 (1), [LNi II 2 (CH 3 OH) 2 Cl 2 ]ClO 4 (2), and [LZn II 2 Cl 2 ]PF 6 (3) supported by the deprotonated form of the ligand 2,6-bis[bis(2-pyridylmethyl) amino-methyl]-4-methylphenol were synthesized, structurally characterized as solids and in solution, and had their electrochemical and spectroscopic behavior established. Species 1–3 had their water reduction ability studied aiming to interrogate the possible cooperative catalytic activity between two neighboring metal centers. Species 1 and 2 reduced H 2 O to H 2 effectively at an applied potential of -1.6 V Ag/AgCl , yielding turnover numbers of 2,820 and 2,290, respectively, after 30 minutes. Species 3 lacked activity and was used as a negative control to eliminate the possibility of ligand-based catalysis. Pre- and post-catalytic data gave evidence of the molecular nature of the process within the timeframe of the experiments. Species 1 showed structural, rather than electronic cooperativity, while species 2 displayed no obvious cooperativity. DFT methods complemented the experimental results determining plausible mechanisms.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Ellipticity controlled dissociative double ionization of ethane by strong fields

We report the yields of all dissociation channels of ethane dications produced by strong field double ionization were measured. It was found that the branching ratios can be controlled by varying the ellipticity of laser pulses. The CH 3 + formation and H + formation channels show a clear competition, producing the highest and lowest branching ratios at ellipticity of ~0.6, respectively. With the help of theoretical calculations, such a control was attributed to the ellipticity dependent yields of different sequential ionization pathways.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Ab Initio Direct Dynamics

The reactivity and dynamics of molecular systems can be explored computationally by classical trajectory calculations. The traditional approach involves fitting a functional form of a potential energy surface (PES) to the energies from a large number of electronic structure calculations and then integrating numerous trajectories on this fitted PES to model the molecular dynamics. The ever-decreasing cost of computing and continuing advances in computational chemistry software have made it possible to use electronic structure calculations directly in molecular dynamics simulations without first having to construct a fitted PES. In this “on-the-fly” approach, every time the energy and its derivatives are needed for the integration of the equations of motion, they are obtained directly from quantum chemical calculations. This approach started to become practical in the mid-1990s as a result of increased availability of inexpensive computer resources and improved computational chemistry software. The application of direct dynamics calculations has grown rapidly over the last 25 years and would require a lengthy review article. The present Account is limited to some of our contributions to methods development and various applications. To improve the efficiency of direct dynamics calculations, we developed a Hessian-based predictor-corrector algorithm for integrating classical trajectories. Hessian updating made this even more efficient. Furthermore, this approach was also used to improve algorithms for following the steepest descent reaction paths. For larger molecular systems, we developed an extended Lagrangian approach in which the electronic structure is propagated along with the molecular structure. Strong field chemistry is a rapidly growing area, and to improve the accuracy of molecular dynamics in intense laser fields, we included the time-varying electric field in a novel predictor-corrector trajectory integration algorithm. Since intense laser fields can excite and ionize molecules, we extended our studies to include electron dynamics. Specifically, we developed code for time-dependent configuration interaction electron dynamics to simulate strong field ionization by intense laser pulses. Our initial application of ab initio direct dynamics in 1994 was to CH 2 O → H 2 + CO; the calculated vibrational distributions in the products were in very good agreement with experiment. In the intervening years, we have used direct dynamics to explore energy partitioning in various dissociation reactions, unimolecular dissociations yielding three fragments, reactions with branching after the transition state, nonstatistical dynamics of chemically activated molecules, dynamics of molecular fragmentation by intense infrared laser pulses, selective activation of specific dissociation channels by aligned intense infrared laser fields, angular dependence of strong field ionization, and simulation of sequential double ionization.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗