Schroedinger Equation Research Articles

Resonant harmonic generation and collective spin rotations in electrically driven quantum dots

Spin rotations induced by an AC electric field in a two-electron double quantum dot are studied by an exact numerical solution of the time dependent Schroedinger equation in the context of recent electric dipole spin resonance experiments based on the Pauli blockade. We demonstrate that the splitting of the main resonance line by the spin exchange coupling is accompanied by the appearance of fractional resonances and that both these effects are triggered by interdot tunnel coupling. We find that the AC driven system generates residual but distinct harmonics of the driving frequency which are amplified when tuned to the main transition frequency. The mechanism is universal for electron systems in electrically driven potentials and works also in the absence of electron-electron interaction or spin-orbit coupling.

Propagating Wave Patterns in a Derivative Nonlinear Schrödinger System with Quintic Nonlinearity

Exact expressions are obtained for a diversity of propagating patterns for a derivative nonlinear Schrodinger equation with the quintic nonlinearity. These patterns include bright pulses, fronts and dark solitons. The evolution of the wave envelope is determined via a pair of integrals of motion, and reduction is achieved to Jacobi elliptic cn and dn function representations. Numerical simulations are performed to establish the existence of parameter ranges for stability. The derivative quintic nonlinear Schrodinger model equations investigated here are relevant in the analysis of strong optical signals propagating in spatial or temporal waveguides.

The concept of quasi-integrability for modified non-linear Schrödinger models

We consider modifications of the nonlinear Schrodinger model (NLS) to look at the recently introduced concept of quasi-integrability. We show that such models possess an infinite number of quasi-conserved charges which present intriguing properties in relation to very specific space-time parity transformations. For the case of two-soliton solutions where the fields are eigenstates of this parity, those charges are asymptotically conserved in the scattering process of the solitons. Even though the charges vary in time their values in the far past and the far future are the same. Such results are obtained through analytical and numerical methods, and employ adaptations of algebraic techniques used in integrable field theories. Our findings may have important consequences on the applications of these models in several areas of non-linear science. We make a detailed numerical study of the modified NLS potential of the form V = |psi|^(2(2+epsilon)), with epsilon being a perturbation parameter. We perform numerical simulations of the scattering of solitons for this model and find a good agreement with the results predicted by the analytical considerations. Our paper shows that the quasi-integrability concepts recently proposed in the context of modifications of the sine-Gordon model remain valid for perturbations of the NLS model.

Extensive Adiabatic Invariants for Nonlinear Chains

We look for extensive adiabatic invariants in nonlinear chains in the thermodynamic limit. Considering the quadratic part of the Klein-Gordon Hamiltonian, by a linear change of variables we transform it into a sum of two parts in involution. At variance with the usual method of introducing normal modes, our constructive procedure allows us to exploit the complete resonance, while keeping the extensive nature of the system. Next we construct a nonlinear approximation of an extensive adiabatic invariant for a perturbation of the discrete nonlinear Schrodinger model. The fluctuations of this quantity are controlled via Gibbs measure estimates independent of the system size, for a large set of initial data at low specific energy. Finally, by numerical calculations we show that our adiabatic invariant is well conserved for times much longer than predicted by our first order theory, with fluctuation much smaller than expected according to standard statistical estimates.

Fano resonance due to discrete breather in nonlinear Klein–Gordon lattice in metamaterials

The richer variety of Klein–Gordon basis is already established for discrete breathers in metamatetrials. Based on this attempt, we show various anomalous Fano resonance behaviors that have been experimentally observed, but cannot be explained by nonlinear Schrodinger model. Certain material parameters of Klein–Gordon lattice in metamaterials are related for the first time with characteristics of Fano resonance, which can be utilized for beam filtering and for high-resolution biological sensing technology. Although relations with coupling and other parameters exist, the most remarkable relation is observed with linear permittivity that could control the wave transmission characteristics in metamaterials for applications in optical engineering.

Nonequilibrium discrete nonlinear Schrödinger equation

We study nonequilibrium steady states of the one-dimensional discrete nonlinear Schrödinger equation. This system can be regarded as a minimal model for the stationary transport of bosonic particles such as photons in layered media or cold atoms in deep optical traps. Due to the presence of two conserved quantities, namely, energy and norm (or number of particles), the model displays coupled transport in the sense of linear irreversible thermodynamics. Monte Carlo thermostats are implemented to impose a given temperature and chemical potential at the chain ends. As a result, we find that the Onsager coefficients are finite in the thermodynamic limit, i.e., transport is normal. Depending on the position in the parameter space, the "Seebeck coefficient" may be either positive or negative. For large differences between the thermostat parameters, density and temperature profiles may display an unusual nonmonotonic shape. This is due to the strong dependence of the Onsager coefficients on the state variables.

Particle motion in a deformed potential using a transformed oscillator basis

The quantum description of a particle moving in a deformed potential is investigated. A pseudostate (PS) basis is used to represent the states of the composite system. This PS basis is obtained by diagonalizing the system Hamiltonian in a family of square integrable functions. In this work the Transformed Harmonic Oscillator (THO) functions, obtained from the solutions of the Harmonic Oscillator using a Local Scale Transformation (LST), are used. The proposed method is applied to the 11Be nucleus, treated in a two-body model (10Be+n). Both structure and reaction observables have been studied. Wavefunctions and energies obtained for the bound states and some low-lying resonances are compared with those obtained by direct integration of the Schroedinger equation. The dipole and quadrupole electric transition probabilities for the low-energy continuum have been calculated in the THO basis, and compared with the exact distributions obtained with the scattering states. Finally, the method is applied to describe the 11Be states in the Coulomb breakup of 11Be+208Pb at 69 MeV/nucleon. The energy and angular distributions of the exclusive breakup have been calculated using the Equivalent Photon Method, including both E1 and E2 contributions. The calculated distributions are found to be in good agreement with the available experimental data from RIKEN [Phys. Rev. C70, 054606]. At the very forward angles, the cross section is completely dominated by the dipole couplings.

The Quadrupole Response of Borromean Bosonic Trimers

The inelastic response of a Borromean trimer composed of three identical bosons is explored. To this end we use the quadrupole approximation as the excitation mechanism of the bosonic system from its ground state. We utilize the hyperspherical-harmonics expansion to solve the Schroedinger equation and the Lorentz integral transform method to calculate the reaction. It is found that the magnitude of the response function and corresponding sum rules increase as a power law when approaching the three-body threshold. It is also found that this increase is governed by unnatural exponents.

Computational study of the rovibrational spectrum of (OCS)2

In this paper, we report a new intermolecular potential energy surface and rovibrational transition frequencies and line strengths computed for the OCS dimer. The potential is made by fitting energies obtained from explicitly correlated coupled-cluster calculations and fit using an interpolating moving least squares method. The rovibrational Schroedinger equation is solved with a symmetry-adapted Lanczos algorithm and an uncoupled product basis set. All four intermolecular coordinates are included in the calculation. On the potential energy surface we find, previously unknown, cross-shaped isomers and also polar and non-polar isomers. The associated wavefunctions and energy levels are presented. To identify polar and cross states we use both calculations of line strengths and vibrational parent analysis. Calculated rotational constants differ from their experimental counterparts by less than 0.001 cm(-1).

Broadband propagation and photoacoustic time reversal imaging using k-space methods

Numerical, time domain, models of broadband acoustic propagation using k-space methods will be described and applied to the problem of image reconstruction in photoacoustic imaging. k-space methods are a subset of time-stepping pseudospectral methods, which use FFTs to calculate field gradients, with an additional correction to make the solutions exact in the homogeneous case. Furthermore, they are closely related to a number of methods in the mathematical, engineering and physical sciences literature, including non-standard finite difference schemes, exponential integrators, wavenumber domain reverse time migration, and wave propagator models developed to solve the time-dependent Schroedinger equation. These connections will be discussed and used to illuminate the advantages of the k-space approach for large scale modelling of broadband acoustic waves. To illustrate the method, a fluid-based k-space model will be applied to the inverse acoustic initial value problem encountered in tomographic photoacoustic imaging. It is well known that this can be solved, in the linear case, using a time-reversed numerical model. Extensions from this basic case to absorbing and non-linear media will be explored.

Semiclassical initial value representation with complex trajectories

Semiclassical time-dependent solutions to the Schroedinger equation can be obtained by use of integral formulas known as initial value representations (IVRs). Over the last few decades IVR formulas have established a sucessful reputation, most notably in the field of physical chemistry, for it is often found that the quantum behavior of molecular systems is very well approximated by semiclassical wavefunctions. The IVR input consists in a swarm of independent classical trajectories, each of which obeys the dynamics dictated by the system’s unquantized Hamiltonian. In a recent work [Aguiar et al, Chem. Phys. 370 (2010)] an IVR expression based on complex classical trajectories was formulated. The complex nature of the classical trajectories in this particular formula is a consequence of employing minimum uncertainty wavepackets (coherent states) to represent the time evolution operator upon which the formula is constructed. Here, an improved and simpler version of this complex IVR formula is presented. In order to test its accuracy, we apply it to a simple model system.

Solitons supported by singular spatial modulation of the Kerr nonlinearity

We introduce a setting based on the one-dimensional (1D) nonlinear Schroedinger equation (NLSE) with the self-focusing (SF) cubic term modulated by a singular function of the coordinate, |x|^{-a}. It may be additionally combined with the uniform self-defocusing (SDF) nonlinear background, and with a similar singular repulsive linear potential. The setting, which can be implemented in optics and BEC, aims to extend the general analysis of the existence and stability of solitons in NLSEs. Results for fundamental solitons are obtained analytically and verified numerically. The solitons feature a quasi-cuspon shape, with the second derivative diverging at the center, and are stable in the entire existence range, which is 0 < a < 1. Dipole (odd) solitons are found too. They are unstable in the infinite domain, but stable in the semi-infinite one. In the presence of the SDF background, there are two subfamilies of fundamental solitons, one stable and one unstable, which exist together above a threshold value of the norm (total power of the soliton). The system which additionally includes the singular repulsive linear potential emulates solitons in a uniform space of the fractional dimension, 0 < D < 1. A two-dimensional extension of the system, based on the quadratic nonlinearity, is formulated too.

Vortex families near a spectral edge in the Gross-Pitaevskii equation with a two-dimensional periodic potential

We examine numerically vortex families near band edges of the Bloch wave spectrum for the Gross-Pitaevskii equation with two-dimensional periodic potentials and for the discrete nonlinear Schrödinger equation. We show that besides vortex families that terminate at a small distance from the band edges via fold bifurcations, there exist vortex families that are continued all the way to the band edges.

InAs/InAsSb strain balanced superlattices for optical detectors: Material properties and energy band simulations

InAsSb/InAs type II strain balanced superlattices lattice matched to GaSb have recently been proposed as an alternative to InAs/(In)GaSb short period superlattices for mid- to long infrared photodetectors. Photoluminescence data at 4 K of OMVPE grown InAsSb (multi-) quantum wells in an InAs matrix on InAs and GaSb substrates is presented for Sb compositions between 4% and 27%. The measured transition energies are simulated with a self-consistent Poisson and Schroedinger equation solver that includes strain and band-offsets. The fitted parameters are then used to predict the type II transition energies of InAsSb/InAs strain balanced superlattice absorber stacks at 77 K for different compositions and periods. The optical matrix element was calculated and compared with InAs/(In)GaSb superlattices. The InAsSb/InAs structures can be designed with higher or equal matrix elements for longer periods. Finally, the initial optical response data of an unoptimized strain balanced InAs0.79Sb0.21/InAs detector with a 40 nm period are shown. Its cutoff wavelength is 0.15 eV (8.5 μm), in good agreement with the predicted type II transition energy of 0.17 eV.

Multiscale perturbative approach toSU(2)-Higgs classical dynamics: Stability of nonlinear plane waves and bounds of the Higgs field mass

We study the classical dynamics of SU(2)-Higgs field theory using multiple scale perturbation theory. In the spontaneously broken phase, assuming small perturbations of the Higgs field around its vacuum expectation value, we derive a nonlinear Schroedinger equation and study the stability of its nonlinear plane wave solutions. The latter, turn out to be stable only if the Higgs amplitude is an order of magnitude smaller than that of the gauge field. In this case, the Higgs field mass possesses some bounds which may be relevant to the search for the Higgs particle at ongoing experiments.

Anderson localization in nonlocal nonlinear media

We theoretically and numerically investigate the effect of focusing and defocusing nonlinearities on Anderson localization in highly nonlocal media. A perturbative approach is developed to solve the nonlocal nonlinear Schrödinger equation in the presence of a random potential, showing that nonlocality stabilizes Anderson states.

Variational determination of approximate bright matter-wave soliton solutions in anisotropic traps

We consider the ground state of an attractively-interacting atomic Bose-Einstein condensate in a prolate, cylindrically symmetric harmonic trap. If a true quasi-one-dimensional limit is realized, then for sufficiently weak axial trapping this ground state takes the form of a bright soliton solution of the nonlinear Schroedinger equation. Using analytic variational and highly accurate numerical solutions of the Gross-Pitaevskii equation we systematically and quantitatively assess how soliton-like this ground state is, over a wide range of trap and interaction strengths. Our analysis reveals that the regime in which the ground state is highly soliton-like is significantly restricted, and occurs only for experimentally challenging trap anisotropies. This result, and our broader identification of regimes in which the ground state is well-approximated by our simple analytic variational solution, are relevant to a range of potential experiments involving attractively-interacting Bose-Einstein condensates.

Electron-impact ionization of Li2using a time-dependent close-coupling method

Electron-impact ionization cross sections for H{sub 2} are calculated using a nonperturbative time-dependent close-coupling method. In a standard frozen-core approximation, the six-dimensional wave function for the valence target electron and the incident projectile electron is expanded in products of rotational functions. The time-dependent Schroedinger equation for the two-electron system is then reduced to a set of close-coupled partial differential equations for the four-dimensional expansion functions in (r{sub 1},{theta}{sub 1},r{sub 2},{theta}{sub 2}) center-of-mass spherical polar coordinates. The nonperturbative close-coupling results are found to be over a factor of 2 lower than perturbative distorted-wave results, but in excellent agreement with experimental measurements, at incident electron energies near the peak of the total integrated cross section.

Controllability of the cubic Schroedinger equation via a low-dimensional source term

We study controllability of $d$-dimensional defocusing cubic Schroe-din-ger equation under periodic boundary conditions. The control is applied additively, via a source term, which is a linear combination of few complex exponentials (modes) with time-variant coefficients - controls. We manage to prove that controlling $2^d$ modes one can achieve controllability of the equation in each finite-dimensional projection of the evolution space $H^{s}(\mathbb{T}^d), s>d/2$, as well as approximate controllability in $H^{s}(\mathbb{T}^d)$. We also present a negative result regarding exact controllability of cubic Schroedinger equation via a finite-dimensional source term.

Transverse instability for periodic waves of KP-I and Schroedinger equations

Transverse instability for periodic waves of KP-I and Schroedinger equations