Search

We present a look at Ookami, a project providing community access to a testbed supercomputer with the ARM-based A64FX processors developed by a collaboration between RIKEN and Fujitsu and deployed in the Japanese supercomputer Fugaku. We describe the project, provide details about the user base and education/training program, and present highlights from performance studies of two astrophysical simulation codes., Comment: 9 pages, 3 figures, submitted to the Proceedings of 15th International Conference on Numerical Modeling of Space Plasma Flows

Neutron grating interferometry provides information on phase and small-angle scatter in addition to attenuation. Previously, phase grating moir\'e interferometers (PGMI) with two- or three-phase gratings have been developed. These phase-grating systems use the moir\'e far-field technique to avoid the need for high-aspect absorption gratings used in Talbot-Lau interferometers (TLI) which reduce the neutron flux reaching the detector. We first demonstrate through theory and simulations a novel phase grating interferometer system for cold neutrons that requires a single modulated phase grating (MPG) for phase-contrast imaging, as opposed to the two- or three-phase gratings in previously employed PGMI systems. The MPG theory was compared to the full Sommerfeld-Rayleigh Diffraction integral simulator. Then we compared the MPG system to experiments in the literature that use a two-phase-grating-based PGMI with best-case visibility of around 39%. An MPG with a modulation period of 300 micron, pitch of 2 micron, and grating heights with a phase modulation of (pi,0), illuminated by a monochromatic beam, produces a visibility of 94.2% with comparable source-to-detector distance (SDD) as the two-phase-grating-based PGMI. Phase sensitivity, another important performance metric of the grating interferometer, was compared to values available in the literature, viz. the conventional TLI with phase sensitivity of 4.5 x 10E+3 for an SDD of 3.5 m and a beam wavelength of 0.44 nm. For a range of modulation periods, the MPG system provides comparable or greater theoretical maximum phase sensitivity of 4.1 x 10E+3 to 10.0 x 10E+3 for SDDs of up to 3.5 m. This proposed MPG system can provide high-performance PGMI that obviates the need to align two phase gratings., Comment: Manuscript accepted in Rev. Sci. Instrum. vol. 94, 045110, (2023), (Published Online: 17 April 2023)

A high sensitivity Cardiac SPECT system using curved crystals with pinhole collimation was proposed previously (Dey, IEEE TNS 2012, Bhusal et al, Med. Phys. 2019). Here, we hypothesize that a high curvature hemi-ellipsoid detector results in measurable differences in light distribution from events at different depths in the crystal. This was tested by analyzing the scintillation light in hemi-ellipsoid detector using Monte-Carlo (Geant4) and evaluation of both the localization error at detector and the back-projected errors in object-space. We used Geant4 to simulate the propagation of scintillation light in a monolithic hemi-ellipsoidal CsI crystal. A look-up table (LUT) was created to map the points inside the crystal to the expected light pattern on the crystal surface using Geant4. In thirteen zones across the crystal, gamma-rays were simulated and the resulting scintillator light intensity on the surface was captured, serving as our experimental interactions. A Poisson-statistics-based algorithm was developed to limit the search of the gamma-ray event locations into small regions of the LUT and fine-tuned by interpolating between selected LUT points by comparing the light distribution of the gamma interactions and LUT light patterns. The localized events were individually back-projected to the object mid-plane, and the errors recorded. Excluding some outliers (up to 2%), the localized errors averaged over all the zones was 0.71 (+/-0.44) mm with a worse case of 1.36 (+/-0.67) mm at the apex. When back-projected to the midplane of the object for Cardiac SPECT, the errors were <1mm and average error was 0.4(+/-0.22) mm, due to the high system magnification afforded by the detector. Thus, for our high sensitivity system, we were also able to achieve high resolution, assuming perfect pinhole collimator resolution recovery for Cardiac SPECT application., Comment: This work is Accepted (on April 2022) for publication by Journal of Instrumentation; Draft manuscript will be updated with the published version; the Abstract is shortened for arXiv, see full document for complete abstract

We report bulk Rashba spin splitting (RSS) and associated Dirac surface state in (Bi$_{0.9}$Sb$_{0.1})_2$Se$_3$, exhibiting dominant $p$-type conductivity. We argue from the synchrotron diffraction studies that origin of the bulk RSS is due to a structural transition to a non-centrosymmetric $R3m$ phase below $\sim$ 30 K. The Shubnikov-de Haas Van (SdH) oscillations observed in the magnetoresistance curves at low temperature and the Landau level fan diagram, as obtained from these oscillations, confirm the presence of nontrivial Dirac surface state. The magnetization data at low temperature exhibit substantial orbital magnetization consistent with the bulk RSS. The existance of both the bulk RSS and Dirac surface states are confirmed by first principles density functional theory calculations. Coexistence of orbital magnetism, bulk RSS, and Dirac surface state is unique for $p$-type (Bi$_{0.9}$Sb$_{0.1})_2$Se$_3$, making it an ideal candidate for spintronic applications.

The completion of the PIP-II project and its superconducting linear accelerator will provide up to 1.2 MW of beam power to the LBNF/DUNE facility for neutrino physics. It will also be able to produce high-power beams directly from the linac that can be used for lower-energy particle physics experiments as well, such as directing beam toward the Muon Campus at Fermilab for example. Any further significant upgrade of the beam power to DUNE, however, will be impeded by the limitations of the present Booster synchrotron at the facility. To increase the power to DUNE by a factor of two would require a new accelerator arrangement to feed the Main Injector that does not include the Booster. In what follows, a path toward upgrading the Fermilab accelerator complex to bring the beam power for DUNE to 2.4 MW is presented, using a new rapid-cycling synchrotron plus an energy upgrade to the PIP-II linac. The path includes the ability to instigate a new lower-energy, very high-power beam delivery system for experiments that can address much of the science program presented by the Booster Replacement Science Working Group. It also allows for the future possibility to go beyond 2.4 MW up to roughly 4 MW from the Main Injector., Comment: 54 pages, 26 figures (also one table with a figure, counting each figure with multiple images), 35 references. Snowmass 21 White Paper

We report a ferroelectric order at ~ 98 K for NiFe2O4, which carries an inverse spinel structure with a centrosymmetric Fd3m structure at room temperature. The value of spontaneous electric polarization is considerably high as ~ 0.29 {\mu}C/cm2 for 5 kV/cm poling field. The electric polarization decreases considerably (~ 17 %) around liquid nitrogen temperature upon application of 50 kOe field, proposing a significant magnetoelectric coupling. The synchrotron diffraction studies confirm a structural transition at ~ 98 K to a noncentrosymmetric structure of P4122 space group. The occurrence of polar order is associated with an ordered occupancy of Ni and Fe atoms at the octahedral sites of the P4122 structure, instead of random occupancies at the octahedral site of the inverse spinel structure. The results propose that NiFe2O4 is a new type-II multiferroic material., Comment: 8 pages, 7 figures

A second harmonic RF cavity which uses perpendicularly biased garnet for frequency tuning is currently being constructed for use in the Fermilab Booster. The cavity will operate at twice the fundamental RF frequency, from ~76 - 106 MHz, and will be turned on only during injection, and transition or extraction. Its main purpose is to reduce beam loss as required by Fermilab's Proton Improvement Plan (PIP). After three years of optimization and study, the cavity design has been finalized and all constituent parts have been received. We discuss the design aspects of the cavity and its associated systems, component testing, and status of the cavity construction., Comment: 6 pp

A perpendicularly biased tuneable 2nd harmonic cavity is being constructed for use in the Fermilab Booster. The cavity's tuner uses National Magnetics AL800 garnet as the tuning media. For quality control, the magnetic properties of the material and the uniformity of the properties within the tuner must be assessed. We describe two tests which are performed on the rings and on their corresponding witness samples., Comment: 3 pp

A perpendicularly biased 2nd harmonic cavity is being designed and built for the Fermilab Booster. Its purpose is to flatten the bucket at injection and thus change the longitudinal beam distribution to decrease space charge effects. It can also help at extraction. The cavity frequency range is 76 - 106 MHz. The power amplifier will be built using the Y567B tetrode, which is also used for the fundamental mode cavities in the Fermilab Booster. We discuss recent progress on the cavity, the biasing solenoid design and plans for testing the tuner's garnet material., Comment: 3 pp

This paper reports the first measurement using the NOvA detectors of $\nu_\mu$ disappearance in a $\nu_\mu$ beam. The analysis uses a 14 kton-equivalent exposure of $2.74 \times 10^{20}$ protons-on-target from the Fermilab NuMI beam. Assuming the normal neutrino mass hierarchy, we measure $\Delta m^{2}_{32}=(2.52^{+0.20}_{-0.18})\times 10^{-3}$ eV$^{2}$ and $\sin^2\theta_{23}$ in the range 0.38-0.65, both at the 68% confidence level, with two statistically-degenerate best fit points at $\sin^2\theta_{23} = $ 0.43 and 0.60. Results for the inverted mass hierarchy are also presented., Comment: 8 pages, 6 figures. Submitted to Phys. Rev. D Rapid Communications

We report results from the first search for $\nu_\mu\to\nu_e$ transitions by the NOvA experiment. In an exposure equivalent to $2.74\times10^{20}$ protons-on-target in the upgraded NuMI beam at Fermilab, we observe 6 events in the Far Detector, compared to a background expectation of $0.99\pm0.11$ (syst.) events based on the Near Detector measurement. A secondary analysis observes 11 events with a background of $1.07\pm0.14$ (syst.). The $3.3\sigma$ excess of events observed in the primary analysis disfavors $0.1\pi < \delta_{CP} < 0.5\pi$ in the inverted mass hierarchy at the 90% C.L., Comment: 7 pages, 4 figures. Minor updates to match version accepted by journal

The Mu2e experiment at Fermilab will search for charged lepton flavor violation via the coherent conversion process mu- N --> e- N with a sensitivity approximately four orders of magnitude better than the current world's best limits for this process. The experiment's sensitivity offers discovery potential over a wide array of new physics models and probes mass scales well beyond the reach of the LHC. We describe herein the preliminary design of the proposed Mu2e experiment. This document was created in partial fulfillment of the requirements necessary to obtain DOE CD-2 approval., Comment: compressed file, 888 pages, 621 figures, 126 tables; full resolution available at http://mu2e.fnal.gov; corrected typo in background summary, Table 3.4

Mu2e at Fermilab will search for charged lepton flavor violation via the coherent conversion process mu- N --> e- N with a sensitivity approximately four orders of magnitude better than the current world's best limits for this process. The experiment's sensitivity offers discovery potential over a wide array of new physics models and probes mass scales well beyond the reach of the LHC. We describe herein the conceptual design of the proposed Mu2e experiment. This document was created in partial fulfillment of the requirements necessary to obtain DOE CD-1 approval, which was granted July 11, 2012., Comment: 562 pages, 339 figures

For Project X, the Fermilab Main Injector will be required to operate with 3 times higher bunch intensity. The plan to study the space charge effects at the injection energy with intense bunches will be discussed., Comment: 3 pp. Particle Accelerator, 24th Conference (PAC'11) 2011. 28 Mar - 1 Apr 2011. New York, USA

For Project X Fermilab Main Injector will be required to provide up to 2.3 MW to a neutrino production target at energies between 60 and 120 GeV. To accomplish the above power levels 3 times the current beam intensity will need to be accelerated. In addition the injection energy of Main Injector will need to be as low as 6 GeV. The current 30 year old Main Injector radio frequency system will not be able to provide the required power and a new system will be required. The specifications of the new system will be described., Comment: 3 pp. Particle Accelerator, 24th Conference (PAC'11) 2011. 28 Mar - 1 Apr 2011. New York, USA

Linear stability and the non-modal transient energy growth in compressible plane Couette flow are investigated for two prototype mean flows: (a) the {\it uniform shear} flow with constant viscosity, and (b) the {\it non-uniform shear} flow with {\it stratified} viscosity. Both mean flows are linearly unstable for a range of supersonic Mach numbers ($M$). For a given $M$, the critical Reynolds number ($Re$) is significantly smaller for the uniform shear flow than its non-uniform shear counterpart. An analysis of perturbation energy reveals that the instability is primarily caused by an excess transfer of energy from mean-flow to perturbations. It is shown that the energy-transfer from mean-flow occurs close to the moving top-wall for ``mode I'' instability, whereas it occurs in the bulk of the flow domain for ``mode II''. For the non-modal analysis, it is shown that the maximum amplification of perturbation energy, $G_{\max}$, is significantly larger for the uniform shear case compared to its non-uniform counterpart. For $\alpha=0$, the linear stability operator can be partitioned into ${\cal L}\sim \bar{\cal L} + Re^2{\cal L}_p$, and the $Re$-dependent operator ${\cal L}_p$ is shown to have a negligibly small contribution to perturbation energy which is responsible for the validity of the well-known quadratic-scaling law in uniform shear flow: $G(t/{\it Re}) \sim {\it Re}^2$. A reduced inviscid model has been shown to capture all salient features of transient energy growth of full viscous problem. For both modal and non-modal instability, it is shown that the {\it viscosity-stratification} of the underlying mean flow would lead to a delayed transition in compressible Couette flow.

Nonmodal transient growth studies and estimation of optimal perturbations have been made for the compressible plane Couette flow with three-dimensional disturbances. The maximum amplification of perturbation energy over time, $G_{\max}$, is found to increase with increasing Reynolds number ${\it Re}$, but decreases with increasing Mach number $M$. More specifically, the optimal energy amplification $G_{\rm opt}$ (the supremum of $G_{\max}$ over both the streamwise and spanwise wavenumbers) is maximum in the incompressible limit and decreases monotonically as $M$ increases. The corresponding optimal streamwise wavenumber, $\alpha_{\rm opt}$, is non-zero at M=0, increases with increasing $M$, reaching a maximum for some value of $M$ and then decreases, eventually becoming zero at high Mach numbers. While the pure streamwise vortices are the optimal patterns at high Mach numbers, the modulated streamwise vortices are the optimal patterns for low-to-moderate values of the Mach number. Unlike in incompressible shear flows, the streamwise-independent modes in the present flow do not follow the scaling law $G(t/{\it Re}) \sim {\it Re}^2$, the reasons for which are shown to be tied to the dominance of some terms in the linear stability operator. Based on a detailed nonmodal energy analysis, we show that the transient energy growth occurs due to the transfer of energy from the mean flow to perturbations via an inviscid {\it algebraic} instability. The decrease of transient growth with increasing Mach number is also shown to be tied to the decrease in the energy transferred from the mean flow ($\dot{\mathcal E}_1$) in the same limit.

We explore the scenario where the core of a neutron star (having experienced a transition to an up and down quark phase) shrinks into the equilibrated quark object after reaching strange quark matter saturation density (where a composition of up, down and strange quarks is the favored state of matter). The overlaying (envelope) material free-falls following the core contraction releasing upto 10^{53} {\rm ergs} in energy as radiation, partly as a result of the conversion of envelope material to quarks. This phenomena, we named Quark-Nova, leads to a wide variety of ejectae ranging form the Newtonian, "dirty" to the ultra-relativistic fireball. The mass range of the corresponding compact remnant (the quark star) ranges from less than 0.3M_{\odot} up to a solar mass. We discuss the connection between Quark-Novae and Gamma ray bursts and suggest the recently studied GRB011211 event as a plausible Quark-Nova candidate., Comment: 4 pages, 2 figures, shortened title. Major changes in the text, version to appear in Astronomy&Astrophysics

We compute numerical models of uniformly rotating strange stars (SS) in general relativity for the recently proposed QCD-based equation of state (EOS) of strange quark matter (Dey et al. 1998). Static models based on this EOS are characterised by a larger surface redshift than strange stars within the MIT bag model. The frequencies of the fastest rotating configurations described by Dey model are much higher than these for neutron stars (NS) and for the simplest SS MIT bag model. We determine a number of physical parameters for such stars and compare them with those obtained for NS. We construct constant baryon mass equilibrium sequences both normal and supramassive. Similarly to the NS a supramassive SS, prior to collapse to a black hole, spins up as it loses angular momentum. We find the upper limits on maximal masses and maximal frequencies of the rotating configurations. We show that the maximal rotating frequency for each of considered evolutionary sequences is never the Keplerian one. A normal and low mass supramassive strange stars gaining angular momentum always slows down just before reaching the Keplerian limit. For a high mass supramassive SS sequence the Keplerian configuration is the one with the lowest rotational frequency in the sequence. The value of $T/W$ for rapidly rotating SS of any mass is significantly higher than those for ordinary NS. For Keplerian configurations it increases as mass decreases. The results are robust for all linear self-bound equations of state., Comment: 8 pages, 5 eps figures, submitted to A&A

Compactness of some stars is explained if they are strange stars (SS) as shown by Dey et al. (1998) (D98) and Li et al. (1999a). One of these compact star candidates is the SAX J1808.4-3658 (SAX in short) believed to be an important link in the genesis of radio pulsars. SS have also been suggested for bursting X-ray pulsars (GRO J1744-28, Cheng et al. 1998), from quasi-periodic oscillations (QPO) of X-ray binaries (4U 1728-34, Li et al. 1999b) and from peculiarity of properties of radio pulsars (PSR 0943+10, Xu et al.1999; Kapoor et al. 2000). We now extend the calculation to include high temperatures upto T = 70 MeV ~ 8 times 10^{11} ^oK and find that the nature of the mass (M) and radius (R), derived from astrophysical data, is still retained. The entropy is calculated and matches onto that calculated from hadronic models thus supporting the idea that the quark-hadron transition may be continuous., Comment: 4 pages, including 4 figures