A Direct Elliptic Solver Based on Hierarchically Low-rank Schur Complements

A parallel fast direct solver for rank-compressible block tridiagonal linear systems is presented. Algorithmic synergies between Cyclic Reduction and Hierarchical matrix arithmetic operations result in a solver with $O(N \log^2 N)$ arithmetic complexity and $O(N \log N)$ memory footprint. We provide a baseline for performance and applicability by comparing with well known implementations of the $\mathcal{H}$-LU factorization and algebraic multigrid with a parallel implementation that leverages the concurrency features of the method. Numerical experiments reveal that this method is comparable with other fast direct solvers based on Hierarchical Matrices such as $\mathcal{H}$-LU and that it can tackle problems where algebraic multigrid fails to converge

Simulating coronal condensation dynamics in 3D

We present numerical simulations in 3D settings where coronal rain phenomena take place in a magnetic configuration of a quadrupolar arcade system. Our simulation is a magnetohydrodynamic simulation including anisotropic thermal conduction, optically thin radiative losses, and parametrised heating as main thermodynamical features to construct a realistic arcade configuration from chromospheric to coronal heights. The plasma evaporation from chromospheric and transition region heights eventually causes localised runaway condensation events and we witness the formation of plasma blobs due to thermal instability, that evolve dynamically in the heated arcade part and move gradually downwards due to interchange type dynamics. Unlike earlier 2.5D simulations, in this case there is no large scale prominence formation observed, but a continuous coronal rain develops which shows clear indications of Rayleigh-Taylor or interchange instability, that causes the denser plasma located above the transition region to fall down, as the system moves towards a more stable state. Linear stability analysis is used in the non-linear regime for gaining insight and giving a prediction of the system's evolution. After the plasma blobs descend through interchange, they follow the magnetic field topology more closely in the lower coronal regions, where they are guided by the magnetic dips.Comment: 47 pages, 59 figure

Parametric survey of longitudinal prominence oscillation simulations

It is found that both microflare-sized impulsive heating at one leg of the loop and a suddenly imposed velocity perturbation can propel the prominence to oscillate along the magnetic dip. An extensive parameter survey results in a scaling law, showing that the period of the oscillation, which weakly depends on the length and height of the prominence, and the amplitude of the perturbations, scales with $\sqrt{R/g_\odot}$, where $R$ represents the curvature radius of the dip, and $g_\odot$ is the gravitational acceleration of the Sun. This is consistent with the linear theory of a pendulum, which implies that the field-aligned component of gravity is the main restoring force for the prominence longitudinal oscillations, as confirmed by the force analysis. However, the gas pressure gradient becomes non-negligible for short prominences. The oscillation damps with time in the presence of non-adiabatic processes. Compared to heat conduction, the radiative cooling is the dominant factor leading to the damping. A scaling law for the damping timescale is derived, i.e., $\tau\sim l^{1.63} D^{0.66}w^{-1.21}v_{0}^{-0.30}$, showing strong dependence on the prominence length $l$, the geometry of the magnetic dip (characterized by the depth $D$ and the width $w$), and the velocity perturbation amplitude $v_0$. The larger the amplitude, the faster the oscillation damps. It is also found that mass drainage significantly reduces the damping timescale when the perturbation is too strong.Comment: 17 PAGES, 8FIGURE

Spectroscopy of reflection-asymmetric nuclei with relativistic energy density functionals

Quadrupole and octupole deformation energy surfaces, low-energy excitation spectra and transition rates in fourteen isotopic chains: Xe, Ba, Ce, Nd, Sm, Gd, Rn, Ra, Th, U, Pu, Cm, Cf, and Fm, are systematically analyzed using a theoretical framework based on a quadrupole-octupole collective Hamiltonian (QOCH), with parameters determined by constrained reflection-asymmetric and axially-symmetric relativistic mean-field calculations. The microscopic QOCH model based on the PC-PK1 energy density functional and $\delta$-interaction pairing is shown to accurately describe the empirical trend of low-energy quadrupole and octupole collective states, and predicted spectroscopic properties are consistent with recent microscopic calculations based on both relativistic and non-relativistic energy density functionals. Low-energy negative-parity bands, average octupole deformations, and transition rates show evidence for octupole collectivity in both mass regions, for which a microscopic mechanism is discussed in terms of evolution of single-nucleon orbitals with deformation.Comment: 36 pages, 21 figures, Accepted for Publication in Physical Review

Characterization of the Soluble Nanoparticles Formed through Coulombic Interaction of Bovine Serum Albumin with Anionic Graft Copolymers at Low pH

A static light scattering (SLS) study of bovine serum albumin (BSA) mixtures with two anionic graft copolymers of poly (sodium acrylate-co-sodium 2-acrylamido-2-methyl-1-propanesulphonate)-graft-poly (N, N-dimethylacrylamide), with a high composition in poly (N, N-dimethylacrylamide) (PDMAM) side chains, revealed the formation of oppositely charged complexes, at pH lower than 4.9, the isoelectric point of BSA. The core-corona nanoparticles formed at pH = 3.00, were characterized. Their molecular weight and radius of gyration were determined by SLS, while their hydrodynamic radius was determined by dynamic light scattering. Small angle neutron scattering measurements were used to determine the radius of the insoluble complexes, comprising the core of the particles. The values obtained indicated that their size and aggregation number of the nanoparticles, were smaller when the content of the graft copolymers in neutral PDMAM side chains was higher. Such particles should be interesting drug delivery candidates, if the gastrointestinal tract was to be used