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.