Departures from standard spherically symmetric solar models, in the form of perturbations such as global and local-scale flows and structural asphericities, result in the splitting of eigenfrequencies in the observed spectrum of solar oscillations. Drawing from prevalent ideas in normal-mode coupling theory in geophysical literature, we devise a procedure that enables the computation of sensitivity kernels for general Lorentz stress fields in the Sun. Mode coupling due to any perturbation requires careful consideration of self- and cross-coupling of multiplets. Invoking the isolated-multiplet approximation allows for limiting the treatment to purely self-coupling, requiring significantly less computational resources. We identify the presence of such isolated multiplets under the effect of Lorentz stresses in the Sun. Currently, solar missions allow precise measurements of self-coupling of multiplets via "$a$-coefficients" and the cross-spectral correlation signal which enables the estimation of the "structure coefficients". We demonstrate the forward problem for both self-coupling ($a$-coefficients) and cross-coupling (structure coefficients). In doing so, we plot the self-coupling kernels and estimate $a$-coefficients arising from a combination of deep-toroidal and surface-dipolar axisymmetric fields. We also compute the structure coefficients for an arbitrary general magnetic field (real and solenoidal) and plot the corresponding "splitting function", a convenient way to visualize the splitting of multiplets under 3D internal perturbations. The results discussed in this paper pave the way to formally pose an inverse problem, and infer solar internal magnetic fields.

Inferring interior properties of the Sun from photospheric measurements of the seismic wavefield constitutes the helioseismic inverse problem. Deviations in seismic measurements (such as wave travel times) from their fiducial values estimated for a given model of the solar interior imply that the model is inaccurate. Contemporary inversions in local helioseismology assume that properties of the solar interior are linearly related to measured travel-time deviations. It is widely known, however, that this assumption is invalid for sunspots and active regions, and likely for supergranular flows as well. Here, we introduce nonlinear optimization, executed iteratively, as a means of inverting for the sub-surface structure of large-amplitude perturbations. Defining the penalty functional as the $L_2$ norm of wave travel-time deviations, we compute the the total misfit gradient of this functional with respect to the relevant model parameters %(only sound speed in this case) at each iteration around the corresponding model. The model is successively improved using either steepest descent, conjugate gradient, or quasi-Newton limited-memory BFGS. Performing nonlinear iterations requires privileging pixels (such as those in the near-field of the scatterer), a practice not compliant with the standard assumption of translational invariance. Measurements for these inversions, although similar in principle to those used in time-distance helioseismology, require some retooling. For the sake of simplicity in illustrating the method, we consider a 2-D inverse problem with only a sound-speed perturbation.

Sunspots are prominent manifestations of solar magnetoconvection and imaging their subsurface structure is an outstanding problem of wide physical importance. Travel times of seismic waves that propagate through these structures are typically used as inputs to inversions. Despite the presence of strongly anisotropic magnetic waveguides, these measurements have always been interpreted in terms of changes to isotropic wavespeeds and flow-advection related Doppler shifts. Here, we employ PDE-constrained optimization to determine the appropriate parameterization of the structural properties of the magnetic interior. Seven different wavespeeds fully characterize helioseismic wave propagation: the isotropic sound speed, a Doppler-shifting flow-advection velocity and an anisotropic magnetic velocity. The structure of magnetic media is sensed by magnetoacoustic slow and fast modes and Alfv\'{e}n waves, each of which propagates at a different wavespeed. We show that even in the case of weak magnetic fields, significant errors may be incurred if these anisotropies are not accounted for in inversions. Translation invariance is demonstrably lost. These developments render plausible the accurate seismic imaging of magnetoconvection in the Sun.

Galaxies observed today are likely to have evolved from density perturbations in the early universe. Perturbations that exceeded some critical threshold are conjectured to have undergone gravitational collapse to form primordial black holes (PBHs) at a range of masses. Such PBHs serve as candidates for cold dark matter and their detection would shed light on conditions in the early universe. Here we propose a mechanism to search for transits of PBHs through/nearby Earth by studying the associated seismic waves. Using a spectral-element method, we simulate and visualize this seismic wave field in Earth's interior. We predict the emergence of two unique signatures, namely, a wave that would arrive almost simultaneously everywhere on Earth's free surface and the excitation of unusual spheroidal modes with a characteristic frequency-spacing in free oscillation spectra. These qualitative characteristics are unaffected by the speed or proximity of the PBH trajectory. The seismic energy deposited by a proximal ${M^{PBH} = 10^{15}}$ g PBH is comparable to a magnitude $M_w=4$ earthquake. The non-seismic collateral damage due to the actual impact of such small PBHs with Earth would be negligible. Unfortunately, the expected collision rate is very low even if PBHs constituted all of dark matter, at ${\sim 10^{-7} {yr}^{-1}}$, and since the rate scales as ${1/M^{PBH}}$, fortunately encounters with larger, Earth-threatening PBHs are exceedingly unlikely. However, the rate at which non-colliding close encounters of PBHs could be detected by seismic activity alone is roughly two orders of magnitude larger --- that is once every hundred thousand years --- than the direct collision rate.

For a given {\it misfit function}, a specified optimality measure of a model, its gradient describes the manner in which one may alter properties of the system to march towards a stationary point. The adjoint method, arising from partial-differential-equation-constrained optimization, describes a means of extracting derivatives of a misfit function with respect to model parameters through finite computation. It relies on the accurate calculation of wavefields that are driven by two types of sources, namely the average wave-excitation spectrum, resulting in the {\it forward wavefield}, and differences between predictions and observations, resulting in an {\it adjoint wavefield}. All sensitivity kernels relevant to a given measurement emerge directly from the evaluation of an interaction integral involving these wavefields. The technique facilitates computation of sensitivity kernels (Fr\'{e}chet derivatives) relative to three-dimensional heterogeneous background models, thereby paving the way for non-linear iterative inversions. An algorithm to perform such inversions using as many observations as desired is discussed.

We report on the properties of pre-main-sequence objects in the Taurus molecular clouds as observed in 7 mid- and far-infrared bands with the Spitzer Space Telescope. There are 215 previously-identified members of the Taurus star-forming region in our ~44 square degree map; these members exhibit a range of Spitzer colors that we take to define young stars still surrounded by circumstellar dust (noting that ~20% of the bonafide Taurus members exhibit no detectable dust excesses). We looked for new objects in the survey field with similar Spitzer properties, aided by extensive optical, X-ray, and ultraviolet imaging, and found 148 candidate new members of Taurus. We have obtained follow-up spectroscopy for about half the candidate sample, thus far confirming 34 new members, 3 probable new members, and 10 possible new members, an increase of 15-20% in Taurus members. Of the objects for which we have spectroscopy, 7 are now confirmed extragalactic objects, and one is a background Be star. The remaining 93 candidate objects await additional analysis and/or data to be confirmed or rejected as Taurus members. Most of the new members are Class II M stars and are located along the same cloud filaments as the previously-identified Taurus members. Among non-members with Spitzer colors similar to young, dusty stars are evolved Be stars, planetary nebulae, carbon stars, galaxies, and AGN.