Mechanical Systems with Symmetry, Variational Principles, and Integration Algorithms

This paper studies variational principles for mechanical systems with symmetry and their applications to integration algorithms. We recall some general features of how to reduce variational principles in the presence of a symmetry group along with general features of integration algorithms for mechanical systems. Then we describe some integration algorithms based directly on variational principles using a discretization technique of Veselov. The general idea for these variational integrators is to directly discretize Hamilton’s principle rather than the equations of motion in a way that preserves the original systems invariants, notably the symplectic form and, via a discrete version of Noether’s theorem, the momentum map. The resulting mechanical integrators are second-order accurate, implicit, symplectic-momentum algorithms. We apply these integrators to the rigid body and the double spherical pendulum to show that the techniques are competitive with existing integrators

Tackling Exascale Software Challenges in Molecular Dynamics Simulations with GROMACS

GROMACS is a widely used package for biomolecular simulation, and over the last two decades it has evolved from small-scale efficiency to advanced heterogeneous acceleration and multi-level parallelism targeting some of the largest supercomputers in the world. Here, we describe some of the ways we have been able to realize this through the use of parallelization on all levels, combined with a constant focus on absolute performance. Release 4.6 of GROMACS uses SIMD acceleration on a wide range of architectures, GPU offloading acceleration, and both OpenMP and MPI parallelism within and between nodes, respectively. The recent work on acceleration made it necessary to revisit the fundamental algorithms of molecular simulation, including the concept of neighborsearching, and we discuss the present and future challenges we see for exascale simulation - in particular a very fine-grained task parallelism. We also discuss the software management, code peer review and continuous integration testing required for a project of this complexity.Comment: EASC 2014 conference proceedin

Extracellular Matrix Aggregates from Differentiating Embryoid Bodies as a Scaffold to Support ESC Proliferation and Differentiation

Embryonic stem cells (ESCs) have emerged as potential cell sources for tissue engineering and regeneration owing to its virtually unlimited replicative capacity and the potential to differentiate into a variety of cell types. Current differentiation strategies primarily involve various growth factor/inducer/repressor concoctions with less emphasis on the substrate. Developing biomaterials to promote stem cell proliferation and differentiation could aid in the realization of this goal. Extracellular matrix (ECM) components are important physiological regulators, and can provide cues to direct ESC expansion and differentiation. ECM undergoes constant remodeling with surrounding cells to accommodate specific developmental event. In this study, using ESC derived aggregates called embryoid bodies (EB) as a model, we characterized the biological nature of ECM in EB after exposure to different treatments: spontaneously differentiated and retinoic acid treated (denoted as SPT and RA, respectively). Next, we extracted this treatment-specific ECM by detergent decellularization methods (Triton X-100, DOC and SDS are compared). The resulting EB ECM scaffolds were seeded with undifferentiated ESCs using a novel cell seeding strategy, and the behavior of ESCs was studied. Our results showed that the optimized protocol efficiently removes cells while retaining crucial ECM and biochemical components. Decellularized ECM from SPT EB gave rise to a more favorable microenvironment for promoting ESC attachment, proliferation, and early differentiation, compared to native EB and decellularized ECM from RA EB. These findings suggest that various treatment conditions allow the formulation of unique ESC-ECM derived scaffolds to enhance ESC bioactivities, including proliferation and differentiation for tissue regeneration applications. © 2013 Goh et al

Thermal Behavior of Benzoic Acid/lsonicotinamide Binary Cocrystals

A comprehensive study of the thermal behavior of the 1:1 and 2:1 benzoic acid/isonicotinamide cocrystals is reported. The 1:1 material shows a simple unit cell expansion followed by melting upon heating. The 2:1 crystal exhibits more complex behavior. Its unit cell first expands upon heating, as a result of C–H···π interactions being lengthened. It then is converted into the 1:1 crystal, as demonstrated by significant changes in its X-ray diffraction pattern. The loss of 1 equiv of benzoic acid is confirmed by thermogravimetric analysis–mass spectrometry. Hot stage microscopy confirms that, as intuitively expected, the transformation begins at the crystal surface. The temperature at which conversion occurs is highly dependent on the sample mass and geometry, being reduced when the sample is under a gas flow or has a greater exposed surface area but increased when the heating rate is elevated

A meta-analytic review of stand-alone interventions to improve body image

Objective Numerous stand-alone interventions to improve body image have been developed. The present review used meta-analysis to estimate the effectiveness of such interventions, and to identify the specific change techniques that lead to improvement in body image. Methods The inclusion criteria were that (a) the intervention was stand-alone (i.e., solely focused on improving body image), (b) a control group was used, (c) participants were randomly assigned to conditions, and (d) at least one pretest and one posttest measure of body image was taken. Effect sizes were meta-analysed and moderator analyses were conducted. A taxonomy of 48 change techniques used in interventions targeted at body image was developed; all interventions were coded using this taxonomy. Results The literature search identified 62 tests of interventions (N = 3,846). Interventions produced a small-to-medium improvement in body image (d+ = 0.38), a small-to-medium reduction in beauty ideal internalisation (d+ = -0.37), and a large reduction in social comparison tendencies (d+ = -0.72). However, the effect size for body image was inflated by bias both within and across studies, and was reliable but of small magnitude once corrections for bias were applied. Effect sizes for the other outcomes were no longer reliable once corrections for bias were applied. Several features of the sample, intervention, and methodology moderated intervention effects. Twelve change techniques were associated with improvements in body image, and three techniques were contra-indicated. Conclusions The findings show that interventions engender only small improvements in body image, and underline the need for large-scale, high-quality trials in this area. The review identifies effective techniques that could be deployed in future interventions

Search for intermediate mass black hole binaries in the first and second observing runs of the Advanced LIGO and Virgo network

Gravitational-wave astronomy has been firmly established with the detection of gravitational waves from the merger of ten stellar-mass binary black holes and a neutron star binary. This paper reports on the all-sky search for gravitational waves from intermediate mass black hole binaries in the first and second observing runs of the Advanced LIGO and Virgo network. The search uses three independent algorithms: two based on matched filtering of the data with waveform templates of gravitational-wave signals from compact binaries, and a third, model-independent algorithm that employs no signal model for the incoming signal. No intermediate mass black hole binary event is detected in this search. Consequently, we place upper limits on the merger rate density for a family of intermediate mass black hole binaries. In particular, we choose sources with total masses M=m1+m2ϵ[120,800] M and mass ratios q=m2/m1ϵ[0.1,1.0]. For the first time, this calculation is done using numerical relativity waveforms (which include higher modes) as models of the real emitted signal. We place a most stringent upper limit of 0.20 Gpc-3 yr-1 (in comoving units at the 90% confidence level) for equal-mass binaries with individual masses m1,2=100 M and dimensionless spins χ1,2=0.8 aligned with the orbital angular momentum of the binary. This improves by a factor of ∼5 that reported after Advanced LIGO's first observing run

Search for Eccentric Binary Black Hole Mergers with Advanced LIGO and Advanced Virgo during Their First and Second Observing Runs

When formed through dynamical interactions, stellar-mass binary black holes (BBHs) may retain eccentric orbits (e > 0.1 at 10 Hz) detectable by ground-based gravitational-wave detectors. Eccentricity can therefore be used to differentiate dynamically formed binaries from isolated BBH mergers. Current template-based gravitational-wave searches do not use waveform models associated with eccentric orbits, rendering the search less efficient for eccentric binary systems. Here we present the results of a search for BBH mergers that inspiral in eccentric orbits using data from the first and second observing runs (O1 and O2) of Advanced LIGO and Advanced Virgo. We carried out the search with the coherent WaveBurst algorithm, which uses minimal assumptions on the signal morphology and does not rely on binary waveform templates. We show that it is sensitive to binary mergers with a detection range that is weakly dependent on eccentricity for all bound systems. Our search did not identify any new binary merger candidates. We interpret these results in light of eccentric binary formation models. We rule out formation channels with rates ⪆100 Gpc-3 yr-1 for e > 0.1, assuming a black hole mass spectrum with a power-law index ≲2

All-sky search for short gravitational-wave bursts in the second Advanced LIGO and Advanced Virgo run

We present the results of a search for short-duration gravitational-wave transients in the data from the second observing run of Advanced LIGO and Advanced Virgo. We search for gravitational-wave transients with a duration of milliseconds to approximately one second in the 32-4096 Hz frequency band with minimal assumptions about the signal properties, thus targeting a wide variety of sources. We also perform a matched-filter search for gravitational-wave transients from cosmic string cusps for which the waveform is well modeled. The unmodeled search detected gravitational waves from several binary black hole mergers which have been identified by previous analyses. No other significant events have been found by either the unmodeled search or the cosmic string search. We thus present the search sensitivities for a variety of signal waveforms and report upper limits on the source rate density as a function of the characteristic frequency of the signal. These upper limits are a factor of 3 lower than the first observing run, with a 50% detection probability for gravitational-wave emissions with energies of ∼10-9 Mc2 at 153 Hz. For the search dedicated to cosmic string cusps we consider several loop distribution models, and present updated constraints from the same search done in the first observing run

Search for gravitational waves from Scorpius X-1 in the second Advanced LIGO observing run with an improved hidden Markov model

We present results from a semicoherent search for continuous gravitational waves from the low-mass X-ray binary Scorpius X-1, using a hidden Markov model (HMM) to track spin wandering. This search improves on previous HMM-based searches of LIGO data by using an improved frequency domain matched filter, the $\mathcal{J}$-statistic, and by analysing data from Advanced LIGO's second observing run. In the frequency range searched, from $60$ to $650\,\mathrm{Hz}$, we find no evidence of gravitational radiation. At $194.6\,\mathrm{Hz}$, the most sensitive search frequency, we report an upper limit on gravitational wave strain (at 95\% confidence) of $h_0^{95\%} = 3.47 \times 10^{-25}$ when marginalising over source inclination angle. This is the most sensitive search for Scorpius X-1, to date, that is specifically designed to be robust in the presence of spin wandering