Testing hydrodynamics schemes in galaxy disc simulations

We examine how three fundamentally different numerical hydrodynamics codes follow the evolution of an isothermal galactic disc with an external spiral potential. We compare an adaptive mesh refinement code (RAMSES), a smoothed particle hydrodynamics code (SPHNG), and a volume-discretised meshless code (GIZMO). Using standard refinement criteria, we find that RAMSES produces a disc that is less vertically concentrated and does not reach such high densities as the SPHNG or gizmo runs. The gas surface density in the spiral arms increases at a lower rate for the RAMSES simulations compared to the other codes. There is also a greater degree of substructure in the SPHNG and GIZMOruns and secondary spiral arms are more pronounced. By resolving the Jeans’ length with a greater number of grid cells we achieve more similar results to the Lagrangian codes used in this study. Other alterations to the refinement scheme (adding extra levels of refinement and refining based on local density gradients) are less successful in reducing the disparity between RAMSES and SPHNG/GIZMO. Although more similar, SPHNG displays different density distributions and vertical mass profiles to all modes of gizmo (including the smoothed particle hydrodynamics version). This suggests differences also arise which are not intrinsic to the particular method but rather due to its implementation. The discrepancies between codes (in particular, the densities reached in the spiral arms) could potentially result in differences in the locations and timescales for gravitational collapse, and therefore impact star formation activity in more complex galaxy disc simulations

Simulations of spiral galaxies with an active potential: molecular cloud formation and gas dynamics

We describe simulations of the response of a gaseous disc to an active spiral potential. The potential is derived from an N-body calculation and leads to a multi-armed time-evolving pattern. The gas forms long spiral arms typical of grand design galaxies, although the spiral pattern is asymmetric. The primary difference from a grand-design spiral galaxy, which has a consistent 2/4-armed pattern, is that instead of passing through the spiral arms, gas generally falls into a developing potential minimum and is released only when the local minimum dissolves. In this case, the densest gas is coincident with the spiral potential, rather than offset as in the grand-design spirals. We would there fore expect no offset between the spiral shock and star formation, and no obvious co-rotation radius. Spurs which occur in grand-design spirals when large clumps are sheared off leaving the spiral arms, are rare in the active, time-evolving spiral reported here. Instead, large branches are formed from spiral arms when the underlying spiral potential is dissolving due to the N-body dynamics. We find that the molecular cloud mass spectrum for the active potential is similar to that for clouds in grand design calculations, depending primarily on the ambient pressure rather than the nature of the potential. The largest molecular clouds occur when spiral arms collide, rather than by agglomeration within a spiral arm.Comment: 11 pages, 7 figures, accepted for publication in MNRA

Effects of bracing in adolescents with idiopathic scoliosis

BACKGROUND: The role of bracing in patients with adolescent idiopathic scoliosis who are at risk for curve progression and eventual surgery is controversial. METHODS: We conducted a multicenter study that included patients with typical indications for bracing due to their age, skeletal immaturity, and degree of scoliosis. Both a randomized cohort and a preference cohort were enrolled. Of 242 patients included in the analysis, 116 were randomly assigned to bracing or observation, and 126 chose between bracing and observation. Patients in the bracing group were instructed to wear the brace at least 18 hours per day. The primary outcomes were curve progression to 50 degrees or more (treatment failure) and skeletal maturity without this degree of curve progression (treatment success). RESULTS: The trial was stopped early owing to the efficacy of bracing. In an analysis that included both the randomized and preference cohorts, the rate of treatment success was 72% after bracing, as compared with 48% after observation (propensity-score–adjusted odds ratio for treatment success, 1.93; 95% confidence interval [CI], 1.08 to 3.46). In the intention-to-treat analysis, the rate of treatment success was 75% among patients randomly assigned to bracing, as compared with 42% among those randomly assigned to observation (odds ratio, 4.11; 95% CI, 1.85 to 9.16). There was a significant positive association between hours of brace wear and rate of treatment success (P CONCLUSIONS: Bracing significantly decreased the progression of high-risk curves to the threshold for surgery in patients with adolescent idiopathic scoliosis. The benefit increased with longer hours of brace wear. (Funded by the National Institute of Arthritis and Musculoskeletal and Skin Diseases and others; BRAIST ClinicalTrials.gov number, NCT00448448opens in new tab.)</p

Next-to-leading order QCD calculations with parton showers II: soft singularities

Programs that calculate observables in quantum chromodynamics at next-to-leading order typically generate events that consist of partons rather than hadrons -- and just a few partons at that. These programs would be much more useful if the few partons were turned into parton showers, which could be given to one of the Monte Carlo event generators to produce hadron showers. In a previous paper, we have seen how to generate parton showers related to the final state collinear singularities of the perturbative calculation for the example of e+ + e- --> 3 jets. This paper discusses the treatment of the soft singularities.Comment: 26 pages with 5 figures. This version is close to the version to be publishe

Incorporating next-to-leading order matrix elements for hadronic diboson production in showering event generators

A method for incorporating information from next-to-leading order QCD matrix elements for hadronic diboson production into showering event generators is presented. In the hard central region (high jet transverse momentum) where perturbative QCD is reliable, events are sampled according to the first order tree level matrix element. In the soft and collinear regions next-to-leading order corrections are approximated by calculating the differential cross section across the phase space accessible to the parton shower using the first order (virtual graphs included) matrix element. The parton shower then provides an all-orders exclusive description of parton emissions. Events generated in this way provide a physical result across the entire jet transverse momentum spectrum, have next-to-leading order normalization everywhere, and have positive definite event weights. The method is generalizable without modification to any color singlet production process.Comment: 13 pages, 9 figure

Sex and the preimplantation embryo: implications of sexual dimorphism in the preimplantation period for maternal programming of embryonic development.

201

A Positive-Weight Next-to-Leading-Order Monte Carlo for Z Pair Hadroproduction

We present a first application of a previously published method for the computation of QCD processes that is accurate at next-to-leading order, and that can be interfaced consistently to standard shower Monte Carlo programs. We have considered Z pair production in hadron-hadron collisions, a process whose complexity is sufficient to test the general applicability of the method. We have interfaced our result to the HERWIG and PYTHIA shower Monte Carlo programs. Previous work on next-to-leading order corrections in a shower Monte Carlo (the MC@NLO program) may involve the generation of events with negative weights, that are avoided with the present method. We have compared our results with those obtained with MC@NLO, and found remarkable consistency. Our method can also be used as a standalone, alternative implementation of QCD corrections, with the advantage of positivity, improved convergence, and next-to-leading logarithmic accuracy in the region of small transverse momentum of the radiated parton.Comment: 33 pages, 10 figure

Classifying the biodiversity of the Great Barrier Reef World Heritage Area for the classification phase of the representative areas program

This technical report outlines the methods that the Great Barrier Reef Marine Park Authority used to classify the biodiversity of the marine environs of the Great Barrier Reef World Heritage Area for the Representative Areas Program. Classification was the first step in the multiphase Representative Areas Program that eventuated in a new network of no-take areas, free from extractive activities, in the Great Barrier Reef Marine Park

The properties, origin and evolution of stellar clusters in galaxy simulations and observations

We investigate the properties and evolution of star particles in two simulations of isolated spiral galaxies, and two galaxies from cosmological simulations. Unlike previous numerical work, where typically each star particle represents one ‘cluster’, for the isolated galaxies we are able to model features we term ‘clusters’ with groups of particles. We compute the spatial distribution of stars with different ages, and cluster mass distributions, comparing our findings with observations including the recent LEGUS survey. We find that spiral structure tends to be present in older (100s Myrs) stars and clusters in the simulations compared to the observations. This likely reflects differences in the numbers of stars or clusters, the strength of spiral arms, and whether the clusters are allowed to evolve. Where we model clusters with multiple particles, we are able to study their evolution. The evolution of simulated clusters tends to follow that of their natal gas clouds. Massive, dense, long-lived clouds host massive clusters, whilst short-lived clouds host smaller clusters which readily disperse. Most clusters appear to disperse fairly quickly, in basic agreement with observational findings. We note that embedded clusters may be less inclined to disperse in simulations in a galactic environment with continuous accretion of gas onto the clouds than isolated clouds and correspondingly, massive young clusters which are no longer associated with gas tend not to occur in the simulations. Caveats of our models include that the cluster densities are lower than realistic clusters, and the simplistic implementation of stellar feedback