Combined R-matrix eigenstate basis set and finite-difference propagation method for the time-dependent Schrodinger equation: The one-electron case

In this work we present the theoretical framework for the solution of the time-dependent Schrödinger equation (TDSE) of atomic and molecular systems under strong electromagnetic fields with the configuration space of the electron’s coordinates separated over two regions; that is, regions I and II. In region I the solution of the TDSE is obtained by an R-matrix basis set representation of the time-dependent wave function. In region II a grid representation of the wave function is considered and propagation in space and time is obtained through the finite-difference method. With this, a combination of basis set and grid methods is put forward for tackling multiregion time-dependent problems. In both regions, a high-order explicit scheme is employed for the time propagation. While, in a purely hydrogenic system no approximation is involved due to this separation, in multielectron systems the validity and the usefulness of the present method relies on the basic assumption of R-matrix theory, namely, that beyond a certain distance (encompassing region I) a single ejected electron is distinguishable from the other electrons of the multielectron system and evolves there (region II) effectively as a one-electron system. The method is developed in detail for single active electron systems and applied to the exemplar case of the hydrogen atom in an intense laser field.

A discrete time-dependent method for metastable atoms and molecules in intense fields

The full-dimensional time-dependent Schrodinger equation for the electronic dynamics of single-electron systems in intense external fields is solved directly using a discrete method. Our approach combines the finite-difference and Lagrange mesh methods. The method is applied to calculate the quasienergies and ionization probabilities of atomic and molecular systems in intense static and dynamic electric fields. The gauge invariance and accuracy of the method is established. Applications to multiphoton ionization of positronium, the hydrogen atom and the hydrogen molecular ion are presented. At very high laser intensity, above the saturation threshold, we extend the method using a scaling technique to estimate the quasienergies of metastable states of the hydrogen molecular ion. The results are in good agreement with recent experiments. (C) 2004 American Institute of Physics.

Age-specific reproductive investment in female burying beetles: independent effects of state and risk of death

1.How much should an individual invest in reproduction as it grows older? Answering this question involves determining whether individuals measure their age as the time left for future reproduction or as the rate of deterioration in their state. Theory suggests that in the former case individuals should increase their allocation of resources to reproduction as opportunities for future breeding dwindle, and terminally invest when they breed for the last time. In the latter case they should reduce their investment in reproduction with age, either through adaptive reproductive restraint or as a passive by-product of senescence.
2.Here we present the results of experiments on female burying beetles, Nicrophorus vespilloides, in which we independently manipulated the perceived risk of death (by activating the immune system) and the extent of deterioration in state (by changing age of first reproduction and/or prior investment in reproduction).
3.We found that the risk of death and state each independently influenced the extent of reproductive investment. Specifically, we found a state-dependent decline in reproductive investment as females grew older that could be attributed to both adaptive reproductive restraint and senescence. A perceived increase in the risk of death, induced by activation of the immune system, caused females to switch from a strategy of reproductive restraint to terminal investment. Nevertheless, absolute reproductive investment was lower in older females, indicating constraints of senescence.
4.Our results show that a decline in reproductive investment with age does not necessarily constitute evidence of reproductive senescence but can also result from adaptive reproductive restraint.
5.Our results further suggest that the extent of reproductive investment is dependent on several different intrinsic cues and that the particular blend of cues available at any given age can yield very different patterns of investment. Perhaps this explains why age-related reproductive investment patterns seen in nature are so diverse.

A multiconfigurational time-dependent Hartree-Fock method for excited electronic states. I. General formalism and application to open-shell states

The solution of the time-dependent Schrodinger equation for systems of interacting electrons is generally a prohibitive task, for which approximate methods are necessary. Popular approaches, such as the time-dependent Hartree-Fock (TDHF) approximation and time-dependent density functional theory (TDDFT), are essentially single-configurational schemes. TDHF is by construction incapable of fully accounting for the excited character of the electronic states involved in many physical processes of interest; TDDFT, although exact in principle, is limited by the currently available exchange-correlation functionals. On the other hand, multiconfigurational methods, such as the multiconfigurational time-dependent Hartree-Fock (MCTDHF) approach, provide an accurate description of the excited states and can be systematically improved. However, the computational cost becomes prohibitive as the number of degrees of freedom increases, and thus, at present, the MCTDHF method is only practical for few-electron systems. In this work, we propose an alternative approach which effectively establishes a compromise between efficiency and accuracy, by retaining the smallest possible number of configurations that catches the essential features of the electronic wavefunction. Based on a time-dependent variational principle, we derive the MCTDHF working equation for a multiconfigurational expansion with fixed coefficients and specialise to the case of general open-shell states, which are relevant for many physical processes of interest. (C) 2011 American Institute of Physics. [doi: 10.1063/1.3600397]

Dissociative ionization of molecules in intense laser fields

Accurate and efficient grid based techniques for the solution of the time-dependent Schrodinger equation for few-electron diatomic molecules irradiated by intense, ultrashort laser pulses are described. These are based on hybrid finite-difference, Lagrange mesh techniques. The methods are applied in three scenarios, namely H-2(+) with fixed internuclear separation, H-2(+) with vibrating nuclei and H-2 with fixed internuclear separation and illustrative results presented.

Dynamic tunnelling ionization of H-+(2) in intense fields

Intense-field ionization of the hydrogen molecular ion by linearly polarized light is modelled by direct solution of the fixed-nuclei time-dependent Schrodinger equation and compared with recent experiments. Parallel transitions are calculated using algorithms which exploit massively parallel computers. We identify and calculate dynamic tunnelling ionization resonances that depend on laser wavelength and intensity, and molecular bond length. Results for lambda similar to 1064 nm are consistent with static tunnelling ionization. At shorter wavelengths lambda similar to 790 nm large dynamic corrections are observed. The results agree very well with recent experimental measurements of the ion spectra. Our results reproduce the single peak resonance and provide accurate ionization rate estimates at high intensities. At lower intensities our results confirm a double peak in the ionization rate as the bond length varies.

Ionization dynamics of laser-driven H-2(+)

We set out aspects of a numerical algorithm used in solving the full-dimensionality time-dependent Schrodinger equation describing the electronic motion of the hydrogen molecular ion driven by an intense, linearly polarized laser pulse aligned along the molecular axis. This algorithm has been implemented within the fixed inter-nuclear separation approximation in a parallel computer code, a brief summary of which is given. Ionization rates are calculated and compared with results from other methods, notably the time-independent Floquet method. Our results compare very favourably with the precise predictions of the Floquet method, although there is some disagreement with other wavepacket calculations. Visualizations of the electron dynamics are also presented in which electron rescattering is observed.

Geometrical dependence in photoionization of H-2(+) in high-intensity, high-frequency, ultrashort laser pulses

The ionization dynamics of H2 + exposed to high-intensity, high-frequency, ultrashort laser pulses is investigated with two theoretical approaches. The time-dependent Schrödinger equation is solved by a direct numerical method, and a simple two-center interference-diffraction model is studied. The energy and angular distributions of the photoelectron for various internuclear distances and relative orientations between the internuclear axis of the molecule and the polarization of the field are calculated. The main features of the photoelectron spectrum pattern are described well by the interference-diffraction model, and excellent quantitative agreement between the two methods is found. The effect of quantal vibration on the photoelectron spectrum is also calculated. We find that vibrational average produces some broadening of the main features, but that the patterns remain clearly distinguishable.

Dissociation of H2+ from a short, intense, infrared laser pulse: proton emission spectra and pulse calibration

The proton energy spectrum from photodissociation of the hydrogen molecular ion by short intense pulses of infrared light is calculated. The time-dependent Schrödinger equation is discretized and integrated. For few-cycle pulses one can resolve vibrational structure, arising from the experimental preparation of the molecular ion. We calculate the corresponding energy spectrum and analyse the dependence on the pulse time delay, pulse length and intensity of the laser for ? ~ 790 nm. We conclude that the proton spectrum is a sensitive probe of both the vibrational populations and phases, and allows us to distinguish between adiabatic and nonadiabatic dissociation. Furthermore, the sensitivity of the proton spectrum from H2+ is a practical means of calibrating the pulse. Our results are compared with recent measurements of the proton spectrum for 65 fs pulses using a Ti:Sapphire laser (? ~ 790 nm) including molecular orientation and focal-volume averaging. Integrating over the laser focal volume, for the intensity I ~ 3 × 1015 W cm-2, we find our results are in excellent agreement with these experiments.

Transient receptor potential melastatin 8 (TRPM8) channel involvement in the regulation of vascular tone

The transient receptor potential melastatin 8 (TRPM8) channel has been characterized as a cold and menthol receptor expressed in a subpopulation of sensory neurons but was recently identified in other tissues, including the respiratory tract, urinary system, and vasculature. Thus TRPM8 may play multiple functional roles, likely to be in a tissue- and activation state-dependent manner. We examined the TRPM8 channel presence in large arteries from rats and the functional consequences of their activation. We also aimed to examine whether these channels contribute to control of conscious human skin blood flow. TRPM8 mRNA and protein were detected in rat tail, femoral and mesenteric arteries, and thoracic aorta. This was confirmed in single isolated vascular myocytes by immunocytochemistry. Isometric contraction studies on endothelium-denuded relaxed rat vessels found small contractions on application of the TRPM8-specific agonist menthol (300 microM). However, both menthol and another agonist icilin (50 microM) caused relaxation of vessels precontracted with KCl (60 mM) or the alpha-adrenoceptor agonist phenylephrine (2 microM) and a reduction in sympathetic nerve-mediated contraction. These effects were antagonized by bromoenol lactone treatment, suggesting the involvement of Ca(2+)-independent phospholipase A(2) activation in TRPM8-mediated vasodilatation. In thoracic aorta with intact endothelium, menthol-induced inhibition of KCl-induced contraction was enhanced. This was unaltered by preincubation with either N(omega)-nitro-l-arginine methyl ester (l-NAME; 100 nM), a nitric oxide synthase inhibitor, or the ACh receptor antagonist atropine (1 microM). Application of menthol (3% solution, topical application) to skin caused increased blood flow in conscious humans, as measured by laser Doppler fluximetry. Vasodilatation was markedly reduced or abolished by prior application of l-NAME (passive application, 10 mM) or atropine (iontophoretic application, 100 nM, 30 s at 70 microA). We conclude that TRPM8 channels are present in rat artery vascular smooth muscle and on activation cause vasoconstriction or vasodilatation, dependent on previous vasomotor tone. TRPM8 channels may also contribute to human cutaneous vasculature control, likely with the involvement of additional neuronal mechanisms.

Ring currents in azulene

We propose a self consistent polarisable ion tight binding theory for the study of push-pull processes in aromatic molecules. We find that the method quantitatively reproduces ab initio calculations of dipole moments and polarisability. We apply the scheme in a simulation which solves the time dependent Schroedinger equation to follow the relaxation of azulene from the second excited to the ground states. We observe rather spectacular oscillating ring currents which we explain in terms of interference between the HOMO and LUMO states.

Doppler cooling of gallium atoms: 2. Simulation in complex multilevel systems

This paper derives a general procedure for the numerical solution of the Lindblad equations that govern the coherences arising from multicoloured light interacting with a multilevel system. A systematic approach to finding the conservative and dissipative terms is derived and applied to the laser cooling of p-block elements. An improved numerical method is developed to solve the time-dependent master equation and results are presented for transient cooling processes. The method is significantly more robust, efficient and accurate than the standard method and can be applied to a broad range of atomic and molecular systems. Radiation pressure forces and the formation of dynamic dark states are studied in the gallium isotope 66Ga.

Double-Electron Above-Threshold Ionization Resonances as Interference Phenomena

We report calculations of double ionization energy spectra and momentum distributions of laser-driven helium due to few-cycle pulses of wavelength 195 nm. The results are obtained from full-dimensional numerical integration of the two electron time-dependent Schr¨odinger equation. A momentum-space analysis of doubly ionizing wavepackets shows that the concentric-ring structure of above-threshold double ionization, together with the associated structure of peaks in the total kinetic energy spectrum, may be attributed to wavepacket interference effects, where at least two doubly-ionizing wavepackets from different recollision events populate the same spatial hemisphere.

From the UV to the static-field limit: rates and scaling laws of intense-field ionization of helium

We present high-accuracy calculations of ionization rates of helium at UV (195 nm) wavelengths. The data are obtained from full-dimensionality integrations of the helium-laser time-dependent Schrödinger equation. Comparison is made with our previously obtained data at 390 nm and 780 nm. We show that scaling laws introduced by Parker et al extend unmodified from the near-infrared limit into the UV limit. Static-field ionization rates of helium are also obtained, again from time-dependent full-dimensionality integrations of the helium Schrödinger equation. We compare the static-field ionization results with those of Scrinzi et al and Themelis et al, who also treat the full-dimensional helium atom, but with time-independent methods. Good agreement is obtained.

The RMT method for many-electron atomic systems in intense short-pulse laser light

We describe a new ab initio method for solving the time-dependent Schrödinger equation for multi-electron atomic systems exposed to intense short-pulse laser light. We call the method the R-matrix with time-dependence (RMT) method. Our starting point is a finite-difference numerical integrator (HELIUM), which has proved successful at describing few-electron atoms and atomic ions in strong laser fields with high accuracy. By exploiting the R-matrix division-of-space concept, we bring together a numerical method most appropriate to the multi-electron finite inner region (R-matrix basis set) and a different numerical method most appropriate to the one-electron outer region (finite difference). In order to exploit massively parallel supercomputers efficiently, we time-propagate the wavefunction in both regions by employing Arnoldi methods, originally developed for HELIUM.