Higher order Schrodinger and Hartree-Fock equations

The domain of validity of the higher-order Schrodinger equations is analyzed for harmonic-oscillator and Coulomb potentials as typical examples. Then the Cauchy theory for higher-order Hartree-Fock equations with bounded and Coulomb potentials is developed. Finally, the existence of associated ground states for the odd-order equations is proved. This renders these quantum equations relevant for physics.Comment: 19 pages, to appear in J. Math. Phy

Models of Distorted and Evolving Dark Matter Halos

We investigate the ability of basis function expansions to reproduce the evolution of a Milky Way-like dark matter halo, extracted from a cosmological zoom-in simulation. For each snapshot, the density of the halo is reduced to a basis function expansion, with interpolation used to recreate the evolution between snapshots. The angular variation of the halo density is described by spherical harmonics, and the radial variation either by biorthonormal basis functions adapted to handle truncated haloes or by splines. High fidelity orbit reconstructions are attainable using either method with similar computational expense. We quantify how the error in the reconstructed orbits varies with expansion order and snapshot spacing. Despite the many possible biorthonormal expansions, it is hard to beat a conventional Hernquist-Ostriker expansion with a moderate number of terms ($\gtrsim15$ radial and $\gtrsim6$ angular). As two applications of the developed machinery, we assess the impact of the time-dependence of the potential on (i) the orbits of Milky Way satellites, and (ii) planes of satellites as observed in the Milky Way and other nearby galaxies. Time evolution over the last 5 Gyr introduces an uncertainty in the Milky Way satellites' orbital parameters of $\sim 15$ per cent, comparable to that induced by the observational errors or the uncertainty in the present-day Milky Way potential. On average, planes of satellites grow at similar rates in evolving and time-independent potentials. There can be more, or less, growth in the plane's thickness, if the plane becomes less, or more, aligned with the major or minor axis of the evolving halo.Comment: MNRAS, submitte

DNA cruciform arms nucleate through a correlated but non-synchronous cooperative mechanism

Inverted repeat (IR) sequences in DNA can form non-canonical cruciform structures to relieve torsional stress. We use Monte Carlo simulations of a recently developed coarse-grained model of DNA to demonstrate that the nucleation of a cruciform can proceed through a cooperative mechanism. Firstly, a twist-induced denaturation bubble must diffuse so that its midpoint is near the centre of symmetry of the IR sequence. Secondly, bubble fluctuations must be large enough to allow one of the arms to form a small number of hairpin bonds. Once the first arm is partially formed, the second arm can rapidly grow to a similar size. Because bubbles can twist back on themselves, they need considerably fewer bases to resolve torsional stress than the final cruciform state does. The initially stabilised cruciform therefore continues to grow, which typically proceeds synchronously, reminiscent of the S-type mechanism of cruciform formation. By using umbrella sampling techniques we calculate, for different temperatures and superhelical densities, the free energy as a function of the number of bonds in each cruciform along the correlated but non-synchronous nucleation pathways we observed in direct simulations.Comment: 12 pages main paper + 11 pages supplementary dat

Quark Masses: An Environmental Impact Statement

We investigate worlds that lie on a slice through the parameter space of the Standard Model over which quark masses vary. We allow as many as three quarks to participate in nuclei, while fixing the mass of the electron and the average mass of the lightest baryon flavor multiplet. We classify as "congenial" worlds that satisfy the environmental constraint that the quark masses allow for stable nuclei with charges one, six, and eight, making organic chemistry possible. Whether a congenial world actually produces observers depends on a multitude of historical contingencies, beginning with primordial nucleosynthesis, which we do not explore. Such constraints may be independently superimposed on our results. Environmental constraints such as the ones we study may be combined with information about the a priori distribution of quark masses over the landscape of possible universes to determine whether the measured values of the quark masses are determined environmentally, but our analysis is independent of such an anthropic approach. We estimate baryon masses as functions of quark masses and nuclear masses as functions of baryon masses. We check for the stability of nuclei against fission, strong particle emission, and weak nucleon emission. For two light quarks with charges 2/3 and -1/3, we find a band of congeniality roughly 29 MeV wide in their mass difference. We also find another, less robust region of congeniality with one light, charge -1/3 quark, and two heavier, approximately degenerate charge -1/3 and 2/3 quarks. No other assignment of light quark charges yields congenial worlds with two baryons participating in nuclei. We identify and discuss the region in quark-mass space where nuclei would be made from three or more baryon species.Comment: 40 pages, 16 figures (in color), 4 tables. See paper for a more detailed abstract. v4: Cleaning up minor typo

On separating propagating and non-propagating dynamics in fluid-flow equations

The ability to separate acoustically radiating and non-radiating components in fluid flow is desirable to identify the true sources of aerodynamic sound, which can be expressed in terms of the non-radiating flow dynamics. These non-radiating components are obtained by filtering the flow field. Two linear filtering strategies are investigated: one is based on a differential operator, the other employs convolution operations. Convolution filters are found to be superior at separating radiating and non-radiating components. Their ability to decompose the flow into non-radiating and radiating components is demonstrated on two different flows: one satisfying the linearized Euler and the other the Navier-Stokes equations. In the latter case, the corresponding sound sources are computed. These sources provide good insight into the sound generation process. For source localization, they are found to be superior to the commonly used sound sources computed using the steady part of the flow

Geomagnetic induction and conductive structures in north-west India

Magnetic disturbance events and quiet daily variation as recorded by the 1979 magnetometer array study in north-west India are analysed for evidence of electrical conductivity structures in the region. Contour maps of Fourier transform parameters are presented, and the disturbance event data are also reduced to sets of real and quadrature Parkinson arrows over a range of periods. A variety of conductive structures in the area are mapped, including some relatively shallow ones thought to be caused by sediments, as in the Ganga basin. More information is obtained on a major conductivity structure which strikes perpendicular to the Ganga basin into the foothills of the Himalayas; a second major conductivity structure is detected to lie to the west of the array area, and may be associated there with some aspect of the suture zone of India and Asia

Non-Equilibrium Reaction Rates in the Macroscopic Chemistry Method for DSMC Calculations

The Direct Simulation Monte Carlo (DSMC) method is used to simulate the flow of rarefied gases. In the Macroscopic Chemistry Method (MCM) for DSMC, chemical reaction rates calculated from local macroscopic flow properties are enforced in each cell. Unlike the standard total collision energy (TCE) chemistry model for DSMC, the new method is not restricted to an Arrhenius form of the reaction rate coefficient, nor is it restricted to a collision cross-section which yields a simple power-law viscosity. For reaction rates of interest in aerospace applications, chemically reacting collisions are generally infrequent events and, as such, local equilibrium conditions are established before a significant number of chemical reactions occur. Hence, the reaction rates which have been used in MCM have been calculated from the reaction rate data which are expected to be correct only for conditions of thermal equilibrium. Here we consider artificially high reaction rates so that the fraction of reacting collisions is not small and propose a simple method of estimating the rates of chemical reactions which can be used in the Macroscopic Chemistry Method in both equilibrium and non-equilibrium conditions. Two tests are presented: (1) The dissociation rates under conditions of thermal non-equilibrium are determined from a zero-dimensional Monte-Carlo sampling procedure which simulates ‘intra-modal’ non-equilibrium; that is, equilibrium distributions in each of the translational, rotational and vibrational modes but with different temperatures for each mode; (2) The 2-D hypersonic flow of molecular oxygen over a vertical plate at Mach 30 is calculated. In both cases the new method produces results in close agreement with those given by the standard TCE model in the same highly nonequilibrium conditions. We conclude that the general method of estimating the non-equilibrium reaction rate is a simple means by which information contained within non-equilibrium distribution functions predicted by the DSMC method can be included in the Macroscopic Chemistry Method

A complex ray-tracing tool for high-frequency mean-field flow interaction effects in jets

This paper presents a complex ray-tracing tool for the calculation of high-frequency Green’s functions in 3D mean field jet flows. For a generic problem, the ray solution suffers from three main deficiencies: multiplicity of solutions, singularities at caustics, and the determining of complex solutions. The purpose of this paper is to generalize, combine and apply existing stationary media methods to moving media scenarios. Multiplicities are dealt with using an equivalent two-point boundary-value problem, whilst non-uniformities at caustics are corrected using diffraction catastrophes. Complex rays are found using a combination of imaginary perturbations, an assumption of caustic stability, and analytic continuation of the receiver curve. To demonstrate this method, the ray tool is compared against a high-frequency modal solution of Lilley’s equation for an off-axis point source. This solution is representative of high-frequency source positions in real jets and is rich in caustic structures. A full utilization of the ray tool is shown to provide excellent results<br/

A fiber injection unit for the Keck Planet Imager and Characterizer (KPIC)

Coupling a high-contrast imaging instrument to a high-resolution spectrograph has the potential to enable the most detailed characterization of exoplanet atmospheres, including spin measurements and Doppler mapping. The high-contrast imaging system serves as a spatial filter to separate the light from the star and the planet while the high-resolution spectrograph acts as a spectral filter, which differentiates between features in the stellar and planetary spectra. The Keck Planet Imager and Characterizer (KPIC) located downstream from the current W. M. Keck II adaptive optics (AO) system will contain a fiber injection unit (FIU) combining a high-contrast imaging system and a fiber feed to Keck’s high resolution infrared spectrograph NIRSPEC. Resolved thermal emission from known young giant exoplanets will be injected into a single-mode fiber linked to NIRSPEC, thereby allowing the spectral characterization of their atmospheres. Moreover, the resolution of NIRSPEC (R = 37,500) is high enough to enable spin measurements and Doppler imaging of atmospheric weather phenomenon. The module will be integrated and tested at Caltech before being transferred to Keck in 2018

LKB1 and SAD kinases define a pathway required for the polarization of cortical neurons

The polarization of axon and dendrites underlies the ability of neurons to integrate and transmit information in the brain. We show here that the serine/threonine kinase LKB1, previously implicated in the establishment of epithelial polarity and control of cell growth, is required for axon specification during neuronal polarization in the mammalian cerebral cortex. LKB1 polarizing activity requires its association with the pseudokinase Stradalpha and phosphorylation by kinases such as PKA and p90RSK, which transduce neurite outgrowth-promoting cues. Once activated, LKB1 phosphorylates and thereby activates SAD-A and SAD-B kinases, which are also required for neuronal polarization in the cerebral cortex. SAD kinases, in turn, phosphorylate effectors such as microtubule-associated proteins that implement polarization. Thus, we provide evidence in vivo and in vitro for a multikinase pathway that links extracellular signals to the intracellular machinery required for axon specification