Source author record

Alberto Castro

Alberto Castro appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

17works
11topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

17 published item(s)

preprint2022arXiv

Floquet engineering the band structure of materials with optimal control theory

We demonstrate that the electronic structure of a material can be deformed into Floquet pseudo-bands with arbitrarily tailored shapes. We achieve this goal with a novel combination of quantum optimal control theory and Floquet engineering. The power and versatility of this framework is demonstrated here by utilizing the independent-electron tight-binding description of the $π$ electronic system of graphene. We show several prototype examples focusing on the region around the K (Dirac) point of the Brillouin zone: creation of a gap with opposing flat valence and conduction bands, creation of a gap with opposing concave symmetric valence and conduction bands -- which would correspond to a material with an effective negative electron-hole mass --, or closure of the gap when departing from a modified graphene model with a non-zero field-free gap. We employ time periodic drives with several frequency components and polarizations, in contrast to the usual monochromatic fields, and use control theory to find the amplitudes of each component that optimize the shape of the bands as desired. In addition, we use quantum control methods to find realistic switch-on pulses that bring the material into the predefined stationary Floquet band structure, i.e. into a state in which the desired Floquet modes of the target bands are fully occupied, so that they should remain stroboscopically stationary, with long lifetimes, when the weak periodic drives are started. Finally, we note that although we have focused on solid state materials, the technique that we propose could be equally used for the Floquet engineering of ultracold atoms in optical lattices, and to other non-equilibrium dynamical and correlated systems.

preprint2019arXiv

Octopus, a computational framework for exploring light-driven phenomena and quantum dynamics in extended and finite systems

Over the last years extraordinary advances in experimental and theoretical tools have allowed us to monitor and control matter at short time and atomic scales with a high-degree of precision. An appealing and challenging route towards engineering materials with tailored properties is to find ways to design or selectively manipulate materials, especially at the quantum level. To this end, having a state-of-the-art ab initio computer simulation tool that enables a reliable and accurate simulation of light-induced changes in the physical and chemical properties of complex systems is of utmost importance. The first principles real-space-based Octopus project was born with that idea in mind, providing an unique framework allowing to describe non-equilibrium phenomena in molecular complexes, low dimensional materials, and extended systems by accounting for electronic, ionic, and photon quantum mechanical effects within a generalized time-dependent density functional theory framework. The present article aims to present the new features that have been implemented over the last few years, including technical developments related to performance and massive parallelism. We also describe the major theoretical developments to address ultrafast light-driven processes, like the new theoretical framework of quantum electrodynamics density-functional formalism (QEDFT) for the description of novel light-matter hybrid states. Those advances, and other being released soon as part of the Octopus package, will enable the scientific community to simulate and characterize spatial and time-resolved spectroscopies, ultrafast phenomena in molecules and materials, and new emergent states of matter (QED-materials).

preprint2016arXiv

Tailored pump-probe transient spectroscopy with time-dependent density-functional theory: controlling absorption spectra

Recent advances in laser technology allow us to follow electronic motion at its natural time-scale with ultra-fast time resolution, leading the way towards attosecond physics experiments of extreme precision. In this work, we assess the use of tailored pumps in order to enhance (or reduce) some given features of the probe absorption (for example, absorption in the visible range of otherwise transparent samples). This type of manipulation of the system response could be helpful for its full characterization, since it would allow us to visualize transitions that are dark when using unshaped pulses. In order to investigate these possibilities, we perform first a theoretical analysis of the non-equilibrium response function in this context, aided by one simple numerical model of the Hydrogen atom. Then, we proceed to investigate the feasibility of using time-dependent density-functional theory as a means to implement, theoretically, this absorption-optimization idea, for more complex atoms or molecules. We conclude that the proposed idea could in principle be brought to the laboratory: tailored pump pulses can excite systems into light-absorbing states. However, we also highlight the severe numerical and theoretical difficulties posed by the problem: large-scale non-equilibrium quantum dynamics are cumbersome, even with TDDFT, and the shortcomings of state-of-the-art TDDFT functionals may still be serious for these out-of-equilibrium situations.

preprint2015arXiv

Real-space grids and the Octopus code as tools for the development of new simulation approaches for electronic systems

Real-space grids are a powerful alternative for the simulation of electronic systems. One of the main advantages of the approach is the flexibility and simplicity of working directly in real space where the different fields are discretized on a grid, combined with competitive numerical performance and great potential for parallelization. These properties constitute a great advantage at the time of implementing and testing new physical models. Based on our experience with the Octopus code, in this article we discuss how the real-space approach has allowed for the recent development of new ideas for the simulation of electronic systems. Among these applications are approaches to calculate response properties, modeling of photoemission, optimal control of quantum systems, simulation of plasmonic systems, and the exact solution of the Schrödinger equation for low-dimensionality systems.

preprint2015arXiv

Shaped electric fields for fast optimal manipulation of electron spin and position in a double quantum dot

We use quantum optimal control theory algorithms to design external electric fields that drive the coupled spin and orbital dynamics of an electron in a double quantum dot, subject to the spin-orbit interaction and Zeeman magnetic fields. We obtain time-profiles of multi-frequency electric pulses which increase the rate of spin-flip transitions by several orders of magnitude in comparison with monochromatic fields, where the spin Rabi oscillations were predicted to be very slow. This precise, with the accuracy higher than $10^{-4}$, and fast, at the timescale of the order of 0.1 ns, simultaneous control of the electron spin orientation and position is achieved while keeping the electron in the four lowest tunneling- and Zeeman-split levels through the duration of the pulse, facilitating high fidelity and suggesting useful applications in spintronics and quantum information devices.

preprint2014arXiv

Enhancing and controlling single-atom high-harmonic generation spectra: a time-dependent density-functional scheme

High harmonic generation (HHG) provides a flexible framework for the development of coherent light sources in the extreme-ultraviolet and soft x-ray regimes. However it suffers from low conversion efficiencies as the control of the HHG spectral and temporal characteristics requires manipulating electron trajectories on attosecond time scale. The phase matching mechanism has been employed to selectively enhance specific quantum paths leading to HHG. A few important fundamental questions remain open, among those how much of the enhancement can be achieved by the single-emitter and what is the role of correlations (or the electronic structure) in the selectivity and control of HHG generation. Here we address those questions by examining computationally the possibility of optimizing the HHG spectrum of isolated Hydrogen and Helium atoms by shaping the slowly varying envelope of a 800 nm, 200-cycles long laser pulse. The spectra are computed with a fully quantum mechanical description, by explicitly computing the time-dependent dipole moment of the systems using a first-principles time-dependent density-functional approach (exact for the case of H). The sought optimization corresponds to the selective enhancement of single harmonics, which we find to be significant. This selectivity is entirely due to the single atom response, and not due to any propagation or phase-matching effect. In fact, this single-emitter enhancement adds to the phase-matching techniques to achieving even larger HHG enhancement factors. Moreover, we see that the electronic correlation plays a role in the determining the degree of optimization that can be obtained.

preprint2013arXiv

Optimal control theory for quantum-classical systems: Ehrenfest Molecular Dynamics based on time-dependent density-functional theory

We derive the fundamental equations of an optimal control theory for systems containing both quantum electrons and classical ions. The system is modeled with Ehrenfest dynamics, a non-adiabatic variant of molecular dynamics. The general formulation, that needs the fully correlated many-electron wave function, can be simplified by making use of time-dependent density-functional theory. In this case, the optimal control equations require some modifications that we will provide. The abstract general formulation is complemented with the simple example of the H$_2^+$ molecule in the presence of a laser field.

preprint2013arXiv

Simulating pump-probe photo-electron and absorption spectroscopy on the attosecond time-scale with time-dependent density-functional theory

Molecular absorption and photo-electron spectra can be efficiently predicted with real-time time-dependent density-functional theory (TDDFT). We show here how these techniques can be easily extended to study time-resolved pump-probe experiments in which a system response (absorption or electron emission) to a probe pulse, is measured in an excited state. This simulation tool helps to interpret the fast evolving attosecond time-resolved spectroscopic experiments, where the electronic motion must be followed at its natural time-scale. We show how the extra degrees of freedom (pump pulse duration, intensity, frequency, and time-delay), which are absent in a conventional steady state experiment, provide additional information about electronic structure and dynamics that improve a system characterization. As an extension of this approach, time-dependent 2D spectroscopies can also be simulated, in principle, for large-scale structures and extended systems.

preprint2011arXiv

On the Combination of TDDFT with Molecular Dynamics: New Developments

In principle, we should not need the time-dependent extension of density-functional theory (TDDFT) for excitations, and in particular not for Molecular Dynamics (MD) studies: the theorem by Hohenberg and Kohn teaches us that for any observable that we wish to look at (including dynamical properties or observables dependent on excited states) there is a corresponding functional of the ground-state density. Yet the unavailability of such magic functionals in many cases (the theorem is a non-constructive existence result) demands the development and use of the alternative exact reformulation of quantum mechanics provided by TDDFT. This theory defines a convenient route to electronic excitations and to the dynamics of a many-electron system subject to an arbitrary time-dependent perturbation. This is, in fact, the main purpose of inscribing TDDFT in a MD framework -the inclusion of the effect of electronic excited states in the dynamics. However, as we will show in this review, it may not be the only use of TDDFT in this context. In this manuscript, we review two recent proposals: In Section 1.2, we show how TDDFT can be used to design efficient gsBOMD algorithms -even if the electronic excited states are in this case not relevant. The work described in Section 1.3 addresses the problem of mixed quantum-classical systems at thermal equilibrium.

preprint2011arXiv

Optimal Control of the Electronic Current Density - An application to one- and two-dimensional one-electron systems

Quantum optimal control theory is a powerful tool for engineering quantum systems subject to external fields such as the ones created by intense lasers. The formulation relies on a suitable definition for a target functional, that translates the intended physical objective to a mathematical form. We propose the use of target functionals defined in terms of the one-particle density and its current. A strong motivation for this is the possibility of using time-dependent density-functional theory for the description of the system dynamics. We exemplify this idea by defining an objective functional that on one hand attempts a large overlap with a target density and on the other hand minimizes the current. The latter requirement leads to optimized states with increased stability, which we prove with a few examples of one- and two-dimensional one-electron systems.

preprint2011arXiv

Optimization Schemes for Selective Molecular Cleavage with Tailored Ultrashort Laser Pulses

We present some approaches to the computation of ultra-fast laser pulses capable of selectively breaking molecular bonds. The calculations are based on a mixed quantum-classical description: The electrons are treated quantum mechanically (making use of time-dependent density-functional theory), whereas the nuclei are treated classically. The temporal shape of the pulses is tailored to maximise a control target functional which is designed to produce the desired molecular cleavage. The precise definition of this functional is a crucial ingredient: we explore expressions based on the forces, on the momenta and on the velocities of the nuclei. The algorithm used to find the optimum pulse is also relevant; we test both direct gradient-free algorithms, as well as schemes based on formal optimal control theory. The tests are performed both on one dimensional models of atomic chains, and on first-principles descriptions of molecules.

preprint2011arXiv

Quantum Optimal Control Theory in the Linear Response Formalism

Quantum optimal control theory (QOCT) aims at finding an external field that drives a quantum system in such a way that optimally achieves some predefined target. In practice this normally means optimizing the value of some observable, a so called merit function. In consequence, a key part of the theory is a set of equations, which provides the gradient of the merit function with respect to parameters that control the shape of the driving field. We show that these equations can be straightforwardly derived using the standard linear response theory, only requiring a minor generalization -- the unperturbed Hamiltonian is allowed to be time-dependent. As a result, the aforementioned gradients are identified with certain response functions. This identification leads to a natural reformulation of QOCT in term of the Keldysh contour formalism of the quantum many-body theory. In particular, the gradients of the merit function can be calculated using the diagrammatic technique for non-equilibrium Green's functions, which should be helpful in the application of QOCT to computationally difficult many-electron problems.

preprint2011arXiv

Scattering of a proton with the Li4 cluster: non-adiabatic molecular dynamics description based on time-dependent density-functional theory

We have employed non-adiabatic molecular dynamics based on time-dependent density-functional theory to characterize the scattering behaviour of a proton with the Li$_4$ cluster. This technique assumes a classical approximation for the nuclei, effectively coupled to the quantum electronic system. This time-dependent theoretical framework accounts, by construction, for possible charge transfer and ionization processes, as well as electronic excitations, which may play a role in the non-adiabatic regime. We have varied the incidence angles in order to analyze the possible reaction patterns. The initial proton kinetic energy of 10 eV is sufficiently high to induce non-adiabatic effects. For all the incidence angles considered the proton is scattered away, except in one interesting case in which one of the Lithium atoms captures it, forming a LiH molecule. This theoretical formalism proves to be a powerful, effective and predictive tool for the analysis of non-adiabatic processes at the nanoscale.

preprint2011arXiv

Water Molecules in Ultrashort Intense Laser Fields

Ionization and excitation of water molecules in intense laser pulses is studied theoretically by solving the three-dimensional time-dependent electronic Schrödinger equation within the single-active-electron approximation. The possibility to image orbital densities by measurement of the orientation-dependent ionization of H2O in few-cycle, 800nm linear-polarized laser pulses is investigated. While the highest-occupied molecular orbital 1b1 is found to dominate the overall ionization behavior, contributions from the energetically lower lying 3a1 orbital dominate the ionization yield in the nodal plane of the 1b1 orbital. The ratio of the ionization yields of the two orbitals depends on the intensity. Furthermore, even for laser pulses as long as 8 cycles the orientation-dependent ion yield depends on the carrier-envelope phase. In the interpretation of the orientation-dependent ionization as an imaging tool these effects have to be considered.

preprint2010arXiv

Alignment-Dependent Ionization of N$_2$, O$_2$, and CO$_2$ in Intense Laser Fields

The ionization probability of N$_2$, O$_2$, and CO$_2$ in intense laser fields is studied theoretically as a function of the alignment angle by solving the time-dependent Schrödinger equation numerically assuming only the single-active-electron approximation. The results are compared to recent experimental data [D.~Pavi{č}i{ć} et al., Phys.\,Rev.\,Lett.\ {\bf 98}, 243001 (2007)] and good agreement is found for N$_2$ and O$_2$. For CO$_2$ a possible explanation is provided for the failure of simplified single-active-electron models to reproduce the experimentally observed narrow ionization distribution. It is based on a field-induced coherent core-trapping effect.

preprint2010arXiv

Controlling the Dynamics of Many-Electron Systems from First Principles: A Marriage of Optimal Control and Time-Dependent Density-Functional Theory

Quantum Optimal Control Theory (QOCT) provides the necessary tools to theoretically design driving fields capable of controlling a quantum system towards a given state or along a prescribed path in Hilbert space. This theory must be complemented with a suitable model for describing the dynamics of the quantum system. Here, we are concerned with many electron systems (atoms, molecules, quantum dots, etc) irradiated with laser pulses. The full solution of the many electron Schr{ö}dinger equation is not feasible in general, and therefore, if we aim to an ab initio description, a suitable choice is time-dependent density-functional theory (TDDFT). In this work, we establish the equations that combine TDDFT with QOCT, and demonstrate their numerical feasibility with examples.

preprint2009arXiv

A modified Ehrenfest formalism for efficient large-scale ab initio molecular dynamics

We present in detail the recently derived ab-initio molecular dynamics (AIMD) formalism [Phys. Rev. Lett. 101 096403 (2008)], which due to its numerical properties, is ideal for simulating the dynamics of systems containing thousands of atoms. A major drawback of traditional AIMD methods is the necessity to enforce the orthogonalization of the wave-functions, which can become the bottleneck for very large systems. Alternatively, one can handle the electron-ion dynamics within the Ehrenfest scheme where no explicit orthogonalization is necessary, however the time step is too small for practical applications. Here we preserve the desirable properties of Ehrenfest in a new scheme that allows for a considerable increase of the time step while keeping the system close to the Born-Oppenheimer surface. We show that the automatically enforced orthogonalization is of fundamental importance for large systems because not only it improves the scaling of the approach with the system size but it also allows for an additional very efficient parallelization level. In this work we provide the formal details of the new method, describe its implementation and present some applications to some test systems. Comparisons with the widely used Car-Parrinello molecular dynamics method are made, showing that the new approach is advantageous above a certain number of atoms in the system. The method is not tied to a particular wave-function representation, making it suitable for inclusion in any AIMD software package.