We present Powderday, a flexible, fast, open-source dust radiative transfer package designed to interface with galaxy formation simulations. Powderday builds on FSPS population synthesis models, Hyperion dust radiative transfer, and employs yt to interface between different software packages. We include our stellar population synthesis modeling on the fly, which allows for significant run-time flexibility in the assumed stellar physics. We include a model for nebular line emission that can employ either precomputed Cloudy lookup tables (for efficiency), or direct photoionization calculations for all young stars (for flexibility). The dust content follows either observationally-motivated prescriptions, direct modeling from galaxy formation simulations, or a novel approach that includes the dust content via learning-based algorithms from the SIMBA cosmological galaxy formation simulation. AGN can additionally be included via a range of prescriptions. The output of these models are broadband SEDs, as well as filter-convolved images. Powderday is designed to eliminate last-mile efforts by researchers that employ different hydrodynamic galaxy formation models, and seamlessly interfaces with GIZMO, AREPO, GASOLINE, CHANGA, and ENZO. We demonstrate the capabilities of the code via three applications: a model for the star formation rate (SFR) - infrared luminosity relation in galaxies (including the impact of AGN); the impact of circumstellar dust around AGB stars on the mid-infrared emission from galaxy SEDs; and the impact of galaxy inclination angle on dust attenuation laws.

This work presents the broad-band time-averaged spectral analysis of neutron star low-mass X-ray binary, XTE J1710-281 by using the Suzaku archival data. The source was in a hard or an intermediate spectral state during this observation. This is the first time that a detailed spectral analysis of the persistent emission spectra of XTE J1710-281 has been done up to 30 keV with improved constraints on its spectral parameters. By simultaneously fitting the XIS (0.6-9.0 keV) and the HXD-PIN (15.0-30.0 keV) data, we have modelled the persistent spectrum of the source with models comprising a soft component from accretion disc and/or neutron star surface/boundary layer and a hard Comptonizing component. The 0.6-30 keV continuum with neutral absorber can be described by a multi-colour disc blackbody with an inner disc temperature of $kT_{\rm disc} = 0.28$ keV, which is significantly Comptonized by the hot electron cloud with electron temperature of $kT_{\rm e} \approx 5$ keV and described by photon index $\Gamma = 1.86$. A more complex three-component model comprising a multi-colour disc blackbody $\approx 0.30$ keV, single temperature blackbody $\approx 0.65$ keV and Comptonization from the disc, partially absorbed (about 38 per cent) by an ionized absorber (log($\xi$) $\approx$ 4) describes the broad-band spectrum equally well.

This paper presents EinsteinPy (version 0.3), a community-developed Python package for gravitational and relativistic astrophysics. Python is a free, easy to use a high-level programming language which has seen a huge expansion in the number of its users and developers in recent years. Specifically, a lot of recent studies show that the use of Python in Astrophysics and general physics has increased exponentially. We aim to provide a very high level of abstraction, an easy to use interface and pleasing user experience. EinsteinPy is developed keeping in mind the state of a theoretical gravitational physicist with little or no background in computer programming and trying to work in the field of numerical relativity or trying to use simulations in their research. Currently, EinsteinPy supports simulation of time-like and null geodesics and calculates trajectories in different background geometries some of which are Schwarzschild, Kerr, and KerrNewmann along with coordinate inter-conversion pipeline. It has a partially developed pipeline for plotting and visualization with dependencies on libraries like Plotly, matplotlib, etc. One of the unique features of EinsteinPy is a sufficiently developed symbolic tensor manipulation utilities which are a great tool in itself for teaching yourself tensor algebra which for many beginner students can be overwhelmingly tricky. EinsteinPy also provides few utility functions for hypersurface embedding of Schwarzschild spacetime which further will be extended to model gravitational lensing simulation.

We use the Romulus25 cosmological simulation volume to identify the largest-ever simulated sample of field ultra-diffuse galaxies (UDGs). At z=0, we find that isolated UDGs have average star formation rates, colors, and virial masses for their stellar masses and environment. Although they tend to have slightly elevated HI masses, UDGs are consistent with the general isolated dwarf galaxy population and make up ~20% of all field galaxies with 10$^7$<M$_\star$/M$_\odot$<10$^{9}$. The properties of our UDGs are consistent with existing observations of field UDGs, but we predict that many isolated UDGs have been missed by current surveys. Despite their isolation at z=0, the UDGs in our sample are the products of major mergers at z>1. Mergers are no more common in UDG than non-UDG progenitors, but mergers that create UDGs tend to happen earlier, produce a temporary boost in spin, and cause star formation to migrate to the outskirts of galaxies, resulting in lower central star formation rates. The centers of the galaxies fade as their central stellar populations age, but their global star formation rates are maintained through bursts of star formation at larger radii. This formation channel is unique relative to other proposals for UDG formation in isolated galaxies, or for UDGs in groups and clusters, demonstrating that UDGs can potentially be formed through multiple mechanisms.

We reconstruct late-time cosmology in a model-independent manner using the technique of Principal Component Analysis (PCA). In particular, we focus on the reconstruction of the dark energy equation of state parameters from two different observational data sets, supernova type Ia data, and the Hubble parameter data. To achieve this reconstruction, we have adopted two different techniques. The first is a derived approach wherein we reconstruct the observable quantities of the data sets, namely the Hubble parameter and the supernova distance modulus from observations using PCA and subsequently reconstruct the allowed equation of state parameter. The other approach is a direct one where dark energy equation of state is reconstructed directly from the data sets. We show that a combination of PCA algorithm and calculation of correlation coefficients can be used as a tool of reconstruction. The derived approach is found to be statistically preferable over the direct approach. We have carried out the analysis with simulated data and observed data sets of Hubble parameter measurements and distance modulus measurements of type Ia supernova. The reconstructed equation of state indicates a slowly varying dark energy equation of state parameter.

We report the discovery of an intermediate-mass transiting brown dwarf, TOI-503b, from the TESS mission. TOI-503b is the first brown dwarf discovered by TESS and orbits a metallic-line A-type star with a period of $P=3.6772 \pm 0.0001$ days. The light curve from TESS indicates that TOI-503b transits its host star in a grazing manner, which limits the precision with which we measure the brown dwarf's radius ($R_b = 1.34^{+0.26}_{-0.15} R_J$). We obtained high-resolution spectroscopic observations with the FIES, Ond\v{r}ejov, PARAS, Tautenburg, and TRES spectrographs and measured the mass of TOI-503b to be $M_b = 53.7 \pm 1.2 M_J$. The host star has a mass of $M_\star = 1.80 \pm 0.06 M_\odot$, a radius of $R_\star = 1.70 \pm 0.05 R_\odot$, an effective temperature of $T_{\rm eff} = 7650 \pm 160$K, and a relatively high metallicity of $0.61\pm 0.07$ dex. We used stellar isochrones to derive the age of the system to be $\sim$180 Myr, which places its age between that of RIK 72b (a $\sim$10 Myr old brown dwarf in the Upper Scorpius stellar association) and AD 3116b (a $\sim$600 Myr old brown dwarf in the Praesepe cluster). We argue that this brown dwarf formed in-situ, based on the young age of the system and the long circularization timescale for this brown dwarf around its host star. TOI-503b joins a growing number of known short-period, intermediate-mass brown dwarfs orbiting main sequence stars, and is the second such brown dwarf known to transit an A star, after HATS-70b. With the growth in the population in this regime, the driest region in the brown dwarf desert ($35-55 M_J \sin{i}$) is reforesting and its mass range shrinking.

We have conducted a survey of 328 protostars in the Orion molecular clouds with ALMA at 0.87 mm at a resolution of $\sim$0.1" (40 au), including observations with the VLA at 9 mm toward 148 protostars at a resolution of $\sim$0.08" (32 au). This is the largest multi-wavelength survey of protostars at this resolution by an order of magnitude. We use the dust continuum emission at 0.87 mm and 9 mm to measure the dust disk radii and masses toward the Class 0, Class I, and Flat Spectrum protostars, characterizing the evolution of these disk properties in the protostellar phase. The mean dust disk radii for the Class 0, Class I, and Flat Spectrum protostars are 44.9$^{+5.8}_{-3.4}$, 37.0$^{+4.9}_{-3.0}$, and 28.5$^{+3.7}_{-2.3}$ au, respectively, and the mean protostellar dust disk masses are 25.9$^{+7.7}_{-4.0}$, 14.9$^{+3.8}_{-2.2}$, 11.6$^{+3.5}_{-1.9}$ Earth masses, respectively. The decrease in dust disk masses is expected from disk evolution and accretion, but the decrease in disk radii may point to the initial conditions of star formation not leading to the systematic growth of disk radii or that radial drift is keeping the dust disk sizes small. At least 146 protostellar disks (35% out of 379 detected 0.87 mm continuum sources plus 42 non-detections) have disk radii greater than 50 au in our sample. These properties are not found to vary significantly between different regions within Orion. The protostellar dust disk mass distributions are systematically larger than that of Class II disks by a factor of $>$4, providing evidence that the cores of giant planets may need to at least begin their formation during the protostellar phase.

SAX J1748.9-2021 is a transient accretion powered millisecond X-ray pulsar located in the Globular cluster NGC 6440. We report on the spectral and timing analysis of SAX J1748.9-2021 performed on AstroSat data taken during its faint and short outburst of 2017. We derived the best-fitting orbital solution for the 2017 outburst and obtained an average local spin frequency of 442.361098(3) Hz. The pulse profile obtained from 3-7 keV and 7-20 keV energy bands suggest constant fractional amplitude ~0.5% for fundamental component, contrary to previously observed energy pulse profile dependence. Our AstroSat observations revealed the source to be in a hard spectral state. The 1-50 keV spectrum from SXT and LAXPC on-board AstroSat can be well described with a single temperature blackbody and thermal Comptonization. Moreover, we found that the combined spectra from XMM-Newton (EPIC-PN) and AstroSat (SXT+LAXPC) indicated the presence of reflection features in the form of iron (Fe K${\alpha}$) line that we modeled with the reflection model xillvercp. One of the two X-ray burst observed during the AstroSat/LAXPC observation showed hard X-ray emission (>30 keV) due to Compton up-scattering of thermal photons by the hot corona. Time resolved analysis performed on the bursts revealed complex evolution in emission radius of blackbody for second burst suggestive of mild photospheric radius expansion.

Generation of magnetic fields during inflation is a promising mechanism for the origin of the observed large scale magnetic fields in the universe. Among several attempts, a popular model is one where the inflaton and the electromagnetic field are coupled through a coupling function $f$ leading to a term in the Lagrangian density of the form, $f^2 F^{\mu \nu} F_{\mu\nu}$. A number of potential difficulties with such models have been raised in the literature. In our earlier work, we have suggested viable models of inflationary magnetogenesis which avoid these problems and at the same time can lead to either nonhelical or helical magnetic fields of astrophysical interest. Our models require a low energy scale for inflation and reheating (reheating temperature, $T_R < 10^4$ GeV) and generate a blue spectrum of electromagnetic (EM) field which peaks around the horizon scale of reheating. We show here that the anisotropic stress associated with these EM fields naturally source the production of a stochastic background of Gravitational waves (GW) with frequencies in the range of tens of nano Hertz to milli Hertz. These two extremes of the range can be probed respectively by pulsar timing arrays (PTA) experiments and the upcoming Laser Interferometric Space Array (LISA). The peak value of the GW spectrum energy represented by $d \Omega_{GW}/d \ln k$ is $ 10^{-6} $ for the models which lead to nonhelical primordial fields and $2 \times 10^{-6} $ for the helical case for $T_R=100$ GeV. In this case the spectrum peaks at a frequency $30 \mu$Hz for non helical case and at $40 \mu$Hz for helical case. These values are obtained when the ratio of EM energy density to the cosmological density at reheating $\epsilon \sim 1$ and decrease approximately as $\epsilon^2$ for smaller values.

We investigate the effects of massive black hole growth on the structural evolution of dwarf galaxies within the Romulus25 cosmological hydrodynamical simulation. We study a sample of 228 central, isolated dwarf galaxies with stellar masses $M_{star} < 10^{10} M_\odot$ and a central BH. We find that the local $M_{BH} - M_{star}$ relation exhibits a high degree of scatter below $M_{star} < 10^{10} M_\odot$, which we use to classify BHs as overmassive or undermassive relative to their host $M_{star}$. Overmassive BHs grow through a mixture of BH mergers and relatively high average accretion rates, while undermassive BHs grow slowly through accretion. We find that isolated dwarf galaxies that host overmassive BHs also follow different evolutionary tracks relative to their undermassive BH counterparts, building up their stars and dark matter earlier and experiencing star formation suppression starting around $z=2$. By $z=0.05$, overmassive BH hosts above $M_{star} > 10^{9} M_\odot$ are more likely to exhibit lower central stellar mass density, lower HI gas content, and lower star formation rates than their undermassive BH counterparts. Our results suggest that overmassive BHs in isolated galaxies above $M_{star} > 10^{9} M_\odot$ are capable of driving feedback, in many cases suppressing and even quenching star formation by late times.

We present ALMA (0.87~mm) and VLA (9~mm) observations toward OMC2-FIR4 and OMC2-FIR3 within the Orion integral-shaped filament that are thought to be the nearest regions of intermediate mass star formation. We characterize the continuum sources within these regions on $\sim$40~AU (0\farcs1) scales and associated molecular line emission at a factor of $\sim$30 better resolution than previous observations at similar wavelengths. We identify six compact continuum sources within OMC2-FIR4, four in OMC2-FIR3, and one additional source just outside OMC2-FIR4. This continuum emission is tracing the inner envelope and/or disk emission on less than 100~AU scales. HOPS-108 is the only protostar in OMC2-FIR4 that exhibits emission from high-excitation transitions of complex organic molecules (e.g., methanol and other lines) coincident with the continuum emission. HOPS-370 in OMC2-FIR3 with L~$\sim$~360~\lsun, also exhibits emission from high-excitation methanol and other lines. The methanol emission toward these two protostars is indicative of temperatures high enough to thermally evaporate methanol from icy dust grains; overall these protostars have characteristics similar to hot corinos. We do not identify a clear outflow from HOPS-108 in \twco, but find evidence of interaction between the outflow/jet from HOPS-370 and the OMC2-FIR4 region. The multitude of observational constraints indicate that HOPS-108 is likely a low to intermediate-mass protostar in its main mass accretion phase and it is the most luminous protostar in OMC2-FIR4. The high resolution data presented here are essential for disentangling the embedded protostars from their surrounding dusty environments and characterizing them.

We present a first proof-of-principle study for using deep neural networks (DNNs) as a novel search method for continuous gravitational waves (CWs) from unknown spinning neutron stars. The sensitivity of current wide-parameter-space CW searches is limited by the available computing power, which makes neural networks an interesting alternative to investigate, as they are extremely fast once trained and have recently been shown to rival the sensitivity of matched filtering for black-hole merger signals. We train a convolutional neural network with residual (short-cut) connections and compare its detection power to that of a fully-coherent matched-filtering search using the WEAVE pipeline. As test benchmarks we consider two types of all-sky searches over the frequency range from $20\,\mathrm{Hz}$ to $1000\,\mathrm{Hz}$: an `easy' search using $T=10^5\,\mathrm{s}$ of data, and a `harder' search using $T=10^6\,\mathrm{s}$. Detection probability $p_\mathrm{det}$ is measured on a signal population for which matched filtering achieves $p_\mathrm{det}=90\%$ in Gaussian noise. In the easiest test case ($T=10^5\,\mathrm{s}$ at $20\,\mathrm{Hz}$) the DNN achieves $p_\mathrm{det}\sim88\%$, corresponding to a loss in sensitivity depth of $\sim5\%$ versus coherent matched filtering. However, at higher-frequencies and longer observation time the DNN detection power decreases, until $p_\mathrm{det}\sim13\%$ and a loss of $\sim 66\%$ in sensitivity depth in the hardest case ($T=10^6\,\mathrm{s}$ at $1000\,\mathrm{Hz}$). We study the DNN generalization ability by testing on signals of different frequencies, spindowns and signal strengths than they were trained on. We observe excellent generalization: only five networks, each trained at a different frequency, would be able to cover the whole frequency range of the search.

In the present paper, we investigate the power-law behaviour of the magnetic field spectra in the Earths magnetosheath region using Cluster spacecraft data under solar minimum condition. The power spectral density of the magnetic field data and spectral slopes at various frequencies are analysed. Propagation angle and compressibility are used to test the nature of turbulent fluctuations. The magnetic field spectra have the spectral slopes between -1.5 to 0 down to spatial scales of 20 ion gyroradius and show clear evidence of a transition to steeper spectra for small scales with a second power-law, having slopes between -2.6 to -1.8. At low frequencies, f_sc<0.3f_ci(where f_ci is ion gyro-frequency), propagation angle approximately 90 degrees to the mean magnetic field, B_0, and compressibility shows a broad distribution, 0.1 < R > 0.9. On the other hand at f_sc>10f_ci, the propagation angle exhibits a broad range between 30-90 degree while 'R' has a small variation: 0.2 < R > 0.5. We conjecture that at high frequencies, the perpendicularly propagating Alfven waves could partly explain the statistical analysis of spectra. To support our prediction of kinetic Alfven wave-dominated spectral slope behaviour at high frequency, we also present a theoretical model and simulate the magnetic field turbulence spectra due to the nonlinear evolution of kinetic Alfven waves. The present study also shows the analogy between the observational and simulated spectra.

With 30 antennas and a maximum baseline length of 25 km, the Giant Metrewave Radio Telescope (GMRT) is the premier low-frequency radio interferometer today. We have carried out a study of possible expansions of the GMRT, via adding new antennas and installing focal plane arrays (FPAs), to improve its point-source sensitivity, surface brightness sensitivity, angular resolution, field of view, and U-V coverage. We have carried out array configuration studies, aimed at minimizing the number of new GMRT antennas required to obtain a well-behaved synthesized beam over a wide range of angular resolutions for full-synthesis observations. This was done via two approaches, tomographic projection and random sampling, to identify the optimal locations for the new antennas. We report results for the optimal locations of the antennas of an expanded array (the "EGMRT"), consisting of the existing 30 GMRT antennas, 30 new antennas at short distances, $\leq 2.5$ km from the array centre, and 26 new antennas at long distances, $\approx 5-25$ km from the array centre. The collecting area and the field of view of the proposed EGMRT array would be larger by factors of, respectively, $\approx 3$ and $\approx 30$, than those of the GMRT. Indeed, the EGMRT continuum sensitivity and survey speed with 550-850 MHz FPAs installed on the 45 antennas within a distance of $\approx 2.5$ km of the array centre would be far better than those of any existing interferometer, and comparable to the sensitivity and survey speed of Phase-1 of the Square Kilometre Array.

Context. Being the most numerous and oldest stars in the galaxy, M dwarfs are objects of great interest for exoplanet searches. The presence of molecules in their atmosphere complicates our understanding of their atmospheric properties. But great advances have recently been made in the modeling of M dwarfs due to the revision of solar abundances. Aims. We aim to determine stellar parameters of M dwarfs using high resolution spectra (R = 90 000) simultaneously in the visible and the near-infrared. The high resolution spectra and broad wavelength coverage provide an unique opportunity to understand the onset of dust and cloud formation at cool temperatures. Furthermore, this study will help in understanding the physical processes which occur in a cool atmospheres, particularly, the redistribution of energy from the optical to the near-infrared. Methods. The stellar parameters of M dwarfs in our sample have been determined by comparing the high resolution spectra both in the optical and in the near-infrared simultaneously observed by CARMENES with the synthetic spectra obtained from the BT-Settl model atmosphere. The detailed spectral synthesis of these observed spectra both in the optical and in the near-infrared helps to understand the missing continuum opacity. Results. For the first time, we derive fundamental stellar parameters of M dwarfs using the high resolution optical and near-infrared spectra simultaneously. We determine Teff , log g and [M/H] for 292 M dwarfs of spectral type M0 to M9, where the formation of dust and clouds are important. The derived Teff for the sample ranges from 2300 to 4000 K, values of log g ranges from 4.5 geq log g leq 5.5 and the resulting metallicity ranges from -0.5 geq [M/H] leq +0.5. We have also explored the possible differences in Teff , log g and [M/H] by comparing them with other studies of the same sample of M dwarfs.

We report the spectral analysis of accretion powered millisecond X-ray pulsar SAX J1748.9-2021 from the NuSTAR observation made during its 2015 outburst. The spectra displayed a broad emission line at $\sim 6.5$ keV with line width of $\sim 0.5$ keV and excess above $\sim 20$ keV due to Compton hump. The persistent emission of SAX J1748.9-2021 is described by a combination of soft thermal component with $kT =0.62^{+0.03}_{-0.05}$ keV and thermally Comptonized component with $kT_e=2.50^{+0.06}_{-0.03}$ keV reflected from the disc with reflection fraction of $0.30^{+0.08}_{-0.04}$. We have modeled the reflection spectrum with self-consistent model relxillCP and have found the inclination of the accretion disc to be $32.3^{^\circ +4.8}_{-4.7}$ and log $\xi$ is equal to $3.05^{+0.40}_{-0.34}$. We have obtained an upper limit of 57 km for inner disc radius; and the estimated magnetic field strength at the poles is less than $3.4 \times 10^9$ G. This kind of detailed investigation of reflection spectrum of SAX J1748.9-2021, especially above 10 keV, has been achieved only because of high statistics NuSTAR data of the source.

The formation and evolution of disks early in the protostellar phase is an area of study in which the ngVLA is poised to make significant breakthroughs. The high-sensitivity and resolution at wavelengths of 3~mm and longer will enable forming disks to be examined with unprecedented detail. The need to observe dust emission at wavelengths of 3~mm and longer is motivated by the fact that dust emission at these wavelengths is more likely to be optically thin, which is essential to understanding the structure of these disks. We explore the feasibility of detecting and resolving protostellar disks with a variety of radii, masses, and distances, out to distances as large as 1.5~kpc using radiative transfer models and simulations with the proposed ngVLA configuration. We also examine the potential for the ngVLA to enable studies of grain growth and radial migration of dust particles early in the protostellar phase with the broad multi-wavelength coverage. Studies of grain growth will require wavelength coverage extending at least to $\sim$4~cm to characterize and quantify the location and intensity of free-free emission, which is expected to be generated at $<$10 AU scales.

Within the framework of resistive MHD, implementing the C7 equilibrium atmosphere model and a 3D potential magnetic field realistic configuration, we simulate the formation of a plasm jet with the morphology, upward velocity up to 130 km/s and timescale formation between 60 and 90 s after beginning of simulation, similar to those expected for Type II spicules. Initial results of this simulation were published in Paper (e.g., Gonz\'alez-Avil\'es et al. 2018) and present paper is devoted to the analysis of transverse displacements and rotational type motion of the jet. Our results suggest that 3D magnetic reconnection may be responsible for the formation of the jet in Paper (Gonz\'alez-Avil\'es et al. 2018). In this paper, by calculating times series of the velocity components $v_x$ and $v_y$ in different points near to the jet for various heights we find transverse oscillations in agreement with spicule observations. We also obtain a time-distance plot of the temperature in a cross-cut at the plane $x=$0.1 Mm and find significant transverse displacements of the jet. By analyzing temperature isosurfaces of $10^{4}$ K with the distribution of $v_x$, we find that if the line-of-sight (LOS) is approximately perpendicular to the jet axis then there is both motion towards and away from the observer across the width of the jet. This red-blue shift pattern of the jet is caused by rotational motion, initially clockwise and anti-clockwise afterwards, which could be interpreted as torsional motion. From a nearly vertical perspective of the jet the LOS velocity component shows a central blue-shift region surrounded by red-shifted plasma.

We report results from a broadband spectral analysis of the low mass X-ray binary MXB 1658-298 with the Swift/XRT and NuSTAR observations made during its 2015-2017 outburst. The source showed different spectral states and accretion rates during this outburst. The source was in low/hard state during 2015; and it was in high/soft state during the 2016 NuSTAR observations. This is the first time that a comparison of the soft and hard spectral states during an outburst is being reported in MXB 1658-298. Compared with the observation of 2015, the X-ray luminosity was about four times higher during 2016. The hard state spectrum can be well described with models consisting of a single-temperature blackbody component along with Comptonized disc emission, or three component model comprising of multi-colour disc, single-temperature blackbody and thermal Comptonization components. The soft state spectrum can be described with blackbody or disc blackbody and Comptonization component, where the neutron star surface (or boundary layer) is the dominant source for the Comptonization seed photons. We have also found a link between the spectral state of the source and the Fe K absorber. The absorption features due to highly ionized Fe were observed only in the soft state. This suggests a probable connection between accretion disc wind (or atmosphere) and spectral state (or accretion state) of MXB 1658-298.

We present the first detection of OH absorption in diffuse gas at $z>0$, along with another eight stringent limits on OH column densities for cold atomic gas in galaxies at $0<z<0.4$. The absorbing gas detected towards Q0248+430 ($z_q$=1.313) originates from a tidal tail emanating from a highly star forming galaxy G0248+430 ($z_g$=0.0519) at an impact parameter of 15 kpc. The measured column density is $N$(OH) = (6.3$\pm$0.8)$\times$10$^{13}$($\frac{T_{\rm ex}}{3.5}$)($\frac{1.0}{f_c^{\rm OH}}$) cm$^{-2}$, where $f_c^{\rm OH}$ and $T_{\rm ex}$ are the covering factor and the excitation temperature of the absorbing gas, respectively. In our Galaxy, the column densities of OH in diffuse clouds are of the order of $N$(OH)$\sim$10$^{13-14}$ cm$^{-2}$. From the incidence (number per unit redshift; $n_{21}$) of HI 21-cm absorbers at $0.5<z<1$ and assuming no redshift evolution, we estimate the incidence of OH absorbers (with log$N$(OH)$>$13.6) to be $n_{\rm OH}$ = $0.008^{+0.018}_{-0.008}$ at $z\sim0.1$. Based on this we expect to detect 10$^{+20}_{-10}$ such OH absorbers from the MeerKAT Absorption Line Survey. Using HI 21-cm and OH 1667 MHz absorption lines detected towards Q0248+430, we estimate ($\Delta F/F$) = (5.2$\pm$4.5)$\times 10^{-6}$, where $F \equiv g_p (\alpha^2/\mu)^{1.57}$, $\alpha$ $-$ the fine structure constant, $\mu$ $-$ the electron-proton mass ratio and $g_p$ $-$ the proton gyromagnetic ratio. This corresponds to $\Delta\alpha/\alpha$($z=0.0519$) = (1.7 $\pm$ 1.4)$\times 10^{-6}$, which is among the stringent constraints on the fractional variation of $\alpha$.

We study the possibility of reheating the universe in its early stages through the evaporation of Primordial Black Holes (PBHs) that are formed due to the collapse of the inhomogeneities that were generated during inflation. By using the current results of the baryon-photon ratio obtained from BBN and CMB observations, we impose constraints on the spectral index of perturbations on those small scales that cannot be estimated through CMB anisotropy and CMB distortions. The masses of the PBHs constrained in this study lie in the range of $10^{9}$ and $10^{11}$g, which corresponds to those PBHs whose maximal evaporation took place during the redshifts $10^{6} < z < 10^{9}$. It is shown that the upper bound on the scalar spectral index, $ n_{s}$ can be constrained for a given threshold value, $ \zeta _{\rm th}$, of the curvature perturbations for PBHs formation. Using Planck results for cosmological parameters we obtained $n_{s} < 1.309 $ for $ \zeta _{\rm th} =0.7 $ and $ n_{s} < 1.334 $ for $ \zeta _{\rm th} = 1.2 $ respectively. The density fraction that has contributed to the formation of Primordial Black Holes has also been estimated.

Eclipsing Binaries (EBs) with one of the companions as very low mass stars (VLMS or M dwarfs) are testbeds to substantiate stellar models and evolutionary theories. Here, we present four EB candidates with F type primaries, namely, SAO 106989, HD 24465, EPIC 211682657 and HD 205403, identified from different photometry missions, SuperWasp, KELT, Kepler 2 (K2) and STEREO. Using the high-resolution spectrograph, PARAS, at the 1.2 m telescope at Mount Abu, Rajasthan, India, we hereby report the detection of four VLMS as companions to the four EBs. We performed spectroscopic analysis and found the companion masses to be $0.256\pm0.005$, $0.233\pm0.002$, $0.599\pm0.017$, and $0.406\pm0.005 M_{\odot}$ for SAO 106989, HD 24465, EPIC 211682657 and HD 205403 respectively. We determined orbital periods of $4.39790\pm0.00001$, $7.19635\pm0.00002$, $3.142023\pm0.000003$ and $2.444949\pm0.000001$~d and eccentricities of $0.248\pm0.005$, $0.208\pm0.002$, $0.0097\pm0.0008$ and $0.002\pm0.002$ for EBs SAO 106989, HD 24465, EPIC 211682657 and HD 205403 respectively. The radius derived by modeling the photometry data are $0.326\pm0.012~R_{\odot}$ for SAO 106989, $0.244\pm0.001 R_{\odot}$ for HD 24465, $0.566\pm0.005 R_{\odot}$ for EPIC 211682657, and $0.444\pm0.014 R_{\odot}$ for HD 205403. The radii of HD 24465B and EPIC 211682657B have been measured by precise KEPLER photometry and are consistent with theory within error bars. However, the radii of SAO 106989B and HD 205403B, measured by KELT and STEREO photometry, are $17-20\%$ higher than those predicted by theory. A brief comparison of the results of the current work is made with the M dwarfs already studied in literature.

We report here strong evidence for a sub-Saturn around EPIC~211945201 and confirm its planetary nature. EPIC~211945201b was found to be a planetary candidate from {\it K2} photometry in Campaigns 5 \& 16, transiting a bright star ($V_{\rm mag}=10.15$, G0 spectral type) in a 19.492 day orbit. However, the photometric data combined with false positive probability calculations using VESPA was not sufficient to confirm the planetary scenario. Here we present high-resolution spectroscopic follow-up of the target using the PARAS spectrograph (19 radial velocity observations) over a time-baseline of 420 days. We conclusively rule out the possibility of an eclipsing binary system and confirm the 2-$\sigma$ detection of a sub-Saturn planet. The confirmed planet has a radius of 6.12$\pm0.1$$~R_{\oplus}$, and a mass of $27_{-12.6}^{+14}$~$M_{\oplus}$. We also place an upper limit on the mass (within the 3-$\sigma$ confidence interval) at 42~$M_{\oplus}$ above the nominal value. This results in the Saturn-like density of $0.65_{-0.30}^{+0.34}$ g~cm$^{-3}$. Based on the mass and radius, we provide a preliminary model-dependent estimate that the heavy element content is 60-70 \% of the total mass. This detection is important as it adds to a sparse catalog of confirmed exoplanets with masses between 10-70 $M_{\oplus}$ and radii between 4-8 $R_{\oplus}$, whose masses and radii are measured to a precision of 50\% or better (only 23 including this work).

LS I +61$^{\circ}$303 is an X-ray binary with a radio outburst every ~27 days. Previous studies of the stellar system revealed radio microflares superimposed on the large radio outburst. We present here new radio observations of LS I +61$^{\circ}$303 at 2.2 GHz with the Westerbork Synthesis Radio Telescope (WSRT). Using various timing analysis methods we find significant Quasi-Periodic Oscillations (QPOs) of 55 minutes stable over the duration of 4 days. We also use archival data obtained from the Suzaku satellite at X-ray wavelengths. We report here for the first time significant X-ray QPOs of about 2 hours present over the time span of 21 hours. We compare our results with the previously reported QPO observations and we conclude that the QPOs seem to be associated with the radio outburst, independent of the amplitude of the outburst. Finally, the different QPO time-scales are discussed in the context of magnetic reconnection.

We study the generation of helical magnetic fields in a model of inflationary magnetogenesis which is free from the strong coupling and back-reaction problems. To generate helical magnetic fields, we add an $f^2 \tilde{F}^{\mu\nu} F_{\mu\nu}$ term to the lagrangian of Ratra model. The strong coupling and back-reaction problems are avoided if we take a particular behaviour of coupling function $f$, in which $f$ increases during inflation and decreases post inflation to reheating. The generated magnetic field is fully helical and has a blue spectrum, $d\rho_B/d\ln k \propto k^4$. This spectrum is obtained when coupling function $f\propto a^2$ during inflation. The scale of reheating in our model has to be lower than $4000$ GeV to avoid back-reaction post inflation. The generated magnetic field spectrum satisfies the $\gamma$-ray bound for all the possible scales of reheating. The comoving magnetic field strength and its correlation length are $\sim 4 \times 10^{-11} $ G and $70$ kpc respectively, if reheating takes place at 100 GeV. For reheating at the QCD scales of $150$ MeV, the field strength increases to $\sim$ nano gauss, with coherence scale of $0.6$ Mpc.

Models of inflationary magnetogenesis with a coupling to the electromagnetic action of the form $f^2 F_{\mu\nu}F^{\mu\nu}$, are known to suffer from several problems. These include the strong coupling problem, the back reaction problem and also strong constraints due to Schwinger effect. We propose a model which resolves all these issues. In our model, the coupling function, $f$, grows during inflation and transits to a decaying phase post inflation. This evolutionary behaviour is chosen so as to avoid the problem of strong coupling. By assuming a suitable power law form of the coupling function, we can also neglect back reaction effects during inflation. To avoid back reaction post-inflation, we find that the reheating temperature is restricted to be below $ \approx 1.7 \times 10^4$ GeV. The magnetic energy spectrum is predicted to be non-helical and generically blue. The estimated present day magnetic field strength and the corresponding coherence length taking reheating at the QCD epoch(150 MeV) are $ 1.4 \times 10^{-12}$ G and $6.1 \times 10^{-4}$ Mpc, respectively. This is obtained after taking account of nonlinear processing over and above the flux freezing evolution after reheating. If we consider also the possibility of a non-helical inverse cascade, as indicated in direct numerical simulations, the coherence length and the magnetic field strength are even larger. In all cases mentioned above, the magnetic fields generated in our models satisfy the $\gamma$-ray bound below a certain reheating temperature.

We investigate particle production \`a la Schwinger mechanism in an expanding, flat de Sitter patch as is relevant for the inflationary epoch of our universe. Defining states and particle content in curved spacetime is certainly not a unique process. There being different prescriptions on how that can be done, we have used the Schr\"odinger formalism to define instantaneous particle content of the state etc. This allows us to go past the adiabatic regime to which the effect has been restricted in the previous studies and bring out its multifaceted nature in different settings. Each of these settings gives rise to contrasting features and behaviour as per the effect of electric field and expansion rate on the instantaneous mean particle number. We also quantify the degree of classicality of the process during its evolution using a "classicality parameter" constructed out of parameters of the Wigner function to obtain information about the quantum to classical transition in this case.

The recent availability of fine grained high sensitivity data from the new generation low radio frequency instruments such as the Murchison Widefield Array (MWA) have opened up opportunities for using novel techniques for characterizing the nature of solar emission at these frequencies. Here we use this opportunity to look for evidence for the presence of weak non-thermal emissions in the 100-240 MHz band, at levels weaker than have usually been probed. The presence of such features is believed to be a necessary consequence of nanoflare-based coronal and chromospheric heating theories. We separate the calibrated MWA solar dynamic spectra into a slowly varying and an impulsive, and hence non-thermal, component. We demonstrate that Gaussian mixtures modeling can be used to robustly model the latter and we estimate the flux density distribution as well as the prevalence of impulsive non-thermal emission in the frequency-time plane. Evidence for presence of non-thermal emission at levels down to $\sim$0.2 SFU (1 SFU = 10$^4$ Jy) is reported, making them the weakest reported emissions of this nature. Our work shows the fractional occupancy of the non-thermal impulsive emission to lie in the 17--45\% range during a period of medium solar activity. We also find that the flux density radiated in the impulsive non-thermal emission is very similar in strength to that of the slowly varying component, dominated by the thermal bremsstrahlung. Such significant prevalence and strength of the weak impulsive non-thermal emission has not been appreciated earlier.

LS I +61${^\circ}$303 is a $\gamma$-ray emitting X-ray binary with periodic radio outbursts with time scales of one month. Previous observations have revealed microflares superimposed on these large outbursts with periods ranging from a few minutes to hours. This makes LS I +61${^\circ}$303, along with Cyg X-1, the only TeV emitting X-ray binary exhibiting radio microflares. To further investigate these microflaring activity in LS I +61${^\circ}$303 we observed the source with the 100-m Effelsberg radio telescope at 4.85, 8.35, and 10.45 GHz and performed timing analysis on the obtained data. Radio oscillations of 15 hours time scales are detected at all three frequencies. We also compare the spectral index evolution of radio data to that of the photon index of GeV data observed by Fermi-LAT. We conclude that the observed QPO could result from multiple shocks in a jet.

We study the role of turbulence in magnetic reconnection, within the framework of magneto-hydrodynamics, using three-dimensional direct numerical simulations. For small turbulent intensity we find that the reconnection rate obeys Sweet-Parker scaling. For large enough turbulent intensity reconnection rate departs significantly from Sweet-Parker behaviour, becomes almost a constant as a function of the Lundquist number. We further study energisation of test-particles in the same setup. We find that the speed of the energised particles obeys a Maxwellian distribution, whose variance also obeys Sweet-Parker scaling for small turbulent intensity but depends weakly on the Lundquist number for large turbulent intensity. Furthermore, the variance is found to increase with the strength of the reconnecting magnetic field.

Low radio frequency solar observations using the Murchison Widefield Array have recently revealed the presence of numerous weak, short-lived and narrow-band emission features, even during moderately quiet solar conditions. These non-thermal features occur at rates of many thousands per hour in the 30.72 MHz observing bandwidth, and hence, necessarily require an automated approach for their detection and characterization. Here, we employ continuous wavelet transform using a mother Ricker wavelet for feature detection from the dynamic spectrum. We establish the efficacy of this approach and present the first statistically robust characterization of the properties of these features. In particular, we examine distributions of their peak flux densities, spectral spans, temporal spans and peak frequencies. We can reliably detect features weaker than 1 SFU, making them, to the best of our knowledge, the weakest bursts reported in literature. The distribution of their peak flux densities follows a power law with an index of -2.23 in the 12-155 SFU range, implying that they can provide an energetically significant contribution to coronal and chromospheric heating. These features typically last for 1-2 seconds and possess bandwidths of about 4-5 MHz. Their occurrence rate remains fairly flat in the 140-210 MHz frequency range. At the time resolution of the data, they appear as stationary bursts, exhibiting no perceptible frequency drift. These features also appear to ride on a broadband background continuum, hinting at the likelihood of them being weak type-I bursts.

We present an X-ray timing analysis of the transient X-ray binary MXB 1658-298, using data obtained from the RXTE and XMM-Newton observatories. We have made 27 new mid eclipse time measurements from observations made during the two outbursts of the source. These new measurements have been combined with the previously known values to study long term changes in orbital period of the binary system. We have found that the mid-eclipse timing record of MXB 1658-298 is quite unusual. The long term evolution of mid-eclipse times indicates an overall orbital period decay with a time scale of -- 6.5(7) x 10^7 year. Over and above this orbital period decay, the O-C residual curve also shows a periodic residual on shorter timescales. This sinusoidal variation has an amplitude of ~9 lt-sec and a period of ~760 d. This is indicative of presence of a third body around the compact X-ray binary. The mass and orbital radius of the third body are estimated to lie in the range, 20.5-26.9 Jupiter mass and 750-860 lt-sec, respectively. If true, then it will be the most massive circumbinary planet and also the smallest period binary known to host a planet.

Sky models have been used in the past to calibrate individual low radio frequency telescopes. Here we generalize this approach from a single antenna to a two element interferometer and formulate the problem in a manner to allow us to estimate the flux density of the Sun using the normalized cross-correlations (visibilities) measured on a low resolution interferometric baseline. For wide field-of-view instruments, typically the case at low radio frequencies, this approach can provide robust absolute solar flux calibration for well characterized antennas and receiver systems. It can provide a reliable and computationally lean method for extracting parameters of physical interest using a small fraction of the voluminous interferometric data, which can be prohibitingly compute intensive to calibrate and image using conventional approaches. We demonstrate this technique by applying it to data from the Murchison Widefield Array and assess its reliability.

The observations of post eruption arcades (PEAs) in X-rays and EUV reveal strong localised brightenings at the loop top regions. The origin of these brightenings and their dynamics is not well understood to date. Here, we study the dynamics of PEAs using one-dimensional hydrodynamic modelling with the focus on the the understanding of the formation of localised brightening.Our findings suggest that these brightenings are the result of collisions between the counter-streaming chromospheric evaporation from both the foot points. We further perform forward modelling of the emission observed in simulated results in various spectral lines observed by the Extreme-Ultraviolet Imaging Telescope aboard Hinode. The forward modelled intensities in various spectral lines are in close agreement with a flare observed in December 17, 2006 by EIS.

We study the impact of the nucleon-$\alpha$ P-wave resonances on neutrino-pair bremsstrahlung. Because of the non-central spin-orbit interaction, these resonances lead to an enhanced contribution to the nucleon spin structure factor for temperatures $T \lesssim 4$ MeV. If the $\alpha$-particle fraction is significant and the temperature is in this range, this contribution is competitive with neutron-neutron bremsstrahlung. This may be relevant for neutrino production in core-collapse supernovae or other dense astrophysical environments. Similar enhancements are expected for resonant non-central nucleon-nucleus interactions.

A sample of 1.53$\times$10$^{9}$ cosmic-ray-induced single muon events has been recorded at 225 meters-water-equivalent using the MINOS Near Detector. The underground muon rate is observed to be highly correlated with the effective atmospheric temperature. The coefficient $\alpha_{T}$, relating the change in the muon rate to the change in the vertical effective temperature, is determined to be 0.428$\pm$0.003(stat.)$\pm$0.059(syst.). An alternative description is provided by the weighted effective temperature, introduced to account for the differences in the temperature profile and muon flux as a function of zenith angle. Using the latter estimation of temperature, the coefficient is determined to be 0.352$\pm$0.003(stat.)$\pm$0.046(syst.).

The Sun has remained a difficult source to image for radio telescopes, especially at the low radio frequencies. Its morphologically complex emission features span a large range of angular scales, emission mechanisms involved and brightness temperatures. In addition, time and frequency synthesis, the key tool used by most radio interferometers to build up information about the source being imaged is not effective for solar imaging, because many of the features of interest are short lived and change dramatically over small fractional bandwidths. Building on the advances in radio frequency technology, digital signal processing and computing, the kind of instruments needed to simultaneously capture the evolution of solar emission in time, frequency, morphology and polarization over a large spectral span with the requisite imaging fidelity, and time and frequency resolution have only recently begun to appear. Of this class of instruments, the Murchison Widefield Array (MWA) is best suited for solar observations. The MWA has now entered a routine observing phase and here we present some early examples from MWA observations.

We present an efficient and general method to compute vortex-pinning interactions - which arise in neutron stars, superconductors, and trapped cold atoms - at arbitrary separations using real-time dynamics. This method overcomes uncertainties associated with matter redistribution by the vortex position and the related choice of ensemble that plague the typical approach of comparing energy differences between stationary pinned and unpinned configurations: uncertainties that prevent agreement in the literature on the sign and magnitude of the vortex-nucleus interaction in the crust of neutron stars. We demonstrate and validate the method with Gross-Pitaevskii-like equations for the unitary Fermi gas, and demonstrate how the technique of adiabatic state preparation with time-dependent simulation can be used to calculate vortex-pinning interactions in fermionic systems such as the vortex-nucleus interaction in the crust of neutron stars.

We formulate a low energy effective theory describing phases of matter that are both solid and superfluid. These systems simultaneously break translational symmetry and the phase symmetry associated with particle number. The symmetries restrict the combinations of terms that can appear in the effective action and the lowest order terms featuring equal number of derivatives and Goldstone fields are completely specified by the thermodynamic free energy, or equivalently by the long-wavelength limit of static correlation functions in the ground state. We show that the underlying interaction between particles that constitute the lattice and the superfluid gives rise to entrainment, and mixing between the Goldstone modes. As a concrete example we discuss the low energy theory for the inner crust of a neutron star, where a lattice of ionized nuclei coexists with a neutron superfluid.

We calculate the screening of the ion-ion potential due to electrons in the presence of a large background magnetic field, at densities of relevance to neutron star crusts. Using the standard approach to incorporate electron screening through the one-loop polarization function, we show that the magnetic field produces important corrections both at short and long distances. In extreme fields, realized in highly magnetized neutron stars called magnetars, electrons occupy only the lowest Landau levels in the relatively low density region of the crust. Here our results show that the screening length for Coulomb interactions between ions can be smaller than the inter-ion spacing. More interestingly, we find that the screening is anisotropic and the screened potential between two static charges exhibits long range Friedel oscillations parallel to the magnetic field. This long-range oscillatory behavior is likely to affect the lattice structure of ions, and can possibly create rod-like structures in the magnetar crusts. We also calculate the imaginary part of the electron polarization function which determines the spectrum of electron-hole excitations and plays a role in damping lattice phonon excitations. We demonstrate that even for modest magnetic fields this damping is highly anisotropic and will likely lead to anisotropic phonon heat transport in the outer neutron star crust.

We report on a new mechanism for heat conduction in the neutron star crust. We find that collective modes of superfluid neutron matter, called superfluid phonons (sPhs), can influence heat conduction in magnetized neutron stars. They can dominate the heat conduction transverse to magnetic field when the magnetic field $B \gsim 10^{13}$ G. At density $\rho \simeq 10^{12}-10^{14} $ g/cm$^3$ the conductivity due to sPhs is significantly larger than that due to lattice phonons and is comparable to electron conductivity when temperature $\simeq 10^8$ K. This new mode of heat conduction can limit the surface anisotropy in highly magnetized neutron stars. Cooling curves of magnetized neutron stars with and without superfluid heat conduction could show observationally discernible differences.

We present a study of the shear modulus of the crystalline color superconducting phase of quark matter, showing that this phase of dense, but not asymptotically dense, quark matter responds to shear stress as a very rigid solid. This phase is characterized by a gap parameter $\Delta$ that is periodically modulated in space and therefore spontaneously breaks translational invariance. We derive the effective action for the phonon fields that describe space- and time-dependent fluctuations of the crystal structure formed by $\Delta$, and obtain the shear modulus from the coefficients of the spatial derivative terms. Within a Ginzburg-Landau approximation, we find shear moduli which are 20 to 1000 times larger than those of neutron star crusts. This phase of matter is thus more rigid than any known material in the universe, but at the same time the crystalline color superconducting phase is also superfluid. These properties raise the possibility that the presence of this phase within neutron stars may have distinct implications for their phenomenology. For example, (some) pulsar glitches may originate in crystalline superconducting neutron star cores.

We calculate the shear modulus of crystalline color superconducting quark matter, showing that this phase of dense, but not asymptotically dense, three-flavor quark matter responds to shear stress like a very rigid solid. To evaluate the shear modulus, we derive the low energy effective Lagrangian that describes the phonons that originate from the spontaneous breaking of translation invariance by the spatial modulation of the gap parameter $\Delta$. These massless bosons describe space- and time-dependent fluctuations of the crystal structure and are analogous to the phonons in ordinary crystals. The coefficients of the spatial derivative terms of the phonon effective Lagrangian are related to the elastic moduli of the crystal; the coefficients that encode the linear response of the crystal to a shearing stress define the shear modulus. We analyze the two particular crystal structures which are energetically favored over a wide range of densities, in each case evaluating the phonon effective action and the shear modulus up to order $\Delta^2$ in a Ginzburg-Landau expansion, finding shear moduli which are 20 to 1000 times larger than those of neutron star crusts. The crystalline color superconducting phase has long been known to be a superfluid -- by picking a phase its order parameter breaks the quark-number $U(1)_B$ symmetry spontaneously. Our results demonstrate that this superfluid phase of matter is at the same time a rigid solid. We close with a rough estimate of the pinning force on the rotational vortices which would be formed embedded within this rigid superfluid upon rotation. Our results raise the possibility that (some) pulsar glitches could originate within a quark matter core deep within a neutron star.

Cold three-flavor quark matter at large (but not asymptotically large) densities may exist as a crystalline color superconductor. We explore this possibility by calculating the gap parameter Delta and free energy Omega(Delta) for possible crystal structures within a Ginzburg-Landau approximation, evaluating Omega(Delta) to order Delta^6. We develop a qualitative understanding of what makes a crystal structure stable, and find two structures with particularly large values of Delta and the condensation energy, within a factor of two of those for the CFL phase known to characterize QCD at asymptotically large densities. The robustness of these phases results in their being favored over wide ranges of density and though it also implies that the Ginzburg-Landau approximation is not quantitatively reliable, previous work suggests that it can be trusted for qualitative comparisons between crystal structures. We close with a look ahead at the calculations that remain to be done in order to make contact with observed pulsar glitches and neutron star cooling.

We analyze and compare candidate crystal structures for the crystalline color superconducting phase that may arise in cold, dense but not asymptotically dense, three-flavor quark matter. We determine the gap parameter Delta and free energy Omega(Delta) for many possible crystal structures within a Ginzburg-Landau approximation, evaluating Omega(Delta) to order Delta^6. In contrast to the two-flavor case, we find a positive Delta^6 term and hence an Omega(Delta) that is bounded from below for all the structures that we analyze. This means that we are able to evaluate Delta and Omega as a function of the splitting between Fermi surfaces for all the structures we consider. We find two structures with particularly robust values of Delta and the condensation energy, within a factor of two of those for the CFL phase which is known to characterize QCD at asymptotically large densities. The robustness of these phases results in their being favored over wide ranges of density. However, it also implies that the Ginzburg-Landau approximation is not quantitatively reliable. We develop qualitative insights into what makes a crystal structure favorable, and use these to winnow the possibilities. The two structures that we find to be most favorable are both built from condensates with face-centered cubic symmetry: in one case, the <ud> and <us> condensates are separately face centered cubic; in the other case <ud> and <us> combined make up a face centered cube.