# Publications

The collective excitation of the conduction electrons in subwavelength structures gives rise to the Localized Surface Plasmon(LSP). The system consisting of two such LSPs, known as the dimer system,is of fundamental interest and is being actively investigated in the literature. Three regimes have been previously identified and they are the photonic regime, the strong coupling regime and the quantum tunneling regime. In this Letter, we propose a new regime for this intriguing systems, the intermediate regime. In this new regime, the quasistatic approximation, which is widely used to study such LSP systems, fails to capture the main physics: the multiple scattering of the electromagnetic waves between the two LSPs, which significantly modifies the properties of the resonant modes in the system. This intermediate regime provides a new route to explore in plasmonics, where controlling both the excited plasmon modes and the damping rates are of paramount significance.

The ability to control electromagnetic fields on the subwavelength scale could open exciting new venues in many fields of science. Transformation optics provides one way to attain such control through the local variation of the permittivity and permeability of a material. Here, we demonstrate another way to shape electromagnetic fields, taking advantage of the enormous size of the configuration space in combinatorial problems and the resonant scattering properties of metallic nanoparticles. Our design does not require the engineering of a material's electromagnetic properties and has relevance to the design of more flexible platforms for probing light-matter interaction and many body physics.

Despite several experiments on chaotic quantum transport in two-dimensional systems such as semiconductor quantum dots, corresponding quantum simulations within a real-space model have been out of reach so far. Here we carry out quantum transport calculations in real space and real time for a two-dimensional stadium cavity that shows chaotic dynamics. By applying a large set of magnetic fields we obtain a complete picture of magnetoconductance that indicates fractal scaling. In the calculations of the fractality we use detrended fluctuation analysis—a widely used method in time-series analysis—and show its usefulness in the interpretation of the conductance curves. Comparison with a standard method to extract the fractal dimension leads to consistent results that in turn qualitatively agree with the previous experimental data.

Increasing fidelity is the ultimate challenge of quantum information technology. In addition to decoherence and dissipation, fidelity is affected by internal imperfections such as impurities in the system. Here we show that the quality of quantum revival, i.e., periodic recurrence in the time evolution, can be restored almost completely by coupling the distorted system to an external field obtained from quantum optimal control theory. We demonstrate the procedure with wave-packet calculations in both one- and two-dimensional quantum wells, and analyze the required physical characteristics of the control field. Our results generally show that the inherent dynamics of a quantum system can be idealized at an extremely low cost.

An interpretation of the probability flux is given, based on a derivation of its eigenstates and relating them to coherent-state projections on a quantum wave function. An extended definition of the flux operator is obtained using coherent states. We present a "processed Husimi" representation, which makes decisions using many Husimi projections at each location. The processed Husimi representation reverse engineers or deconstructs the wave function, yielding the underlying classical ray structure. Our approach makes possible interpreting the dynamics of systems where the probability flux is uniformly zero or strongly misleading. The new technique is demonstrated by the calculation of particle flow maps of the classical dynamics underlying a quantum wave function in simple model systems such as a circular billiard with and without a magnetic field.

We investigate the transport of electrons in disordered and pristine graphene devices. Fano shot noise, a standard metric to assess the mechanism for electronic transport in mesoscopic devices, has been shown to produce almost the same magnitude (≈1/3) in ballistic and diffusive graphene devices and is therefore of limited applicability. We consider a two-terminal geometry where the graphene flake is contacted by narrow metallic leads. We propose that the dependence of the conductance on the position of one of the leads, a conductance profile, can give us insight into the charge flow, which can in turn be used to analyze the transport mechanism. Moreover, we simulate scanning probe microscopy (SPM) measurements for the same devices, which can visualize the flow of charge inside the device, thus complementing the transport calculations. From our simulations, we find that both the conductance profile and SPM measurements are excellent tools to assess the transport mechanism differentiating ballistic and diffusive graphene systems.

The dephasing relation (DR), a linearization of semiclassical fidelity, is generalized to include the overlap of “off-diagonal” elements. The accuracy of the formulation is tested in integrable and chaotic systems and its scaling with dimensionality is studied in a Caldeira-Leggett model with many degrees of freedom. It is shown that the DR is often in very good agreement with numerically analytic quantum results and frequently outperforms an alternative semiclassical treatment. Most importantly, since there is no computationally expensive prefactor, and Monte Carlo Metropolis sampling is used to facilitate the calculation, the DR is found to scale remarkably well with increasing dimension. We further demonstrate that a propagator based on the DR can include more quantum coherence and outperform other popular linearized semiclassical methods, such as forward-backward semiclassical dynamics (FBSD) and the linearized semiclassical initial value representation (LSC-IVR).

We apply quantum optimal control theory to establish a local voltage-control scheme that operates in conjunction with the numerically exact solution of the time-dependent Schrödinger equation. The scheme is demonstrated for high-fidelity coherent control of electronic charge in semiconductor double quantum dots. We find tailored gate voltages in the viable gigahertz regime that drive the system to a desired charge configuration with >99% yield. The results could be immediately verified in experiments and would play an important role in applications towards solid-state quantum computing.

We present a method for bridging the gap between the Dirac effective field theory and atomistic simulations in graphene based on the Husimi projection, allowing us to depict phenomena in graphene at arbitrary scales. This technique takes the atomistic wave function as an input, and produces semiclassical pictures of quasiparticles in the two Dirac valleys. We use the Husimi technique to produce maps of the scattering behavior of boundaries, giving insight into the properties of wave functions at energies both close to and far from the Dirac point. Boundary conditions play a significant role to the rise of Fano resonances, which we examine using the processed Husimi map to deepen our understanding of bond currents near resonance.

We present theory and calculations for coherent high-fidelity quantum control of many-particle states in semiconductor quantum wells. We show that coupling a two-electron double quantum dot to a terahertz optical source enables targeted excitations that are one to two orders of magnitude faster and significantly more accurate than those obtained with electric gates. The optical fields subject to physical constraints are obtained through quantum optimal control theory that we apply in conjunction with the numerically exact solution of the time-dependent Schrödinger equation. Our ability to coherently control arbitrary two-electron states, and to maximize the entanglement, opens up further perspectives in solid-state quantum information.

We introduce a method for classical trajectory calculations to simulate collisions between atoms and large rigid asymmetric-top molecules. We investigate the formation of molecule-helium complexes in buffer-gas cooling experiments at a temperature of 6.5 K for molecules as large as naphthalene. Our calculations show that the mean lifetime of the naphthalene-helium quasi-bound collision complex is not long enough for the formation of stable clusters under the experimental conditions. Our results suggest that it may be possible to improve the efficiency of the production of cold molecules in buffer-gas cooling experiments by increasing the density of helium. In addition, we find that the shape of molecules is important for the collisiondynamics when the vibrational motion of molecules is frozen. For some molecules, it is even more crucial than the number of accessible degrees of freedom. This indicates that by selecting molecules with suitable shape for buffer-gas cooling, it may be possible to cool molecules with a very large number of degrees of freedom.

Classical atom–diatom collisions at low velocities can be considered as a transient perturbation to the (integrable) diatomic system. We present an analysis that makes explicit the contributions of the terms of the Fourier expansion of the interaction potential to the changes in the molecular actions due to the collision process. Each term is associated with a resonance condition between the vibrational and rotational molecular frequencies, and leads to a vibrational, rotational or vibrotational contribution to the total action changes. The analysis is applied to the system Li2⁎Ne.

With the growth in interest in graphene, controlled nanoscale device geometries with complex form factors are now being studied and characterized. There is a growing need to understand new techniques to handle efficient electronic transport calculations for these systems. We present an algorithm that dramatically reduces the computational time required to find the local density of states and transmission matrix for open systems regardless of their topology or boundary conditions. We argue that the algorithm, which generalizes the recursive Green's function method by incorporating the reverse Cuthill-McKee algorithm for connected graphs, is ideal for calculating transmission through devices with multiple leads of unknown orientation and becomes a computational necessity when the input and output leads overlap in real space. This last scenario takes the Landauer-Buttiker formalism to general scattering theory in a computational framework that makes it tractable to perform full-spectrum calculations of the quantum scattering matrix in mesoscopic systems. We demonstrate the efficacy of these approaches on graphene stadiums, a system of recent scientific interest, and contribute to a physical understanding of Fano resonances which appear in these systems.

We study conductance fluctuations (CF) and the sensitivity of the conductance to the motion of a single scatterer in two-dimensional massless Dirac systems. Our extensive numerical study finds limits to the predicted universal value of CF. We find that CF are suppressed for ballistic systems near the Dirac point and approach the universal value at sufficiently strong disorder. The conductance of massless Dirac fermions is sensitive to the motion of a single scatterer. CF of order *e*^{2}/*h* result from the motion of a single impurity by a distance comparable to the Fermi wavelength. This result applies to graphene systems with a broad range of impurity strength and concentration while the dependence on the Fermi wavelength can be explored *via* gate voltages. Our prediction can be tested by comparing graphene samples with varying amounts of disorder and can be used to understand interference effects in mesoscopic graphene devices.

Motivated by recent experiments by the Westervelt group, which used a mobile tip to probe the electronic state of a segmented nanowire, we calculate shifts in Coulomb blockade peak positions, as a function of tip location, which we term “Coulomb blockade microscopy.” We show that if the tip can be brought sufficiently close to the nanowire, one can distinguish a high-density electronic liquid state from a Wigner-crystal state by microscopy with a weak-tip potential. In the opposite limit of a strongly negative tip potential, the potential depletes the electronic density under it and divides the quantum wire into two partitions. There the tip can push individual electrons from one partition to the other and the Coulomb blockade micrograph can clearly track such transitions. We show that this phenomenon can be used to qualitatively estimate the relative importance of the electron interaction compared to one-particle potential and kinetic energies. Finally, we propose that a weak-tip Coulomb blockade micrograph focusing on the transition between electron number N=0 and N=1 states may be used to experimentally map the one-particle potential landscape produced by impurities and inhomogeneities.

Graphene provides a fascinating testbed for new physics and exciting opportunities for future applications based on quantum phenomena. To understand the coherent flow of electrons through a graphene device, we employ a nanoscale probe that can access the relevant length scales—the tip of a liquid-He-cooled scanning probe microscope (SPM) capacitively couples to the graphene device below, creating a movable scatterer for electron waves. At sufficiently low temperatures and small size scales, the diffusive transport of electrons through graphene becomes coherent, leading to universal conductance fluctuations (UCF). By scanning the tip over a device, we map these conductance fluctuations versus scatterer position. We find that the conductance is highly sensitive to the tip position, producing δ*G* ~ *e*^{2}/*h* fluctuations when the tip is displaced by a distance comparable to half the Fermi wavelength. These measurements are in good agreement with detailed quantum simulations of the imaging experiment and demonstrate the value of a cooled SPM for probing coherent transport in graphene.

We study matter-wave scattering from an ultracold, many-body atomic system trapped in an optical lattice. The angular cross section of the target lattice for a matter wave is determined and is demonstrated to have a strong dependence on the many-body phase, superfluid, or Mott insulator. Analytical approaches are employed deep in the superfluid and Mott-insulator regimes, while intermediate points in the phase transition are treated numerically. Matter-wave scattering offers a convenient method for nondestructively probing the quantum many-body phase transition of atoms in an optical lattice.

Using a first-principles classical many-body simulation of a Hall bar, we study the necessary conditions for the formation of the Hall potential: (i) Ohmic contacts with metallic reservoirs, (ii) electron-electron interactions, and (iii) confinement to a finite system. By propagating thousands of interacting electrons over million time-steps we capture the build-up of the self-consistent potential. The microscopic model sheds light on the current injection process and directly links the Hall effect to specific boundary conditions at the particle reservoirs.

We study the quantum Hall effect (QHE) in graphene based on the current injection model, which takes into account the finite rectangular geometry with source and drain electrodes. In our model, the presence of disorder, the edge-state picture, extended states, and localized states, which are believed to be indispensable ingredients in describing the QHE, do not play an important role. Instead the boundary conditions during the injection into the graphene sheet, which are enforced by the presence of the Ohmic contacts, determine the current-voltage characteristics.