Publications

Dynamic manipulation of magnetism in topological materials is demonstrated here via a Floquet engineering approach using circularly polarized light. Increasing the strength of the laser field, besides the expected topological phase transition, the magnetically doped topological insulator thin film also undergoes a magnetic phase transition from ferromagnetism to paramagnetism, whose critical behavior strongly depends on the quantum quenching. In sharp contrast to the equilibrium case, the non-equilibrium Curie temperatures vary for different time scale and experimental setup, not all relying on change of topology. Our discoveries deepen the understanding of the relationship between topology and magnetism in the non-equilibrium regime and extend optoelectronic device applications to topological materials.

Floquet engineering consists in the modification of physical systems by the application of periodic time-dependent perturbations. The search for the shape of the periodic perturbation that best modifies the properties of a system in order to achieve some predefined metastable target behavior can be formulated as an optimal control problem. We discuss several ways to formulate and solve this problem. We present, as examples, some applications in the context of material science, although the methods discussed here are valid for any quantum system (from molecules and nanostructures to extended periodic and non periodic quantum materials). In particular, we show how one can achieve the manipulation of the Floquet pseudo-bandstructure of a transition metal dichalcogenide monolayer (MoS2).

The coupling of light to electrical charge carriers in semiconductors is the foundation of many technological applications. Attosecond transient absorption spectroscopy measures simultaneously how excited electrons and the vacancies they leave behind dynamically react to the applied optical fields. In compound semiconductors, these dynamics can be probed via any of their atomic constituents. Often, the atomic species forming the compound contribute comparably to the relevant electronic properties of the material. One therefore expects to observe similar dynamics, irrespective of the choice of atomic species via which it is probed. Here, we show in the two-dimensional transition metal dichalcogenide semiconductor MoSe

Quantum computing has emerged as a promising platform for simulating strongly correlated systems in chemistry, for which the standard quantum chemistry methods are either qualitatively inaccurate or too expensive. However, due to the hardware limitations of the available noisy near-term quantum devices, their application is currently limited only to small chemical systems. One way for extending the range of applicability can be achieved within the quantum embedding approach. Herein, we employ the projection-based embedding method for combining the variational quantum eigensolver (VQE) algorithm, although not limited to, with density functional theory (DFT). The developed VQE-in-DFT method is then implemented efficiently on a real quantum device and employed for simulating the triple bond breaking process in butyronitrile. The results presented herein show that the developed method is a promising approach for simulating systems with a strongly correlated fragment on a quantum computer. The developments as well as the accompanying implementation will benefit many different chemical areas including the computer aided drug design as well as the study of metalloenzymes with a strongly correlated fragment.

We present a semilocal exchange-correlation energy functional for noncollinear spin density functional theory based on short-range expansions of the spin-resolved exchange hole and the two-body density matrix. Our functional is explicitly derived for noncollinear magnetism, U(1) and SU(2) gauge invariant, and gives rise to nonvanishing exchange-correlation torques. Testing the functional for frustrated antiferromagnetic chromium clusters, the exchange part is shown to perform favorably compared to the far more expensive Slater and optimized effective potentials, and a delicate interplay between exchange and correlation torques is uncovered.

The excitonic insulator is an electronically driven phase of matter that emerges upon the spontaneous formation and Bose condensation of excitons. Detecting this exotic order in candidate materials is a subject of paramount importance, as the size of the excitonic gap in the band structure establishes the potential of this collective state for superfluid energy transport. However, the identification of this phase in real solids is hindered by the coexistence of a structural order parameter with the same symmetry as the excitonic order. Only a few materials are currently believed to host a dominant excitonic phase, Ta2NiSe5 being the most promising. Here, we test this scenario by using an ultrashort laser pulse to quench the broken-symmetry phase of this transition metal chalcogenide. Tracking the dynamics of the material's electronic and crystal structure after light excitation reveals spectroscopic fingerprints that are compatible only with a primary order parameter of phononic nature. We rationalize our findings through state-of-the-art calculations, confirming that the structural order accounts for most of the gap opening. Our results suggest that the spontaneous symmetry breaking in Ta2NiSe5 is mostly of structural character, hampering the possibility to realize quasi-dissipationless energy transport.

Living systems are intrinsically nonequilibrium: They use metabolically derived chemical energy to power their emergent dynamics and self-organization. A crucial driver of these dynamics is the cellular cytoskeleton, a defining example of an active material where the energy injected by molecular motors cascades across length scales, allowing the material to break the constraints of thermodynamic equilibrium and display emergent nonequilibrium dynamics only possible due to the constant influx of energy. Notwithstanding recent experimental advances in the use of local probes to quantify entropy production and the breaking of detailed balance, little is known about the energetics of active materials or how energy propagates from the molecular to emergent length scales. Here, we use a recently developed picowatt calorimeter to experimentally measure the energetics of an active microtubule gel that displays emergent large-scale flows. We find that only approximately one-billionth of the system’s total energy consumption contributes to these emergent flows. We develop a chemical kinetics model that quantitatively captures how the system’s total thermal dissipation varies with ATP and microtubule concentrations but that breaks down at high motor concentration, signaling an interference between motors. Finally, we estimate how energy losses accumulate across scales. Taken together, these results highlight energetic efficiency as a key consideration for the engineering of active materials and are a powerful step toward developing a nonequilibrium thermodynamics of living systems.

Lipid molecules such as cholesterol interact with the surface of integral membrane proteins (IMP) in a mode different from drug-like molecules in a protein binding pocket. These differences are due to the lipid molecule’s shape, the membrane’s hydrophobic environment, and the lipid’s orientation in the membrane. We can use the recent increase in experimental structures in complex with cholesterol to understand protein-cholesterol interactions. We developed the RosettaCholesterol protocol consisting of (1) a prediction phase using an energy grid to sample and score native-like binding poses and (2) a specificity filter to calculate the likelihood that a cholesterol interaction site may be specific. We used a multi-pronged benchmark (self-dock, flip-dock, cross-dock, and global-dock) of protein-cholesterol complexes to validate our method. RosettaCholesterol improved sampling and scoring of native poses over the standard RosettaLigand baseline method in 91% of cases and performs better regardless of benchmark complexity. On the β2AR, our method found one likely-specific site, which is described in the literature. The RosettaCholesterol protocol quantifies cholesterol binding site specificity. Our approach provides a starting point for high-throughput modeling and prediction of cholesterol binding sites for further experimental validation.

Transmembrane helix folding and self-association play important roles in biological signaling and transportation pathways across biomembranes. With molecular simulations, studies to explore the structural biochemistry of this process have been limited to focusing on individual fragments of this process – either helix formation or dimerization. While at an atomistic resolution, it can be prohibitive to access long spatio-temporal scales, at the coarse grained (CG) level, current methods either employ additional constraints to prevent spontaneous unfolding or have a low resolution on sidechain beads that restricts the study of dimer disruption caused by mutations. To address these research gaps, in this work, we apply our recent, in-house developed CG model (ProMPT) to study the folding and dimerization of Glycophorin A (GpA) and its mutants in the presence of Dodecyl-phosphocholine (DPC) micelles. Our results first validate the two-stage model that folding and dimerization are independent events for transmembrane helices and found a positive correlation between helix folding and DPC-peptide contacts. The wild type (WT) GpA is observed to be a right-handed dimer with specific GxxxG contacts, which agrees with experimental findings. Specific point mutations reveal several features responsible for the structural stability of GpA. While the T87L mutant forms anti-parallel dimers due to an absence of T87 interhelical hydrogen bonds, a slight loss in helicity and a hinge-like feature at the GxxxG region develops for the G79L mutant. We note that the local changes in the hydrophobic environment, affected by the point mutation, contribute to the development of this helical bend. This work presents a holistic overview of the structural stability of GpA in a micellar environment, while taking secondary structural fluctuations into account. Moreover, it presents opportunities for applications of computationally efficient CG models to study conformational alterations of transmembrane proteins that have physiological relevance.

Many active particles, such as swimming micro-organisms or motor proteins, do work on their environment by going though a periodic sequence of shapes. Interactions between particles can lead to synchronization of their duty cycles. Here, we study the collective dynamics of a suspension of active particles coupled through hydrodynamics. We find that at high enough density the system transitions to a state of collective motion by a mechanism that is distinct from other instabilities in active matter systems. Second, we demonstrate that the emergent nonequilibrium states feature stationary chimera patterns in which synchronized and phase-isotropic regions coexist. Third, we show that in confinement, oscillatory flows and robust unidirectional pumping states exist, and can be selected by choice of alignment boundary conditions. These results point toward a new route to collective motion and pattern formation and could guide the design of new active materials.