Coupling of Coronal and Heliospheric Magnetohydrodynamic Models: Solution Comparisons and Verification

Two well-established magnetohydrodynamic (MHD) codes are coupled to model the solar corona and the inner heliosphere. The corona is simulated using the MHD algorithm outside a sphere (MAS) model. The Lyon–Fedder–Mobarry (LFM) model is used in the heliosphere. The interface between the models is placed in a spherical shell above the critical point and allows both models to work in either a rotating or an inertial frame. Numerical tests are presented examining the coupled model solutions from 20 to 50 solar radii. The heliospheric simulations are run with both LFM and the MAS extension into the heliosphere, and use the same polytropic coronal MAS solutions as the inner boundary condition. The coronal simulations are performed for idealized magnetic configurations, with an out-of-equilibrium flux rope inserted into an axisymmetric background, with and without including the solar rotation. The temporal evolution at the inner boundary of the LFM and MAS solutions is shown to be nearly identical, as are the steady-state background solutions, prior to the insertion of the flux rope. However, after the coronal mass ejection has propagated through the significant portion of the simulation domain, the heliospheric solutions diverge. Additional simulations with different resolution are then performed and show that the MAS heliospheric solutions approach those of LFM when run with progressively higher resolution. Following these detailed tests, a more realistic simulation driven by the thermodynamic coronal MAS is presented, which includes solar rotation and an azimuthally asymmetric background and extends to the Earth’s orbit

Nonlinear force-free reconstruction of the global solar magnetic field: methodology

We present a novel numerical method that allows the calculation of nonlinear force-free magnetostatic solutions above a boundary surface on which only the distribution of the normal magnetic field component is given. The method relies on the theory of force-free electrodynamics and applies directly to the reconstruction of the solar coronal magnetic field for a given distribution of the photospheric radial field component. The method works as follows: we start with any initial magnetostatic global field configuration (e.g. zero, dipole), and along the boundary surface we create an evolving distribution of tangential (horizontal) electric fields that, via Faraday's equation, give rise to a respective normal field distribution approaching asymptotically the target distribution. At the same time, these electric fields are used as boundary condition to numerically evolve the resulting electromagnetic field above the boundary surface, modelled as a thin ideal plasma with non-reflecting, perfectly absorbing outer boundaries. The simulation relaxes to a nonlinear force-free configuration that satisfies the given normal field distribution on the boundary. This is different from existing methods relying on a fixed boundary condition - the boundary evolves toward the a priori given one, at the same time evolving the three-dimensional field solution above it. Moreover, this is the first time a nonlinear force-free solution is reached by using only the normal field component on the boundary. This solution is not unique, but depends on the initial magnetic field configuration and on the evolutionary course along the boundary surface. To our knowledge, this is the first time that the formalism of force-free electrodynamics, used very successfully in other astrophysical contexts, is applied to the global solar magnetic field.Comment: 18 pages, 5 figures, Solar Physic

Revisiting the origin of the bending in group 2 metallocenes AeCp2 (Ae = Be–Ba)

Metallocenes are well-established compounds in organometallic chemistry, and can exhibit either a coplanar structure or a bent structure according to the nature of the metal center (E) and the cyclopentadienyl ligands (Cp). Herein, we re-examine the chemical bonding to underline the origins of the geometry and stability observed experimentally. To this end, we have analysed a series of group 2 metallocenes [Ae(C5R5)2] (Ae = Be–Ba and R = H, Me, F, Cl, Br, and I) with a combination of computational methods, namely energy decomposition analysis (EDA), polarizability model (PM), and dispersion interaction densities (DIDs). Although the metal–ligand bonding nature is mainly an electrostatic interaction (65–78%), the covalent character is not negligible (33–22%). Notably, the heavier the metal center, the stronger the d-orbital interaction with a 50% contribution to the total covalent interaction. The dispersion interaction between the Cp ligands counts only for 1% of the interaction. Despite that orbital contributions become stronger for heavier metals, they never represent the energy main term. Instead, given the electrostatic nature of the metallocene bonds, we propose a model based on polarizability, which faithfully predicts the bending angle. Although dispersion interactions have a fair contribution to strengthen the bending angle, the polarizability plays a major role

Neuroanatomic Correlates of Female Sexual Dysfunction in Multiple Sclerosis

OBJECTIVE: This study intended to determine associations between alterations of female sexual arousal as well as vaginal lubrication and the site of cerebral multiple sclerosis (MS) lesions. METHODS: In 44 women with MS (mean age: 36.5 ± 9.9 years), we assessed their medical history and evaluated sexual function using the Female Sexual Function Index scores for arousal and vaginal lubrication. We determined potential confounding factors of sexual dysfunction: age; disease duration; physical disability; depression; bladder or urinary dysfunction; and total volume of cerebral lesions. Arousal and lubrication scores were correlated with one another and with potential confounding factors. Cerebral MS lesions were recorded on imaging scans. A voxel-based lesion symptom mapping (VLSM) analysis adjusted for confounding variables was performed correlating cerebral sites of MS lesions with arousal and lubrication scores. RESULTS: Decreased arousal scores correlated with decreased lubrication scores; decreased lubrication scores were associated with bladder or urinary symptoms. Arousal and lubrication scores were not associated with any other variables. Multivariate VLSM analysis, including arousal and lubrication scores as covariables of interest, showed right occipital lesions associated with impaired arousal and left insular lesions associated with decreased lubrication. Impaired lubrication remained associated with left insular lesions after adjustment for bladder or urinary dysfunction. INTERPRETATION: Our data indicate that impaired female sexual arousal is associated with MS lesions in the occipital region, integrating visual information and modulating attention toward visual input. Impaired lubrication correlated with lesions in the left insular region, contributing to mapping and generating visceral arousal states