Three parallel algorithms for classical molecular dynamics are presented. The first assigns each processor a fixed subset of atoms; the second assigns each a fixed subset of inter-atomic forces to compute; the third assigns each a fixed spatial region. The algorithms are suitable for molecular dynamics models which can be difficult to parallelize efficiently—those with short-range forces where the neighbors of each atom change rapidly. They can be implemented on any distributed-memory parallel machine which allows for message-passing of data between independently executing processors. The algorithms are tested on a standard Lennard-Jones benchmark problem for system sizes ranging from 500 to 100,000,000 atoms on several parallel supercomputers--the nCUBE 2, Intel iPSC/860 and Paragon, and Cray T3D. Comparing the results to the fastest reported vectorized Cray Y-MP and C90 algorithm shows that the current generation of parallel machines is competitive with conventional vector supercomputers even for small problems. For large problems, the spatial algorithm achieves parallel efficiencies of 90% and a 1840-node Intel Paragon performs up to 165 faster than a single Cray C9O processor. Trade-offs between the three algorithms and guidelines for adapting them to more complex molecular dynamics simulations are also discussed.

Fast parallel algorithms for short-range molecular dynamics

Electricity-driven water splitting can facilitate the storage of electrical energy in the form of hydrogen gas. As a half-reaction of electricity-driven water splitting, the oxygen evolution reaction (OER) is the major bottleneck due to the sluggish kinetics of this four-electron transfer reaction. Developing low-cost and robust OER catalysts is critical to solving this efficiency problem in water splitting. The catalyst design has to be built based on the fundamental understanding of the OER mechanism and the origin of the reaction overpotential. In this article, we summarize the recent progress in understanding OER mechanisms, which include the conventional adsorbate evolution mechanism (AEM) and lattice-oxygen-mediated mechanism (LOM) from both theoretical and experimental aspects. We start with the discussion on the AEM and its linked scaling relations among various reaction intermediates. The strategies to reduce overpotential based on the AEM and its derived descriptors are then introduced. To further reduce the OER overpotential, it is necessary to break the scaling relation of HOO* and HO* intermediates in conventional AEM to go beyond the activity limitation of the volcano relationship. Strategies such as stabilization of HOO*, proton acceptor functionality, and switching the OER pathway to LOM are discussed. The remaining questions on the OER and related perspectives are also presented at the end.

/pdf/a-review-on-fundamentals-for-designing-oxygen-evolution-3mw3rcvtbe.pdf

A review on fundamentals for designing oxygen evolution electrocatalysts

The Journal of Physical Chemistry Letters

Summary The vision of a hydrogen economy relies on efficient utilization and production of hydrogen in a hydrogen fuel cell and a water electrolyzer. In both technologies, the sluggish kinetics of oxygen reduction reaction (ORR) and oxygen evolution reaction (OER) account for the most efficiency loss because the reactions on catalytic sites are constrained by adsorption-energy scaling relations involving intermediates of *OOH, *O, and *OH (where * denotes the active site). Therefore, a novel paradigm for catalyst design is required by overcoming or circumventing the adsorption-energy scaling relations. In this Review, the fundamentals of oxygen electrocatalysis, including reaction of the mechanism and origin of the adsorption-energy scaling relationship, are first introduced. Crucial strategies to overcome the scaling relations are then summarized. Finally, future research directions in this area are proposed. This work provides guidelines for the rational design of efficient catalysts for oxygen electrocatalysis beyond the limitation posed by scaling relations.

Strategies to Break the Scaling Relation toward Enhanced Oxygen Electrocatalysis

DOI: 10.1002/adom.201900681 reason for this interest relies on the fact that THz radiation can couple resonantly to numerous fundamental motions of ions, electrons, and electron spins in all phases of matter. For example, in solids, the THz range overlaps with the frequency of lattice vibrations (phonons), the collision rates of conduction electrons, the binding energy of bound electron–hole pairs (excitons), and the precession frequency of spin waves (magnons). Consequently, THz radiation, both continuous-wave and pulsed, has been used for characterization of and gaining insight into elementary processes in complex materials. The majority of these studies used relatively weak THz fields and, thus, probed the linear response of the material, without inducing notable material modifications. Only recently, however, completely new avenues in THz science were opened up by triggering nonlinear THz responses of materials.[3–13] Instead of using weak fields to primarily observe selected THz modes such as phonons or magnons, strong fields allow one to actively drive them to unprecedentedly large amplitudes, potentially thereby resulting in novel states of matter.[6,8,11] For example, simulations suggest that exciting matter with intense THz transients may lead to massive modifications of electrically[14] or magnetically[15] ordered domains and enable the acceleration of free ions to ≈1 MeV,[16] and postacceleration to 50–100 MeV energies.[17] Remarkable experimental results such as switching of magnetic order,[18,19] parametric amplification of optical phonons,[20] novel insights into spin-lattice coupling,[21,22] and acceleration of free electrons in a THz linear accelerator[23] were achieved only recently. This progress has been made possible by the development of laser-driven table-top THz sources routinely providing pulses with unprecedented energies and peak electric and magnetic field strengths throughout the entire THz spectral range. Different laser-based THz pulse generation techniques can be used to access different parts of the spectral range extending from 0.1 to 10 THz. Some of the recently developed technologies enable the generation of radiation with even larger bandwidth or tuning range up to 100 THz and beyond, which lead to an extension of what is called the THz spectral range. An overview of the approximate spectral coverage and the achieved highest pulse energies and peak electric-field strengths of various laserdriven technologies is given in Figure 1.

/pdf/laser-driven-strong-field-terahertz-sources-36tf7wv2ib.pdf

Laser‐Driven Strong‐Field Terahertz Sources

"Hidden phases" are metastable collective states of matter that are typically not accessible on equilibrium phase diagrams. These phases can host exotic properties in otherwise conventional materials and hence may enable novel functionality and applications, but their discovery and access are still in early stages. Using intense terahertz electric field excitation, we found that an ultrafast phase transition into a hidden ferroelectric phase can be dynamically induced in quantum paraelectric strontium titanate (SrTiO3). The induced lowering in crystal symmetry yields substantial changes in the phonon excitation spectra. Our results demonstrate collective coherent control over material structure, in which a single-cycle field drives ions along the microscopic pathway leading directly to their locations in a new crystalline phase on an ultrafast time scale.

Terahertz field–induced ferroelectricity in quantum paraelectric SrTiO3

"Hidden phases" are metastable collective states of matter that are typically not accessible on equilibrium phase diagrams. These phases can host exotic properties in otherwise conventional materials and hence may enable novel functionality and applications, but their discovery and access are still in early stages. Using intense terahertz electric field excitation, we show that an ultrafast phase transition into a hidden ferroelectric phase can be dynamically induced in the quantum paraelectric SrTiO$_3$. The induced lowering in crystal symmetry yields dramatic changes in the phonon excitation spectra. Our results demonstrate collective coherent control over material structure, in which a single-cycle field drives ions along the microscopic pathway leading directly to their locations in a new crystalline phase on an ultrafast timescale.

Terahertz-Field-Induced Ferroelectricity in Quantum Paraelectric SrTiO$_3$

We report the optical conductivity in high-quality crystals of the chiral topological semimetal CoSi, which hosts exotic quasiparticles known as multifold fermions. We find that the optical response is separated into several distinct regions as a function of frequency, each dominated by different types of quasiparticles. The low-frequency intraband response is captured by a narrow Drude peak from a high-mobility electron pocket of double Weyl quasiparticles, and the temperature dependence of the spectral weight is consistent with its Fermi velocity. By subtracting the low-frequency sharp Drude and phonon peaks at low temperatures, we reveal two intermediate quasilinear interband contributions separated by a kink at 0.2 eV. Using Wannier tight-binding models based on first-principle calculations, we link the optical conductivity above and below 0.2 eV to interband transitions near the double Weyl fermion and a threefold fermion, respectively. We analyze and determine the chemical potential relative to the energy of the threefold fermion, revealing the importance of transitions between a linearly dispersing band and a flat band. More strikingly, below 0.1 eV our data are best explained if spin-orbit coupling is included, suggesting that at these energies, the optical response is governed by transitions between a previously unobserved fourfold spin-3/2 node and a Weyl node. Our comprehensive combined experimental and theoretical study provides a way to resolve different types of multifold fermions in CoSi at different energy. More broadly, our results provide the necessary basis to interpret the burgeoning set of optical and transport experiments in chiral topological semimetals.

https://www.pnas.org/content/pnas/117/44/27104.full.pdf

Optical signatures of multifold fermions in the chiral topological semimetal CoSi.

The properties of a material are often strongly influenced by its surfaces. Depending on the nature of the chemical bonding in a material, its surface can undergo a variety of stabilizing reconstructions that dramatically alter the chemical reactivity, light absorption, and electronic band offsets. For decades, ab initio thermodynamics has been the method of choice for computationally determining the surface phase diagram of a material under different conditions. The surfaces considered for these studies, however, are often human-selected and too few in number, leading both to insufficient exploration of all possible surfaces and to biases toward portions of the composition–structure phase space that often do not encompass the most stable surfaces. To overcome these limitations and automate the discovery of realistic surfaces, we combine density functional theory and grand canonical Monte Carlo (GCMC) into “ab initio GCMC.” This paper presents the successful application of ab initio GCMC to the study of o...

Automatic Prediction of Surface Phase Diagrams Using Ab Initio Grand Canonical Monte Carlo

Experimental results have shown promising catalytic activity for the oxygen evolution reaction (OER) on the perovskite-type material CaMnO3. Through density functional theory investigations, we study the OER mechanism on CaMnO3, on the basis of a thermodynamic stability approach. Our results reveal that the formation of Mn vacancies caused by the solubility of Mn enhances lattice oxygen activity, which then reduces the energy of the adsorbate *OOH and therefore loosens the lower overpotential limit predicted from the “scaling relationship”. This effect suggests that, by doping manganite oxides with soluble elements, we could enhance their OER catalytic activities due to the high surface lattice oxygen activity induced by surface metal ion vacancies.

Tian Qiu

Papers

Terahertz field–induced ferroelectricity in quantum paraelectric SrTiO3

Terahertz-Field-Induced Ferroelectricity in Quantum Paraelectric SrTiO$_3$

Optical signatures of multifold fermions in the chiral topological semimetal CoSi.

Automatic Prediction of Surface Phase Diagrams Using Ab Initio Grand Canonical Monte Carlo

Ab initio Simulation Explains the Enhancement of Catalytic Oxygen Evolution on CaMnO3