Nebular gas drag and co-orbital system dynamics

Aims. We study trajectories of planetesimals whose orbits decay due to gas drag in a primordial solar nebula and are perturbed by the gravity of the secondary body on an eccentric orbit whose mass ratio takes values from mu(2) = 10(-7) to mu(2) = 10(-3) increasing ten times at each step. Each planetesimal ultimately suffers one of the three possible fates: (1) trapping in a mean motion resonance with the secondary body; (2) collision with the secondary body and consequent increase of its mass; or (3) diffusion after crossing the orbit of the secondary body.Methods. We take the Burlirsh-Stoer numerical algorithm in order to integrate the Newtonian equations of the planar, elliptical restricted three-body problem with the secondary body and the planetesimal orbiting the primary. It is assumed that there is no interaction among planetesimals, and also that the gas does not affect the orbit of the secondary body.Results. The results show that the optimal value of the gas drag constant k for the 1: 1 resonance is between 0.9 and 1.25, representing a meter size planetesimal for each AU of orbital radius. In this study, the conditions of the gas drag are such that in theory, L4 no longer exists in the circular case for a critical value of k that defines a limit size of the planetesimal, but for a secondary body with an eccentricity larger than 0.05 when mu(2) = 10(-6), it reappears. The decrease of the cutoff collision radius increase the difusions but does not affect the distribution of trapping. The contribution to the mass accretion of the secondary body is over 40% with a collision radius 0.05R(Hill) and less than 15% with 0.005R(Hill) for mu(2) = 10(-7). The trappings no longer occur when the drag constant k reachs 30. That means that the size limit of planetesimal trapping is 0.2 m per AU of orbital radius. In most cases, this accretion occurs for a weak gas drag and small secondary eccentricity. The diffusions represent most of the simulations showing that gas drag is an efficient process in scattering planetesimals and that the trapping of planetesimals in the 1: 1 resonance is a less probable fate. These results depend on the specific drag force chosen.

Classical Langevin dynamics of a charged particle moving on a sphere and diamagnetism: A surprise

It is generally known that the orbital diamagnetism of a classical system of charged particles in thermal equilibrium is identically zero —the Bohr-van Leeuwen theorem. Physically, this null result derives from the exact cancellation of the orbital diamagnetic moment associated with the complete cyclotron orbits of the charged particles by the paramagnetic moment subtended by the incomplete orbits skipping the boundary in the opposite sense. Motivated by this crucial but subtle role of the boundary, we have simulated here the case of a finite but unbounded system, namely that of a charged particle moving on the surface of a sphere in the presence of an externally applied uniform magnetic field. Following a real space-time approach based on the classical Langevin equation, we have computed the orbital magnetic moment that now indeed turns out to be non-zero and has the diamagnetic sign. To the best of our knowledge, this is the first report of the possibility of finite classical diamagnetism in principle, and it is due to the avoided cancellation.

Embracing Integrability in Stellar and Planetary Dynamics

Hamiltonian systems in stellar and planetary dynamics are typically near integrable. For example, Solar System planets are almost in two-body orbits, and in simulations of the Galaxy, the orbits of stars seem regular. For such systems, sophisticated numerical methods can be developed through integrable approximations. Following this theme, we discuss three distinct problems. We start by considering numerical integration techniques for planetary systems. Perturbation methods (that utilize the integrability of the two-body motion) are preferred over conventional "blind" integration schemes. We introduce perturbation methods formulated with Cartesian variables. In our numerical comparisons, these are superior to their conventional counterparts, but, by definition, lack the energy-preserving properties of symplectic integrators. However, they are exceptionally well suited for relatively short-term integrations in which moderately high positional accuracy is required. The next exercise falls into the category of stability questions in solar systems. Traditionally, the interest has been on the orbital stability of planets, which have been quantified, e.g., by Liapunov exponents. We offer a complementary aspect by considering the protective effect that massive gas giants, like Jupiter, can offer to Earth-like planets inside the habitable zone of a planetary system. Our method produces a single quantity, called the escape rate, which characterizes the system of giant planets. We obtain some interesting results by computing escape rates for the Solar System. Galaxy modelling is our third and final topic. Because of the sheer number of stars (about 10^11 in Milky Way) galaxies are often modelled as smooth potentials hosting distributions of stars. Unfortunately, only a handful of suitable potentials are integrable (harmonic oscillator, isochrone and Stäckel potential). This severely limits the possibilities of finding an integrable approximation for an observed galaxy. A solution to this problem is torus construction; a method for numerically creating a foliation of invariant phase-space tori corresponding to a given target Hamiltonian. Canonically, the invariant tori are constructed by deforming the tori of some existing integrable toy Hamiltonian. Our contribution is to demonstrate how this can be accomplished by using a Stäckel toy Hamiltonian in ellipsoidal coordinates.

Dynamics of prolate spheroidal mass distributions with varying eccentricity

In this paper we calculate the potential for a prolate spheroidal distribution as in a dark matter halo with a radially varying eccentricity. This is obtained by summing up the shell-by-shell contributions of isodensity surfaces, which are taken to be concentric and with a common polar axis and with an axis ratio that varies with radius. Interestingly, the constancy of potential inside a shell is shown to be a good approximation even when the isodensity contours are dissimilar spheroids, as long as the radial variation in eccentricity is small as seen in realistic systems. We consider three cases where the isodensity contours are more prolate at large radii, or are less prolate or have a constant eccentricity. Other relevant physical quantities like the rotation velocity, the net orbital and vertical frequency due to the halo and an exponential disc of finite thickness embedded in it are obtained. We apply this to the kinematical origin of Galactic warp, and show that a prolate-shaped halo is not conducive to making long-lived warps - contrary to what has been proposed in the literature. The results for a prolate mass distribution with a variable axis ratio obtained are general, and can be applied to other astrophysical systems, such as prolate bars, for a more realistic treatment.

Quantum chemical study on influence of intermolecular hydrogen bonding on the geometry, the atomic charges and the vibrational dynamics of 2,6-dichlorobenzonitrile

FT-IR (4000-400 cm(-1)) and FT-Raman (4000-200 cm(-1)) spectral measurements on solid 2,6-dichlorobenzonitrile (2,6-DCBN) have been done. The molecular geometry, harmonic vibrational frequencies and bonding features in the ground state have been calculated by density functional theory at the B3LYP/6-311++G (d,p) level. A comparison between the calculated and the experimental results covering the molecular structure has been made. The assignments of the fundamental vibrational modes have been done on the basis of the potential energy distribution (PED). To investigate the influence of intermolecular hydrogen bonding on the geometry, the charge distribution and the vibrational spectrum of 2,6-DCBN; calculations have been done for the monomer as well as the tetramer. The intermolecular interaction energies corrected for basis set superposition error (BSSE) have been calculated using counterpoise method. Based on these results, the correlations between the vibrational modes and the structure of the tetramer have been discussed. Molecular electrostatic potential (MEP) contour map has been plotted in order to predict how different geometries could interact. The Natural Bond Orbital (NBO) analysis has been done for the chemical interpretation of hyperconjugative interactions and electron density transfer between occupied (bonding or lone pair) orbitals to unoccupied (antibonding or Rydberg) orbitals. UV spectrum was measured in methanol solution. The energies and oscillator strengths were calculated by Time Dependent Density Functional Theory (TD-DFT) and matched to the experimental findings. TD-DFT method has also been used for theoretically studying the hydrogen bonding dynamics by monitoring the spectral shifts of some characteristic vibrational modes involved in the formation of hydrogen bonds in the ground and the first excited state. The C-13 nuclear magnetic resonance (NMR) chemical shifts of the molecule were calculated by the Gauge independent atomic orbital (GIAO) method and compared with experimental results. Standard thermodynamic functions have been obtained and changes in thermodynamic properties on going from monomer to tetramer have been presented. (C) 2013 Elsevier B.V. All rights reserved.

QM/MM and Classical Molecular Dynamics Simulation of Histidine-Tagged Peptide Immobilization on Nickel Surface

The hybrid quantum mechanics (QM) and molecular mechanics (MM) method is employed to simulate the His-tagged peptide adsorption to ionized region of nickel surface. Based on the previous experiments, the peptide interaction with one Ni ion is considered. In the QM/MM calculation, the imidazoles on the side chain of the peptide and the metal ion with several neighboring water molecules are treated as QM part calculated by “GAMESS”, and the rest atoms are treated as MM part calculated by “TINKER”. The integrated molecular orbital/molecular mechanics (IMOMM) method is used to deal with theQMpart with the transitional metal. By using the QM/MM method, we optimize the structure of the synthetic peptide chelating with a Ni ion. Different chelate structures are considered. The geometry parameters of the QM subsystem we obtained by QM/MM calculation are consistent with the available experimental results. We also perform a classical molecular dynamics (MD) simulation with the experimental parameters for the synthetic peptide adsorption on a neutral Ni(1 0 0) surface. We find that half of the His-tags are almost parallel with the substrate, which enhance the binding strength. Peeling of the peptide from the Ni substrate is simulated in the aqueous solvent and in vacuum, respectively. The critical peeling forces in the two environments are obtained. The results show that the imidazole rings are attached to the substrate more tightly than other bases in this peptide.

QM/MM and Classical Molecular Dynamics Simulation of Histidine-Tagged Peptide Immobilization on Nickel Surface

The hybrid quantum mechanics (QM) and molecular mechanics (MM) method is employed to simulate the His-tagged peptide adsorption to ionized region of nickel surface. Based on the previous experiments, the peptide interaction with one Ni ion is considered. In the QM/MM calculation, the imidazoles on the side chain of the peptide and the metal ion with several neighboring water molecules are treated as QM part calculated by "GAMESS", and the rest atoms are treated as MM part calculated by "TINKER". The integrated molecular orbital/molecular mechanics (IMOMM) method is used to deal with the QM part with the transitional metal. By using the QM/MM method, we optimize the structure of the synthetic peptide chelating with a Ni ion. Different chelate structures are considered. The geometry parameters of the QM subsystem we obtained by QM/MM calculation are consistent with the available experimental results. We also perform a classical molecular dynamics (MD) simulation with the experimental parameters for the synthetic peptide adsorption on a neutral Ni(100) surface. We find that half of the His-tags are almost parallel with the substrate, which enhance the binding strength. Peeling of the peptide from the Ni substrate is simulated in the aqueous solvent and in vacuum, respectively. The critical peeling forces in the two environments are obtained. The results show that the in-tidazole rings are attached to the substrate more tightly than other bases in this peptide.

Dynamics of rotationally resolved multiphoton ionization processes in molecules

This dissertation presents the results of studies of several rotationally- resolved resonance enhanced multiphoton ionization (REMPI) processes in some simple molecular systems. The objective of these studies is to quantitatively identify the underlying dynamics of this highly state-specific process which utilizes the narrow bandwidth radiation of a laser to ionize a molecule by first preparing an excited state via multiphoton absorption and subsequently ionizing that state before it can decay. Coupled with high-resolution photoelectron spectroscopy, REMPI is clearly an important probe of molecular excited states and their photoioniza tion dynamics.

A key feature of our studies is that they are carried out using accurate Hartree-Fock orbitals to describe the photoelectron orbitals of the molecular ions. The use of such photoelectron orbitals is important in rotationally-resolved studies where the angular momentum coupling in the photoelectron orbital plays a significant role in the photoionization dynamics. In these studies the Hartree-Fock molecular molecular photoelectron orbitals are obtained by numerical solution of a Lippmann-Schwinger integral equation.

Studies reported here include investigations of (i) ionic rotational branching ratios and their energy dependence for REMPI via the A^2Σ^+(3sσ) and D^2Σ^+(3pσ)states of NO, (ii) the influence of angular momentum constraints on branching ratios at low photoelectron energies for REMPI via low-J levels of the resonant intermediate state, (iii) the strong dependence of photoelectron angular distributions on final ionic rotational state and on the alignment in REMPI of the A^2Σ^+ state of NO, (iv) vibrational state dependence of ionic rotational branching ratios arising from rapid orbital evolution in resonant states (E'^2Σ^+(3pσ) of CH), (v) the influence of rovibronic interactions on the rotational branching ratios seen in REMPI via the D^2Σ^+(3pσ) state of NO, and (vi) effects of laser intensity on the photoionization dynamics of REMPI.

Multi-objective optimization for multiphase orbital rendezvous missions

A multi-objective optimization approach was proposed for multiphase orbital rendezvous missions and validated by application to a representative numerical problem. By comparing the Pareto fronts obtained using the proposed method, the relationships between the three objectives considered are revealed, and the influences of other mission parameters, such as the sensors' field of view, can also be analyzed effectively. For multiphase orbital rendezvous missions, the tradeoff relationships between the total velocity increment and the trajectory robustness index as well as between the total velocity increment and the time of flight are obvious and clear. However, the tradeoff relationship between the time of flight and the trajectory robustness index is weak, especially for the four- and five-phase missions examined. The proposed approach could be used to reorganize a stable rendezvous profile for an engineering rendezvous mission, when there is a failure that prevents the completion of the nominal mission.

Half-precessional dynamics of monsoon rainfall near the East African Equator

External climate forcings-such as long-term changes in solar insolation-generate different climate responses in tropical and high latitude regions(1). Documenting the spatial and temporal variability of past climates is therefore critical for understanding how such forcings are translated into regional climate variability. In contrast to the data-richmiddle and high latitudes, high-quality climate-proxy records from equatorial regions are relatively few(2-4), especially from regions experiencing the bimodal seasonal rainfall distribution associated with twice-annual passage of the Intertropical Convergence Zone. Here we present a continuous and well-resolved climate-proxy record of hydrological variability during the past 25,000 years from equatorial East Africa. Our results, based on complementary evidence from seismic-reflection stratigraphy and organic biomarker molecules in the sediment record of Lake Challa near Mount Kilimanjaro, reveal that monsoon rainfall in this region varied at half-precessional (similar to 11,500-year) intervals in phase with orbitally controlled insolation forcing. The southeasterly and northeasterly monsoons that advect moisture from the western Indian Ocean were strengthened in alternation when the inter-hemispheric insolation gradient was at a maximum; dry conditions prevailed when neither monsoon was intensified and modest local March or September insolation weakened the rain season that followed. On sub-millennial timescales, the temporal pattern of hydrological change on the East African Equator bears clear high-northern-latitude signatures, but on the orbital timescale it mainly responded to low-latitude insolation forcing. Predominance of low-latitude climate processes in this monsoon region can be attributed to the low-latitude position of its continental regions of surface air flow convergence, and its relative isolation from the Atlantic Ocean, where prominent meridional overturning circulation more tightly couples low-latitude climate regimes to high-latitude boundary conditions.

Palaeoenvironmental changes since the Last Glacial Maximum: Patterns, timing and dynamics throughout South America

The vast diversity of present vegetation and environments that occur throughout South America (12°N to 56°S) is the result of diverse processes that have been operating and interacting at different spatial and temporal scales. Global factors, such as the concentration of CO2 in the atmosphere, may have been significant for high altitude vegetation during times of lower abundance, while lower sea levels of glacial stages potentially opened areas of continental shelf for colonisation during a substantial portion of the Quaternary. Latitudinal variation in orbital forcing has operated on a regional scale. The pace of climate change in the tropics is dominated by precessional oscillations of c. 20 kyr, while the high latitudes of the south are dominated by obliquity oscillations of c. 40 kyr. In particular, seasonal insolation changes forced by precessional oscillations must have had important consequences for the distribution limits of species, with potentially different effects depending on the latitude. The availability of taxa, altitude and human impact, among other events, have locally influenced the environments. Disentangling the different forcing factors of environmental change that operate on different timescales, and understanding the underlying mechanisms leads to considerable challenges for palaeoecologists. The papers in this Special Issue present a selection of palaeoecological studies throughout South America on vegetation changes and other aspects of the environment, providing a window on the possible complexity of the nature of transitions and timings that are potentially available.

Probing spin-orbit-interaction-induced electron dynamics in the carbon atom by multiphoton ionization

We use R-matrix theory with time dependence (RMT) to investigate multiphoton ionization of ground-state atomic carbon with initial orbital magnetic quantum number M_L=0 and M_L=1 at a laser wavelength of 390 nm and peak intensity of 10(14) W/cm(2). Significant differences in ionization yield and ejected-electron momentum distribution are observed between the two values for M_L. We use our theoretical results to model how the spin-orbit interaction affects electron emission along the laser polarization axis. Under the assumption that an initial C atom is prepared at zero time delay with M_L=0, the dynamics with respect to time delay of an ionizing probe pulse modeled by using RMT theory is found to be in good agreement with available experimental data.

ORBITAL MIXING IN CO CHEMISORPTION ON TRANSITION-METAL SURFACES

The chemisorption of CO on metal surfaces is widely accepted as a model for understanding chemical bonding between molecules and solid surfaces, but is nevertheless still a controversial subject. Ab initio total energy calculations using density functional theory with gradient corrections for CO chemisorption on an extended Pd{110} slab yield good agreement with experimental adsorption energies. Examination of the spatial distribution of individual Bloch states demonstrates that the conventional model for CO chemisorption involving charge donation from CO 5 sigma states to metal states and back-donation from metal states into CO 2 pi states is too simplistic, but the computational results provide direct insight into the chemical bonding within the framework of orbital mixing (or hybridisation). The results provide a sound basis for understanding the bonding between molecules and metal surfaces.

Ultrafast excited state dynamics controlling photochemical isomerization of N-methyl-4-[trans-2-(4-pyridyl)ethenyl]pyridinium coordinated to a {Re-1(CO)(3)(2,2 '-bipyridine)} chromophore

Two Multifunctional photoactive complexes [Re(Cl)(CO)(3)-(MeDpe(+))(2)](2+) and [Re(MeDpe(+))(CO)(3)(bpy)](2+) (MeDpe(+) = N-methyl-4-[trans-2-(4-pyridyl)ethenyl]pyridinium, bpy = 2,2'-bipyridine) were synthesized. characterized. and their redox and photonic properties were investigated by cyclic voltammetry: ultraviolet-visible-infrared (UV/Vis/IR) spectroelectrochemistry, stationary UV/Vis and resonance Raman spectroscopy; photolysis; picosecond time-resolved absorption spectroscopy in the visible and infrared regions: and time-resolved resonance Raman spectroscopy. The first reduction step of either complex Occurs at about -1.1 V versus Fc/Fc(+) and is localized at MeDpe(+). Reduction alone does not induce a trans -> cis isomerization of MeDpe(+). [Re(Cl)(CO)(3)(MeDPe(+))(2)](2+) is photostable, while [Re(MeDpe(+))(CO)(3)(bpy)](2+) and free MeDpe(+) isomerize under near-UV irradiation. The lowest excited state of [Re(Cl)(CO)(3)(MeDPe(+))(2)](2+) has been identified as the Re(Cl)(CO)(3) -> MeDpe(+) (MLCT)-M-3 (MLCT = metal-to-ligand charge transfer), decaying directly to the ground state with lifetimes of approximate to 42 (73%) and approximate to 430ps (27%). Optical excitation of [Re(MeDpe(+))(CO)(3)(bpy)](2+) leads to population of Re(CO)(3) -> MeDpe(+) and Re(CO)(3) -> bpy (MLCT)-M-3 states, from which a MeDpe(+) localized intraligand 3 pi pi* excited state ((IL)-I-3) is populated with lifetimes of approximate to 0.6 and approximate to 10 ps, respectively. The 3IL state undergoes a approximate to 21 ps internal rotation, which eventually produces the cis isomer on a much longer timescale. The different excited-state behavior of the two complexes and the absence of thermodynamically favorable interligand electron transfer in excited [Re(MeDpe(+))(CO)(3)(bpy)](2+) reflect the fine energetic balance between excited states of different orbital origin, which can be tuned by subtle Structural variations. The complex [Re(MeDpe+)(CO)(3)(bpy)](2+) emerges as a prototypical, multifunctional species with complementary redox and photonic behavior.