Search

Variational empirical Bayes (VEB) methods provide a practically attractive approach to fitting large, sparse, multiple regression models. These methods usually use coordinate ascent to optimize the variational objective function, an approach known as coordinate ascent variational inference (CAVI). Here we propose alternative optimization approaches based on gradient-based (quasi-Newton) methods, which we call gradient-based variational inference (GradVI). GradVI exploits a recent result from Kim et. al. [arXiv:2208.10910] which writes the VEB regression objective function as a penalized regression. Unfortunately the penalty function is not available in closed form, and we present and compare two approaches to dealing with this problem. In simple situations where CAVI performs well, we show that GradVI produces similar predictive performance, and GradVI converges in fewer iterations when the predictors are highly correlated. Furthermore, unlike CAVI, the key computations in GradVI are simple matrix-vector products, and so GradVI is much faster than CAVI in settings where the design matrix admits fast matrix-vector products (e.g., as we show here, trendfiltering applications) and lends itself to parallelized implementations in ways that CAVI does not. GradVI is also very flexible, and could exploit automatic differentiation to easily implement different prior families. Our methods are implemented in an open-source Python software, GradVI (available from https://github.com/stephenslab/gradvi ).

The integration of machine learning techniques into Inertial Confinement Fusion (ICF) simulations has emerged as a powerful approach for enhancing computational efficiency. By replacing the costly Non-Local Thermodynamic Equilibrium (NLTE) model with machine learning models, significant reductions in calculation time have been achieved. However, determining how to optimize machine learning-based NLTE models in order to match ICF simulation dynamics remains challenging, underscoring the need for physically relevant error metrics and strategies to enhance model accuracy with respect to these metrics. Thus, we propose novel physics-informed transformations designed to emphasize energy transport, use these transformations to establish new error metrics, and demonstrate that they yield smaller errors within reduced principal component spaces compared to conventional transformations., Comment: 10 pages, 7 figures

Software updates are essential to enhance security, fix bugs, and add better features to existing software. However, while some users comply and update their systems upon notification, non-compliance is common. Delaying or ignoring updates leaves systems exposed to security vulnerabilities. Despite research efforts, users' noncompliance behavior with software updates is still prevalent. In this study, we explored how psychological factors influence users' perception and behavior toward software updates. In addition, we proposed a model to assess the security risk score associated with delaying software updates. We conducted a user study with Windows OS users to explore how information about potential vulnerabilities and risk scores influence their behavior. Furthermore, we also studied the influence of demographic factors such as gender on the users' decision-making process for software updates. Our results showed that psychological traits, such as knowledge, awareness, and experience, impact users' decision-making about software updates. To increase users' compliance, providing a risk score for not updating their systems and information about vulnerabilities statistically significantly increased users' willingness to update their systems. Additionally, our results indicated no statistically significant difference in male and female users' responses in terms of concerns about securing their systems. The implications of this study are relevant for software developers and manufacturers as they can use this information to design more effective software update notification messages. Highlighting potential risks and corresponding risk scores in future software updates can motivate users to act promptly to update the systems in a timely manner, which can ultimately improve the overall security of the system., Comment: 16 pages, 10 tables, 7 figures

In this note, we present examples showing that several natural ways of constructing lattices from error-correcting codes do not in general yield a correspondence between minimum-weight non-zero codewords and shortest non-zero lattice vectors. From these examples, we conclude that the main results in two works of Vl\u{a}du\c{t} (Moscow J. Comb. Number Th., 2019 and Discrete Comput. Geom., 2021) on constructing lattices with exponential kissing number from error-correcting codes are invalid. Exhibiting a family of lattices with exponential kissing number therefore remains an open problem.

Earth's magnetotail, and the plasma sheet embedded in it, is a highly dynamic region that is coupled to both the solar wind and to the inner magnetosphere. As a consequence of this coupling, the plasma sheet undergoes explosive energy releases in the form of substorms. One consequence of this energy release is heating of thermal electrons and acceleration of energetic (non-thermal) electrons. The upper-energy limit as well as the spatial scale size of the electron acceleration regions remain mysteries in magnetotail physics because current missions can effectively only offer us a single-point glimpse into the numerous magnetotail phenomena ranging from electron- to global-scales. These energetic electrons can provide a significant source of seed electrons for the Van Allen Radiation belts. Here we use a unique approach to study relativistic plasma sheet electron acceleration. We combine high-altitude Magnetospheric Multiscale (MMS) mission observations with low-altitude Electron Losses and Fields Investigation (ELFIN) observations, to quantify the upper-energy extent and radial scale of a burst of plasma sheet electrons that mapped to 33 Earth radii. The plasma sheet locally accelerated an intense mesoscale burst of 3 MeV electrons -- far higher and more intense than the outer Van Allen radiation belt -- and scattered them into the atmospheric loss cone. High-altitude observations Earthward of the burst at 17 Earth radii showed only the usual substorm activity signatures -- demonstrating that this burst was 1) intense, 2) localized to the far magnetotail, and 3) likely accelerated by a very efficient and rapid mechanism.

Magnetic fields (B-fields) are ubiquitous in the interstellar medium (ISM), and they play an essential role in the formation of molecular clouds and subsequent star formation. However, B-fields in interstellar environments remain challenging to measure, and their properties typically need to be inferred from dust polarization observations over multiple physical scales. In this work, we seek to use a recently proposed approach called the Velocity Gradient Technique (VGT) to study B-fields in star-forming regions and compare the results with dust polarization observations in different wavelengths. The VGT is based on the anisotropic properties of eddies in magnetized turbulence to derive B-field properties in the ISM. We investigate that this technique is synergistic with dust polarimetry when applied to a turbulent diffused medium for the purpose of measuring its magnetization. Specifically, we use the VGT on molecular line data toward the NGC~1333 star-forming region ($\rm ^{12}CO$, $\rm ^{13}CO$, $\rm C^{18}O$, and $\rm N_{2}H^{+}$), and we compare the derived B-field properties with those inferred from 214 and 850~$\mu$m dust polarization observations of the region using SOFIA/HAWC+ and JCMT/POL-2, respectively. We estimate both the inclination angle and the 3D Alfv\'enic Mach Number $M_A$ from the molecular line gradients. Crucially, testing this technique on gravitationally bound, dynamic, and turbulent regions, and comparing the results with those obtained from polarization observations at different wavelength, such as the plane-of-the-sky field orientation, is an important test on the applicability of the VGT in various density regimes of the ISM., Comment: 16 pages, 6 figures, Accepted for publication in ApJ

As the United Kingdom wind fleet increases in size, wind lulls and slews will increasingly challenge the stability of its electricity system. The paper describes the use of models based on real time records and including solar slews, to investigate the most extreme wind variations likely to be encountered in future, enabling strategies to be devised to mitigate them. Wind lulls are surprisingly frequent, occasionally lasting a week or more, and are always likely to be beyond the capabilities of stored or imported electrical energy to mitigate them. The models indicate that there will be a continuing need for gas powered generation to mitigate wind lulls. Currently, Combined Cycle Gas Turbines (CCGTs) provide most of the dispatchable generation. However, CCGTs are not sufficiently fast acting to cope with the wind and solar slews anticipated in future. The paper suggests that a range of already proven fast-acting sources of dispatchable generation, including Open Cycle Gas Turbines (OCGTs), Internal Combustion Gas-Fired Reciprocating engines (ICGRs) and stored electrical energy systems, should be capable of coping with the largest wind and solar slews likely to be encountered up to the year 2035. Examples are given of the recent introduction of these fast-acting sources of generation which, it is suggested, will progressively replace CCGTs as the wind and solar fleets increase in size. Moreover, we see the pattern of recent investments, summarised in the paper, as a good indication of likely future investments, with OCGT investments mainly serving the 440 kV grid, and ICGRs and stored electrical energy more local networks., Comment: 13 pages, 8 figures, 3 tables

Across multiple age groups and academic disciplines, research has shown that incorporating multimodal signs into instruction can enhance student learning. For instance, gesture might be synchronized with speech or written signs. The present study reports on the semiotic resources used by a teacher-researcher to support kindergarten students' developing understanding of the equal sign. An analysis of seven classroom lessons designed to advance students' understanding of the equal sign and equations in various forms revealed four categories of language and two categories of gesture used by the teacher-researcher that were specifically related to mathematical equivalence. Our findings contribute to understanding the role of multimodal signs in supporting instruction on mathematical equivalence. [This is the online first version of an article published in "Mathematics Education Research Journal."]

We use 1-4" (300-1200 au) resolution 12CO(2-1) data from the MASSES (Mass Assembly of Stellar Systems and their Evolution with the SMA) project to measure the projected opening angles of 46 protostellar outflows in the Perseus Molecular Cloud, 37 of which are measured with sufficiently high confidence to use in further analysis. We find that there is a statistically significant difference in the distributions of outflow opening angles for Class 0 and Class I outflows, with a distinct lack of both wide-angle Class 0 outflows and highly collimated Class I outflows. Synthesizing our results with several previous studies, we find that outflows widen with age through the Class 0 stage but do not continue to widen in the Class I stage. The maximum projected opening angle reached is approximately 90 degrees +/- 20 degrees, with the transition between widening and remaining constant occurring near the boundary between the Class 0 and Class I phases of evolution. While the volume fractions occupied by these outflows are no more than a few tens of percent of the total core volume, at most, recent theoretical work suggests outflows may still be capable of playing a central role in setting the low star formation efficiencies of 25%-50% observed on core scales., Comment: 35 pages, accepted by MNRAS

The formation of the massive stars, and in particular, the role that the magnetic fields play in their early evolutionary phase is still far from being completely understood. Here, we present Atacama Large Millimeter/Submillimeter Array (ALMA) 1.2 mm full polarized continuum, and H$^{13}$CO$^+$(3$-$2), CS(5$-$4), and HN$^{13}$C(3$-$2) line observations with a high angular resolution ($\sim$0.4$''$ or 1100 au). In the 1.2 mm continuum emission, we reveal a dusty envelope surrounding the massive protostars, IRAS16547-E and IRAS16547-W, with dimensions of $\sim$10,000 au. This envelope has a bi-conical structure likely carved by the powerful thermal radio jet present in region. The magnetic fields vectors follow very-well the bi-conical envelope. The polarization fraction is $\sim$2.0\% in this region. Some of these vectors seem to converge to IRAS 16547-E, and IRAS 16547-W, the most massive protostars. Moreover, the velocity fields revealed from the spectral lines H$^{13}$CO$^+$(3$-$2), and HN$^{13}$C(3$-$2) show velocity gradients with a good correspondence with the magnetic fields, that maybe are tracing the cavities of molecular outflows or maybe in some parts infall. We derived a magnetic field strength in some filamentary regions that goes from 2 to 6.1\,mG. We also find that the CS(5$-$4) molecular line emission reveals multiple outflow cavities or bow-shocks with different orientations, some of which seem to follow the NW-SE radio thermal jet., Comment: Accepted by the Astrophysical Journal, 13 pages

We expand on recent exciting work of Debris-Alazard, Ducas, and van Woerden [Transactions on Information Theory, 2022], which introduced the notion of basis reduction for codes, in analogy with the extremely successful paradigm of basis reduction for lattices. We generalize DDvW's LLL algorithm and size-reduction algorithm from codes over $\mathbb{F}_2$ to codes over $\mathbb{F}_q$, and we further develop the theory of proper bases. We then show how to instantiate for codes the BKZ and slide-reduction algorithms, which are the two most important generalizations of the LLL algorithm for lattices. Perhaps most importantly, we show a new and very efficient basis-reduction algorithm for codes, called full backward reduction. This algorithm is quite specific to codes and seems to have no analogue in the lattice setting. We prove that this algorithm finds vectors as short as LLL does in the worst case (i.e., within the Griesmer bound) and does so in less time. We also provide both heuristic and empirical evidence that it outperforms LLL in practice, and we give a variant of the algorithm that provably outperforms LLL (in some sense) for random codes. Finally, we explore the promise and limitations of basis reduction for codes. In particular, we show upper and lower bounds on how ``good'' of a basis a code can have, and we show two additional illustrative algorithms that demonstrate some of the promise and the limitations of basis reduction for codes.

We present a study connecting the physical properties of protostellar envelopes to the morphology of the envelope-scale magnetic field. We use the ALMA polarization observations of 55 young prtostars at 0.87 mm on $\sim400-3000$ au scales from the {\em B}-field Orion Protostellar Survey (BOPS) to infer the envelope-scale magnetic field and both dust and gas emission on comparable scales to measure the envelope properties. We find that the protostellar envelopes with compact polarized dust emission tend to have lower envelope masses, than the sources with more extended envelopes. We also find that protostars showing hourglass-field morphologies tend to have lower velocity dispersions in their envelopes, whereas systems with spiral-field morphologies have higher velocity dispersion. Combining with the disk properties taken from the Orion VLA/ALMA Nascent Disk and Multiplicity (VANDAM) survey, we connect envelope properties to fragmentation. Our results suggest that envelope mass may not correlate with fragmentation, whereas turbulence appears to promote fragmentation. On the other hand, we find that fragmentation is suppressed in systems with pinched magnetic fields, suggesting that the magnetic field play a role on providing additional support against gravitational collapse, and the formation of an hourglass-like field may coincide with enhanced magnetic braking that removes angular momentum and hinders the formation of embedded disks. Nevertheless, significant misalignment between magnetic field and outflow axes tends to reduce magnetic braking, leading to the formation of larger disks., Comment: 25 pages, 6 figures, submitted to ApJ

The contribution of the magnetic field to the formation of high-mass stars is poorly understood. We report the high-angular resolution ($\sim0.3^{\prime\prime}$, 870 au) map of the magnetic field projected on the plane of the sky (B$_\mathrm{POS}$) towards the high-mass star forming region G333.46$-$0.16 (G333), obtained with the Atacama Large Millimeter/submillimeter Array (ALMA) at 1.2 mm as part of the Magnetic Fields in Massive Star-forming Regions (MagMaR) survey. The B$_\mathrm{POS}$ morphology found in this region is consistent with a canonical ``hourglass'' which suggest a dynamically important field. This region is fragmented into two protostars separated by $\sim1740$ au. Interestingly, by analysing H$^{13}$CO$^{+}$ ($J=3-2$) line emission, we find no velocity gradient over the extend of the continuum which is consistent with a strong field. We model the B$_\mathrm{POS}$, obtaining a marginally supercritical mass-to-flux ratio of 1.43, suggesting an initially strongly magnetized environment. Based on the Davis-Chandrasekhar-Fermi method, the magnetic field strength towards G333 is estimated to be 5.7 mG. The absence of strong rotation and outflows towards the central region of G333 suggests strong magnetic braking, consistent with a highly magnetized environment. Our study shows that despite being a strong regulator, the magnetic energy fails to prevent the process of fragmentation, as revealed by the formation of the two protostars in the central region.

We classify the finite dimensional representations of the quantum symmetric pair coideal subalgebra $B_{\mathbf c}$ of type $DII$ corresponding to the symmetric pair $(so(2N),so(2N-1))$. For $B_{\mathbf c}$ defined over an arbitrary field $k$ and $q\in k$ not a root of unity we establish a one-to-one correspondence between finite dimensional, simple $B_{\mathbf c}$-modules and dominant integral weights for $so(2N-1)$. We use specialisation to show that the category of finite dimensional $B_{\mathbf c}$-modules is semisimple if $\mathrm{char}(k)=0$ and $q$ is transcendental over ${\mathbb Q}$. In this case the characters of simple $B_{\mathbf c}$-modules are given by Weyl's character formula. This means in particular that the quantum symmetric pair of type $DII$ can be used to obtain Gelfand-Tsetlin bases for irreducible representations of the Drinfeld-Jimbo quantum group $U_q(so(2N))$., Comment: 36 pages

Recent (sub)millimeter polarization observations of protoplanetary disks reveal toroidally aligned, effectively prolate dust grains large enough (at least ~100 $\mu$m) to efficiently scatter millimeter light. The alignment mechanism for these grains remains unclear. We explore the possibility that gas drag aligns grains through gas-dust relative motion when the grain's center of mass is offset from its geometric center, analogous to a badminton birdie's alignment in flight. A simple grain model of two non-identical spheres illustrates how a grain undergoes damped oscillations from flow-induced restoring torques which align its geometric center in the flow direction relative to its center of mass. Assuming specular reflection and subsonic flow, we derive an analytical equation of motion for spheroids where the center of mass can be shifted away from the spheroid's geometric center. We show that a prolate or an oblate grain can be aligned with the long axis parallel to the gas flow when the center of mass is shifted along that axis. Both scenarios can explain the required effectively prolate grains inferred from observations. Application to a simple disk model shows that the alignment timescales are shorter than or comparable to the orbital time. The grain alignment direction in a disk depends on the disk (sub-)structure and grain Stokes number (St) with azimuthal alignment for large St grains in sub-Keplerian smooth gas disks and for small St grains near the gas pressure extrema, such as rings and gaps., Comment: 21 pages, 12 figures, accepted by MNRAS

The central regions of the Milky Way constitute a unique laboratory for a wide swath of astrophysical studies, consequently the inner $\sim$400 pc has been the target of numerous large surveys at all accessible wavelengths. In this paper we present a catalog of sources at 25 and 37 $\mu$m located within all of the regions observed with the SOFIA/FORCAST instrument in the inner $\sim$200 pc of the Galaxy. The majority of the observations were obtained as part of the SOFIA Cycle 7 Galactic Center Legacy program survey, which was designed to complement the Spitzer/MIPS 24 $\mu$m catalog in regions saturated in the MIPS observations. Due to the wide variety of source types captured by our observations at 25 and 37 $\mu$m, we do not limit the FORCAST source catalog to unresolved point sources, or treat all sources as if they are point-like sources. The catalog includes all detectable sources in the regions, resulting in a catalog of 950 sources, including point sources, compact sources, and extended sources. We also provide the user with metrics to discriminate between the source types., Comment: 29 pages, 13 figures, Accepted to ApJ

Global numerical weather models are starting to resolve atmospheric moist convection which comes with a critical need for observational constraints. One avenue for such constraints is spaceborne radar which tend to operate at three wavelengths, Ku-, Ka- and W-band. Many studies of deep convection in the past have primarily leveraged Ku-band because it is less affected by attenuation and multiple scattering. However, future spaceborne radar missions might not contain a Ku-band radar and thus considering the view of convection from Ka-band or W-band compared to the Ku-band would be useful. This study examines a coincident dataset between the Global Precipitation Measurement (GPM) Mission and CloudSat as well as the entire GPM record to compare convective characteristics across various wavelengths within deep convection. We find that W-band reflectivity (Z) tends to maximize near the Ku-band defined echo-top while Ka-band often maximizes 4-5 km below. The height of the maximum Z above the melting level for W-band does not linearly relate to the Ku-band maximum. However, using the full GPM record the Ka-band 30 dBZ echo-tops can be linearly related to the Ku-band 40 dBZ echo-top with an $R^2$ of 0.62 and a root mean squared error of about 1 km. The spatial distribution of echo-tops from Ka-band corresponds well to the Ku-band echo-tops, highlighting regions of relatively large ice water path. This paper suggests that Ka-band only missions, like NASA's Investigation for Convective Updrafts, should be able to characterize global convection in a similar manner to a Ku-band system.

Since ZMap's debut in 2013, networking and security researchers have used the open-source scanner to write hundreds of research papers that study Internet behavior. In addition, ZMap has been adopted by the security industry to build new classes of enterprise security and compliance products. Over the past decade, much of ZMap's behavior -- ranging from its pseudorandom IP generation to its packet construction -- has evolved as we have learned more about how to scan the Internet. In this work, we quantify ZMap's adoption over the ten years since its release, describe its modern behavior (and the measurements that motivated changes), and offer lessons from releasing and maintaining ZMap for future tools., Comment: 10 pages, 7 figures, in submission at Internet Measurement Conference 2024

We report on ALMA observations of polarized dust emission at 1.2 mm from NGC6334I, a source known for its significant flux outbursts. Between five months, our data show no substantial change in total intensity and a modest 8\% variation in linear polarization, suggesting a phase of stability or the conclusion of the outburst. The magnetic field, inferred from this polarized emission, displays a predominantly radial pattern from North-West to South-East with intricate disturbances across major cores, hinting at spiral structures. Energy analysis of CS$(J=5 \rightarrow 4)$ emission yields an outflow energy of approximately $3.5\times10^{45}$ ergs, aligning with previous interferometric studies. Utilizing the Davis-Chandrasekhar-Fermi method, we determined magnetic field strengths ranging from 1 to 11 mG, averaging at 1.9 mG. This average increases to 4 $\pm 1$ mG when incorporating Zeeman measurements. Comparative analyses using gravitational, thermal, and kinetic energy maps reveal that magnetic energy is significantly weaker, possibly explaining the observed field morphology. We also find that the energy in the outflows and the expanding cometary {\HII} region is also larger than the magnetic energy, suggesting that protostellar feedback maybe the dominant driver behind the injection of turbulence in NGC6334I at the scales sampled by our data. The gas in NGC6334I predominantly exhibits supersonic and trans-Alfvenic conditions, transitioning towards a super-Alfvenic regime, underscoring a diminished influence of the magnetic field with increasing gas density. These observations are in agreement with prior polarization studies at 220 GHz, enriching our understanding of the dynamic processes in high-mass star-forming regions., Comment: Accepted for Publication at the Astrophysical Journal

A global reconfiguration of the magnetotail characterizes substorms. Current sheet thinning, intensification, and magnetic field stretching are defining features of the substorm growth phase and their spatial distributions control the timing and location of substorm onset. Presently, sparse in-situ observations cannot resolve these distributions. A promising approach is to use new substorm magnetic field reconstruction methods based on data mining, termed SST19. Here we compare the SST19 reconstructions to low-altitude ELFIN measurements of energetic particle precipitations to probe the radial profile of the equatorial magnetic field curvature during a 19~August 2022 substorm. ELFIN and SST19 yield a consistent dynamical picture of the magnetotail during the growth phase and capture expected features such as the formation of a thin current sheet and its earthward motion. Furthermore, they resolve a V-like pattern of isotropic electron precipitation boundaries in the time-energy plane, consistent with earlier observations but now over a broad energy range.

We present ALMA Band 7 molecular line observations of the protostars within the VLA 1623 system. We map C$^{17}$O (3 - 2) in the circumbinary disk around VLA 1623A and the outflow cavity walls of the collimated outflow. We further detect red-shifted and blue-shifted velocity gradients in the circumstellar disks around VLA 1623B and VLA 1623W that are consistent with Keplerian rotation. We use the radiative transfer modeling code, pdspy, and simple flared disk models to measure stellar masses of $0.27 \pm 0.03$ M$_\odot$, $1.9^{+0.3}_{-0.2}$ M$_\odot$, and $0.64 \pm 0.06$ M$_\odot$ for the VLA 1623A binary, VLA 1623B, and VLA 1623W, respectively. These results represent the strongest constraints on stellar mass for both VLA 1623B and VLA 1623W, and the first measurement of mass for all stellar components using the same tracer and methodology. We use these masses to discuss the relationship between the young stellar objects (YSOs) in the VLA 1623 system. We find that VLA 1623W is unlikely to be an ejected YSO, as has been previously proposed. While we cannot rule out that VLA 1623W is a unrelated YSO, we propose that it is a true companion star to the VLA 1623A/B system and that the these stars formed in situ through turbulent fragmentation and have had only some dynamical interactions since their inception., Comment: Accepted to A&A; 16 pages, 12 figures

Estimating the sharing of genetic effects across different conditions is important to many statistical analyses of genomic data. The patterns of sharing arising from these data are often highly heterogeneous. To flexibly model these heterogeneous sharing patterns, Urbut et al. (2019) proposed the multivariate adaptive shrinkage (MASH) method to jointly analyze genetic effects across multiple conditions. However, multivariate analyses using MASH (as well as other multivariate analyses) require good estimates of the sharing patterns, and estimating these patterns efficiently and accurately remains challenging. Here we describe new empirical Bayes methods that provide improvements in speed and accuracy over existing methods. The two key ideas are: (1) adaptive regularization to improve accuracy in settings with many conditions; (2) improving the speed of the model fitting algorithms by exploiting analytical results on covariance estimation. In simulations, we show that the new methods provide better model fits, better out-of-sample performance, and improved power and accuracy in detecting the true underlying signals. In an analysis of eQTLs in 49 human tissues, our new analysis pipeline achieves better model fits and better out-of-sample performance than the existing MASH analysis pipeline. We have implemented the new methods, which we call ``Ultimate Deconvolution'', in an R package, udr, available on GitHub.

Center-based early care and education (ECE) is important for promoting positive early development and supporting families by providing child care so parents can work. However, the center-based supply varies substantially in terms of funding sources, indicators of quality, and services offered; and many families experience a lack of equitable access to providers that meet their needs. Using a nationally representative sample of the center-based supply from the 2019 National Survey of Early Care and Education, this study examined patterns in the services center-based providers offer that may facilitate or pose barriers to families' access from a multidimensional lens. Two research questions were investigated using latent class analysis: First, are there patterns of provider accessibility within the center-based supply, and second, are there covariates (i.e., funding, community poverty/ urban density, geographic region) that explain patterns of ECE access? Center-based features selected as indicators were informed by a definition of access comprised of multiple dimensions: availability, affordability, supports children's development, and meets parent's needs. Class enumeration suggested that a three-class model provided the best fit, indicating three patterns of program accessibility: providers that were: "flexible at cost" (25.9%), "free but inflexible" (29.2%), and "somewhat accessible" (44.9%). Additional covariate analyses revealed providers were significantly more or less likely to be classified into the latent classes of access according to community and funding characteristics. Notably, programs receiving Head Start and public preschool funding and those in communities with high poverty and urban density were more likely to be assigned to the "free but inflexible" class. A key component driving ECE utilization is the accessibility of providers, and these findings inform recommendations for future research and policy to better serve the needs of families with young children. [This paper will be published in "Children and Youth Services Review."]

We present the results from the first two years of the Planet Hunters NGTS citizen science project, which searches for transiting planet candidates in data from the Next Generation Transit Survey (NGTS) by enlisting the help of members of the general public. Over 8,000 registered volunteers reviewed 138,198 light curves from the NGTS Public Data Releases 1 and 2. We utilize a user weighting scheme to combine the classifications of multiple users to identify the most promising planet candidates not initially discovered by the NGTS team. We highlight the five most interesting planet candidates detected through this search, which are all candidate short-period giant planets. This includes the TIC-165227846 system that, if confirmed, would be the lowest-mass star to host a close-in giant planet. We assess the detection efficiency of the project by determining the number of confirmed planets from the NASA Exoplanet Archive and TESS Objects of Interest (TOIs) successfully recovered by this search and find that 74% of confirmed planets and 63% of TOIs detected by NGTS are recovered by the Planet Hunters NGTS project. The identification of new planet candidates shows that the citizen science approach can provide a complementary method to the detection of exoplanets with ground-based surveys such as NGTS., Comment: 42 pages, 20 figures, 17 tables. To be published in AJ

Millimeter and sub-millimeter observations of continuum linear dust polarization provide insight into dust grain growth in protoplanetary disks, which are the progenitors of planetary systems. We present the results of the first survey of dust polarization in protoplanetary disks at 870 $\mu$m and 3 mm. We find that protoplanetary disks in the same molecular cloud at similar evolutionary stages can exhibit different correlations between observing wavelength and polarization morphology and fraction. We explore possible origins for these differences in polarization, including differences in dust populations and protostar properties. For RY Tau and MWC 480, which are consistent with scattering at both wavelengths, we present models of the scattering polarization from several dust grain size distributions. These models aim to reproduce two features of the observational results for these disks: (1) both disks have an observable degree of polarization at both wavelengths and (2) the polarization fraction is higher at 3 mm than at 870 $\mu$m in the centers of the disks. For both disks, these features can be reproduced by a power-law distribution of spherical dust grains with a maximum radius of 200 $\mu$m and high optical depth. In MWC 480, we can also reproduce features (1) and (2) with a model containing large grains ($a_{max}$ = 490 $\mu$m ) near the disk midplane and small grains ($a_{max}$ = 140 $\mu$m) above and below the midplane., Comment: 22 pages, 12 figures

We describe a system called Qr-Hint that, given a (correct) target query Q* and a (wrong) working query Q, both expressed in SQL, provides actionable hints for the user to fix the working query so that it becomes semantically equivalent to the target. It is particularly useful in an educational setting, where novices can receive help from Qr-Hint without requiring extensive personal tutoring. Since there are many different ways to write a correct query, we do not want to base our hints completely on how Q* is written; instead, starting with the user's own working query, Qr-Hint purposefully guides the user through a sequence of steps that provably lead to a correct query, which will be equivalent to Q* but may still "look" quite different from it. Ideally, we would like Qr-Hint's hints to lead to the "smallest" possible corrections to Q. However, optimality is not always achievable in this case due to some foundational hurdles such as the undecidability of SQL query equivalence and the complexity of logic minimization. Nonetheless, by carefully decomposing and formulating the problems and developing principled solutions, we are able to provide provably correct and locally optimal hints through Qr-Hint. We show the effectiveness of Qr-Hint through quality and performance experiments as well as a user study in an educational setting., Comment: SIGMOD 2024

The Electrocardiography Brain Perfusion index (EBPi) is a novel electrocardiography (ECG)-based metric that may function as a proxy for cerebral blood flow (CBF). We investigated the spatio-temporal correlation between EBPi and epileptic seizure events. EBPi was computed retrospectively from clinical EEG and ECG data captured previously from 30 epilepsy patients during seizures. Significant EBPi changes were compared temporally with clinically defined ground-truth seizure onset and offset times. We also assessed the spatial correlation between EBPi metrics and clinically defined ground-truth seizure locations. A significant increase in EBPi was detected 10.5 s [-6, 53] (median [95% confidence interval (CI)]) after ground-truth seizure onset, and a significant decrease in EBPi was detected 5 s [-42, 74] (median [95% CI]) after ground-truth seizure offset. EBPi demonstrated a positive predictive value of 61.5% [33.3, 75] (median [95% CI]) and a sensitivity of 57.1% [38.5, 66.7] (median [95% CI]) for the detection of ground truth seizure locations. EBPi signals exhibited a temporal sensitivity to seizure events and in some cases were correlated spatially to the seizure location. Therefore, EBPi, which has been linked to CBF, appears to contain spatio-temporal information related to seizure activity and might have useful application in augmenting EEG data captured during seizures.

Despite the rise of mobile robot deployments in home and work settings, perceived safety of users and bystanders is understudied in the human-robot interaction (HRI) literature. To address this, we present a study designed to identify elements of a human-robot encounter that correlate with observed stress response. Stress is a key component of perceived safety and is strongly associated with human physiological response. In this study a Boston Dynamics Spot and a Unitree Go1 navigate autonomously through a shared environment occupied by human participants wearing multimodal physiological sensors to track their electrocardiography (ECG) and electrodermal activity (EDA). The encounters are varied through several trials and participants self-rate their stress levels after each encounter. The study resulted in a multidimensional dataset archiving various objective and subjective aspects of a human-robot encounter, containing insights for understanding perceived safety in such encounters. To this end, acute stress responses were decoded from the human participants' ECG and EDA and compared across different human-robot encounter conditions. Statistical analysis of data indicate that on average (1) participants feel more stress during encounters compared to baselines, (2) participants feel more stress encountering multiple robots compared to a single robot and (3) participants stress increases during navigation behavior compared with search behavior., Comment: 8 pages, 7 figs, 5 tables