First Principles Simulations and Modeling

Computer modeling has grown over the years into a scientific discipline on its own. Models are utilized to assess real-world phenomena maybe too complex to be analyzed in the laboratory or under hypotheses at a fraction of the cost of undertaking the actual activities. Models in industry, government, and educational institutions shorten design cycles, reduce costs, and enhance knowledge.

First Principles Simulations and Modeling Figure 1
Figure 1: Different types of chiral systems such as those in the figure may exhibit the CISS effect. From ACS Nano, 17, 6452 (2023).

In Physics, the modeling of materials, through what is known as “first-principles”, has become a major research field. By “first-principles” one understands the use of the fundamental quantum mechanical laws of nature without any assumptions. The properties of the materials should emerge from the numerical solution of these laws. The models here are in fact a faithful representation of reality, but in a controlled environment.

Researchers at IFIMAC have a long-standing and well-deserved international reputation on computer modeling and, in particular, on the development of efficient first-principles codes for the simulation of the optoelectronic and structural properties of molecules, materials, and systems in general described down to the atomic level (SIESTA, ANT. Gaussian, FIREBALL, MOLCAS, EDUS, etc.). All the other research lines at IFIMAC benefit one way or another from this expertise. This knowledge not only benefits fundamental research, but can also be transferred directly into the society through spin-offs.

IFIMAC researchers exploit these simulation codes to explore a wide variety of problems, including electronic transport, catalysis on reducible oxides and out-of-equilibrium electron dynamics.

Spin-dependent electronic transport is an essential feature not only in engineered devices for spintronics, but also in chemical and biological processes involving the propagation of a current through molecules. In recent years, it has been experimentally verified that an initially unpolarized beam of electrons will emerge polarized (in some cases, significantly) upon traversing a chiral molecule such as a DNA-like helicene; a realization of the so-called chirality-induced spin selectivity (CISS) effect. In a recent article published in ACS Nano IFIMAC researchers show based on the use of representation theory within the scattering formalism for transport, that the appearance of such spin polarization is fundamentally allowed in a much wider family of systems, namely those that lack symmetry planes or axes containing the propagation direction (Figure 1). The role of the contacts is hence as qualitatively important as that of the molecules, to the point that the presence of the latter is not generally needed to observe spin polarization. These predictions are illustrated by DFT calculations where they also show that the polarization vector is accompanied by a net spin accumulation in the system for finite bias, which may be detected in magneto-conductance setups.

First Principles Simulations and Modeling Figure 2
Figure 2: An electrically driven switch can be formed by a cytochrome C monolayer, sandwiched between a gold layer and a gold nanowire. From Angewandte Chemie, 58, 11923 (2019).

Proteins are key biological molecules that are responsible for numerous energy conversion processes such as photosynthesis or respiration. In recent years, proteins have been investigated in a new setting, namely in solid-state electronic junctions, with the goal of understanding the charge transfer mechanisms in these biomolecules, but also with the hope of developing a new generation of bio-inspired nanoscale electronic devices. A new step towards this goal was reported in Angewandte Chemie by a collaboration between Weizmann Institute of Science (Israel) and IFIMAC researchers. In this work, the authors show that a redox protein, cytochrome C, can behave as an electrically driven switch when incorporated in a solid-state junction with gold electrodes (Figure 2). By changing the external bias voltage in the junction, it was shown that the relevant molecular orbitals of the protein can be brought in and out of resonance with the chemical potential of the electrodes, which leads to the current-switch behavior. Showing transition from off- to on- resonance can be very challenging and this is the first time it has been achieved for proteins within the same working junction. Extensive ab initio DFT calculations revealed that the charge transport proceeds through the heme unit in these proteins and that the coupling between the protein’s frontier orbitals and the electrodes is sufficiently weak to prevent Fermi level pinning. The on-off change in the electrical current was shown to persist up to room temperature, demonstrating reversible, bias-controlled switching of a protein ensemble, which provides a realistic path to protein-based bioelectronics.

First Principles Simulations and Modeling Figure 3
Figure 3: Front cover artistic view of the catalytic process for acetylene hydrogenation. From J. Phys. Chem. C, 118, 5352 (2014).

The unique catalytic properties of ceria for the partial hydrogenation of alkynes have been examined for acetylene hydrogenation (Figure 3). Catalytic tests over polycrystalline CeO2 at different temperatures and H2/C2H2 ratios reveal ethylene selectivities in the range of 75–85% at high degrees of acetylene conversion and hint at the crucial role of hydrogen dissociation on the overall process. DFT is applied to CeO2(111) in order to investigate reaction intermediates and to calculate the enthalpy and energy barrier for each elementary step. At a high hydrogen coverage, b-C2H2 radicals adsorbed on-top of surface oxygen atoms are the initial reactive species forming C2H3 species effectively barrierless. The high alkene selectivity is owed to the lower activation barrier for subsequent hydrogenation leading to gas-phase C2H4 compared to that for the formation of b-C2H4 radical species. These findings rationalize for the first time the applicability of CeO2 as a catalyst for olefin production and potentially broaden its use for the hydrogenation of polyunsaturated and polyfunctionalized substrates containing triple bonds.

First Principles Simulations and Modeling Figure 4
Figure 4: Schematics of a pump-probe experiment where a first laser pulse excites the material while a second one explores the time evolution of the excitation. From Journal of Chemical Theory and Computation,19, 333 (2023).

In a recent manuscript published in Journal of Chemical Theory and Computation a collaboration of IFIMAC researchers have presented a theoretical framework and its numerical implementation to simulate the out-of-equilibrium electron dynamics induced by the interaction of ultrashort laser pulses in condensed-matter systems (Figure 4). Their approach is based on evolving in real time the density matrix of the system in reciprocal space. It considers excitonic and nonperturbative light–matter interactions. They show some relevant examples that illustrate the efficiency and flexibility of the approach to describe realistic ultrafast spectroscopy experiments. The approach is suitable for modeling the promising and emerging ultrafast studies at the attosecond time scale that aim at capturing the electron dynamics and the dynamical electron–electron correlations via X-ray absorption spectroscopy.