Most molecular clouds are filamentary or elongated. For those forming low-mass stars (<8 solar masses), the competition between self-gravity and turbulent pressure along the dynamically dominant intercloud magnetic field (10 to 100 parsecs) shapes the clouds to be elongated either perpendicularly or parallel to the fields. A recent study also suggested that on the scales of 0.1 to 0.01 parsecs, such fields are dynamically important within cloud cores forming massive stars (>8 solar masses). But whether the core field morphologies are inherited from the intercloud medium or governed by cloud turbulence is unknown, as is the effect of magnetic fields on cloud fragmentation at scales of 10 to 0.1 parsecs. Here we report magnetic-field maps inferred from polarimetric observations of NGC 6334, a region forming massive stars, on the 100 to 0.01 parsec scale. NGC 6334 hosts young star-forming sites where fields are not severely affected by stellar feedback, and their directions do not change much over the entire scale range. This means that the fields are dynamically important. The ordered fields lead to a self-similar gas fragmentation: at all scales, there exist elongated gas structures nearly perpendicular to the fields. Many gas elongations have density peaks near the ends, which symmetrically pinch the fields. The field strength is proportional to the 0.4th power of the density, which is an indication of anisotropic gas contractions along the field. We conclude that magnetic fields have a crucial role in the fragmentation of NGC 6334.

A full-dimensional quantum dynamical study of the rotational excitation in para-para H2 + H2 collisions using the potential-energy surface of Boothroyd et al. [J. Chem. Phys. 116, 666 (2002)] is reported. The multiconfiguration time-dependent Hartree algorithm is utilized to propagate wave packets and the cross sections for collision energies up to 1.2 eV are determined by a flux analysis through the interaction of the wave packet with a complex absorbing potential. Calculations for a collection of total angular momenta up to J = 70 are performed; the missing channels are obtained with a J-interpolation algorithm.

The multi-layer multi-configuration time-dependent Hartree method (ML-MCTDH) is a highly efficient scheme for studying the dynamics of high-dimensional quantum systems. Its use is greatly facilitated if the Hamiltonian of the system possesses a particular structure through which the multi-dimensional matrix elements can be computed efficiently. In the field of quantum molecular dynamics, the effective interaction between the atoms is often described by potential energy surfaces (PES), and it is necessary to fit such PES into the desired structure. For high-dimensional systems, the current approaches for this fitting process either lead to fits that are too large to be practical, or their accuracy is difficult to predict and control.This article introduces multi-layer Potfit (MLPF), a novel fitting scheme that results in a PES representation in the hierarchical tensor (HT) format. The scheme is based on the hierarchical singular value decomposition, which can yield a near-optimal fit and give strict bounds for the obtained accuracy. Here, a recursive scheme for using the HT-format PES within ML-MCTDH is derived, and theoretical estimates as well as a computational example show that the use of MLPF can reduce the numerical effort for ML-MCTDH by orders of magnitude, compared to the traditionally used Potfit representation of the PES. Moreover, it is shown that MLPF is especially beneficial for high-accuracy PES representations, and it turns out that MLPF leads to computational savings already for comparatively small systems with just four modes.

Articles you may be interested inRotational excitations in collisions: Full-and reduced-dimensional quantum wave packet studies comparing different potential energy surfaces

We study the process of rotational excitation in the collisions of para-H(2) with para-H(2) by propagating wave packets with the multiconfiguration time-dependent Hartree (MCTDH) algorithm. Transition probabilities are then calculated by the method of Tannor and Weeks based on time-correlation functions. Calculations were carried out up to a total angular momentum of J=70 to compute integral cross sections up to 1.2 eV in collision energy and thermal rate coefficients from 100 to 3000 K. The process is studied on the full-dimensional potential energy surface of Boothroyd-Martin-Keogh-Peterson (BMKP) as well as on the rigid rotor surface of Diep and Johnson. We test the validity of the rigid rotor approximation by also considering two rigid rotor restrictions of the BMKP potential energy surface (PES). Additionally, we investigate a variant of the BMKP PES suggested by Pogrebnya and Clary [Chem. Phys. Lett. 363, 523 (2002)] with reduced anisotropy. We compare our results with previous theoretical data for the cross sections and with experimental data for the rate coefficients at low temperatures.

We present the results of a full-dimensional quantum mechanical study of the rovibrational energy transfer in the collision between ortho-H2 and para-H2 in the energy range of 0.1-1.0 eV. The multiconfiguration time-dependent Hartree algorithm has been used to propagate the wave packets on the global potential energy surface by Boothroyd et al. [J. Chem. Phys. 116, 666 (2002)] and on a modification of this surface where the short range anisotropy is reduced. State-to-state attributes such as probabilities or integral cross sections are obtained using the formalism of Tannor and Weeks [J. Chem. Phys. 98, 3884 (1993)] by Fourier transforming the correlation functions. The effect of initial rotation of the diatoms on the inelastic and de-excitation processes is investigated.

We investigate the exact wavefunction as a single product of electronic and nuclear wavefunction for a model conical intersection system. Exact factorized spiky potentials and nodeless nuclear wavefunctions are found. The exact factorized potential preserves the symmetry breaking effect when the coupling mode is present. Additionally the nodeless wavefunctions are found to be closely related to the adiabatic nuclear eigenfunctions. This phenomenon holds even for the regime where the non-adiabatic coupling is relevant, and sheds light on the relation between the exact wavefunction factorization and the adiabatic approximation. 1 arXiv:1311.4459v1 [quant-ph]

The complex interplay between turbulence, magnetic fields, and self-gravity leads to the formation of molecular clouds out of the diffuse interstellar medium (ISM). One avenue of studying this interplay is by analyzing statistical features derived from observations, where the interpretation of these features is greatly facilitated by comparisons with numerical simulations. Here we focus on the statistical anisotropy present in synthetic maps of velocity centroid data, which we derive from three-dimensional magnetohydrodynamic simulations of a turbulent, magnetized, self-gravitating patch of ISM. We study how the orientation and magnitude of the velocity anisotropy correlate with the magnetic field and with the structures generated by gravitational collapse. Motivated by recent observational constraints, our simulations focus on the supersonic (sonic Mach number M ≈ 2 − 17) but sub-to trans-alfvénic (alfvénic Mach number M A ≈ 0.2 − 1.2) turbulence regime, and we consider clouds which are barely to mildly magnetically supercritical (mass-to-flux ratio equal to once or twice the critical value). Additionally we explore the impact of the turbulence driving mechanism (solenoidal or compressive) on the velocity anisotropy. While we confirm previous findings that the velocity anisotropy generally aligns well with the plane-of-sky magnetic field, our inclusion of the effects of self-gravity reveals that in regions of higher column density, the velocity anisotropy may be destroyed or even reoriented to align with the gravitationally formed structures. We provide evidence that this effect is not necessarily due to the increase of M A inside the high-density regions.