Search

A new coordinate system on the tangent space to planar configurations is introduced to simplify some calculations on central configurations and relative equilibria in the $N$-body problem with a homogeneous potential, which includes the case of Newtonian gravity. These coordinates are applied to some problems on four-body central configurations to illustrate their utility.

Hyper-spectral imaging has recently gained increasing attention for use in different applications, including agricultural investigation, ground tracking, remote sensing and many other. However, the high cost, large physical size and complicated operation process stop hyperspectral cameras from being employed for various applications and research fields. In this paper, we introduce a cost-efficient, compact and easy to use active illumination camera that may benefit many applications. We developed a fully functional prototype of such camera. With the hope of helping with agricultural research, we tested our camera for plant root imaging. In addition, a U-Net model for spectral reconstruction was trained by using a reference hyperspectral camera's data as ground truth and our camera's data as input. We demonstrated our camera's ability to obtain additional information over a typical RGB camera. In addition, the ability to reconstruct hyperspectral data from multi-spectral input makes our device compatible to models and algorithms developed for hyperspectral applications with no modifications required.

The integration of technology into exercise regimens has emerged as a strategy to enhance normal human capabilities and return human motor function after injury or illness by enhancing motor learning and retention. Much research has focused on how active devices, whether confined to a lab or made into a wearable format, can apply forces at set times and conditions to optimize the process of learning. However, the focus on active force production often forces devices to either be confined to simple movements or interventions. As such, in this paper, we investigate how passive device behaviors can contribute to the process of motor learning by themselves. Our approach involves using a wearable resistance (WR) device, which is outfitted with elastic bands, to apply a force field that changes in response to a person's movements while performing exercises. We develop a method to measure the produced forces from the device without impeding the function and we characterize the device's force generation abilities. We then present a study assessing the impact of the WR device on motor learning of proper squat form compared to visual or no feedback. Biometrics such as knee and hip angles were used to monitor and assess subject performance. Our findings indicate that the force fields produced while training with the WR device can improve performance in full-body exercises similarly to a more direct visual feedback mechanism, though the improvement is not consistent across all performance metrics. Through our research, we contribute important insights into the application of passive wearable resistance technology in practical exercise settings., Comment: 8 pages, 9 figures, To be published in IEEE International Conference on Biomedical Robotics and Biomechatronics (BioRob) 2024

Microstate geometries are proposed microstates of black holes which can be described within supergravity. Even though their number may not reproduce the full entropy of black holes with finite-sized horizons, they still offer a glimpse into the microscopic structure of black holes. In this paper we construct a new set of microstate geometries of the supersymmetric D1-D5-P black hole, where the momentum charge is carried by a vector field, as seen from the perspective of six-dimensional supergravity. To aid our construction, we develop an algorithm which solves a complicated partial differential equation using the regularity of the geometries. The new solutions are asymptotically AdS$_3\times S^3$, and have a long, but finite AdS$_2$ throat that caps off without ever developing a horizon. These microstate geometries have a holographic interpretation as coherent superpositions of heavy states in the boundary D1-D5 CFT. We identify the states which are dual to our newly constructed solutions and carry out some basic consistency checks to support our identification., Comment: 27+19 pages, 3 figures, comments welcome

We use families of circular null geodesics as probes of a family of microstate geometries, known as $(1,0,n)$ superstrata. These geometries carry a left-moving momentum wave and the behavior of some of the geodesic probes is very sensitive to this background wave. The left-moving geodesics behave like BPS particles and so can be placed in circular orbits anywhere in the geometry and actually "float" at fixed radius and angle in the three-dimensional "capped BTZ" geometry. The right-moving geodesics behave like non-BPS particles. We show that they provide a simple geometric characterization of the black-hole bound: when the momentum charge of the geometry is below this bound, such geodesics can be placed anywhere, but exceeding the bound, even by a small amount, means these geodesics are restricted to the deep interior of the geometry. We also show that for left-moving string probes, the tidal forces remain comparable with those of global AdS$_3$. Nevertheless, for some of these probes, the "bumps" in the geometry induce an oscillatory mass term and we discuss how this can lead to chaotic scrambling of the state of the string., Comment: 29 pages, 8 figures

We investigate the effects of the twist-2 operator in 2D symmetric orbifold CFTs. The twist operator can join together a twist-$M$ state and a twist-$N$ state, creating a twist-$(M+N)$ state. This process involves three effects: pair creation, propagation, and contraction. We study these effects by using a Bogoliubov ansatz and conformal symmetry. In this multi-wound scenario, pair creation no longer decouples from propagation, in contrast to the previous study where $M=N=1$. We derive equations for these effects, which organize themselves into recursion relations and constraints. Using the recursion relations, we can determine the infinite number of coefficients in the effects through a finite number of inputs. Moreover, the number of required inputs can be further reduced by applying constraints., Comment: 38 pages

We construct a two-parameter four-dimensional non-BPS NS-NS smooth microstate solution that asymptotes to flat spacetime with a linear dilaton in type II superstring theory. From the microscopic point of view, the background is made out of a certain number of decoupled (i.e. $g_s\to 0$) NS5 branes wrapping $T^3\times S^1\times S^1$ with fundamental strings wrapping non-contractable cycles of $S^1\times S^1$ with integer momentum modes along them. We show that perturbative worldsheet theory in this background is given by a null-gauged WZW model. We also show that the consistency of the worldsheet theory imposes non-trivial constraints on the supergravity background., Comment: 32 pages

We previously showed that three-dimensional quadratic diffeomorphisms have anti-integrable (AI) limits that correspond to a quadratic correspondence; a pair of one-dimensional maps. At the AI limit the dynamics is conjugate to a full shift on two symbols. Here we consider a more general AI limit, allowing two parameters of the map to go to infinity. We prove the existence of AI states for each symbol sequence for three cases of the quadratic correspondence: parabolas, ellipses and hyperbolas. A contraction argument gives parameter domains such that this is a bijection, but the correspondence also is observed to apply more generally. We show that orbits of the original map can be obtained by numerical continuation for a volume-contracting case. These results show that periodic AI states evolve into the observed periodic attractors of the diffeomorphism. We also continue a periodic AI state with a symbol sequence chosen so that it continues to an orbit resembling a chaotic attractor that is a 3D version of the classical 2D H\'enon attractor., Comment: 27 pages, 9 figures, link to movie

We argue that the fundamental "atomic objects" in string theory are themelia: extended objects that have 16 supersymmetries locally. We show that all existing smooth horizonless microstate geometries can be seen as bound states of themelia, and we conjecture that all such bound states with suitable KKM charges will give rise to microstate geometries. We also construct the most general themelion with a three-torus isometry and show that it interpolates between superstrata and the super-maze., Comment: 6 pages

The entropy of the three-charge NS5-F1-P black hole in Type IIA string theory comes from the breaking of $N_1$ F1 strings into $N_1 N_5$ little strings, which become independent momentum carriers. In M theory, the little strings correspond to strips of M2 brane that connect pairs of parallel M5 branes separated along the M-theory direction. We show that if one takes into account the backreaction of the M-theory little strings on the M5 branes one obtains a maze-like structure, to which one can add momentum waves. We also show that adding momentum waves to the little strings gives rise to a momentum-carrying brane configuration -- a super-maze -- which locally preserves 16 supercharges. We therefore expect the backreaction of the super-maze to give rise to a new class of horizonless black-hole microstate solutions, which preserve the rotational symmetry of the black-hole horizon and carry $\sqrt{5/6}$ of its entropy., Comment: 29 pages, 7 figures, LaTeX

In the D1D5 CFT the twist operator of order 2 can twist together two copies in the untwisted sector into a single joined copy in the twisted sector. Traditionally, this effect is computed by using the covering map method. Recently, a new method was developed using the Bogoliubov ansatz and conformal symmetry to compute this effect in a toy model of one free boson. In this paper, we use this method with superconformal symmetry to compute the effect of the twist operator in the D1D5 CFT. This may provide more effective tools for computing correlation functions of twist operators in this system., Comment: 33 pages

We study dynamics of a generic quadratic diffeomorphism, a 3D generalization of the planar H\'{e}non map. Focusing on the dissipative, orientation preserving case, we give a comprehensive parameter study of codimension-one and two bifurcations. Periodic orbits, born at resonant, Neimark-Sacker bifurcations, give rise to Arnold tongues in parameter space. Aperiodic attractors include invariant circles and chaotic orbits; these are distinguished by rotation number and Lyapunov exponents. Chaotic orbits include H\'{e}non-like and Lorenz-like attractors, which can arise from period-doubling cascades, and those born from the destruction of invariant circles. The latter lie on paraboloids near the local unstable manifold of a fixed point., Comment: 19 pages, 13 figures

We study the 2D symmetric orbifold CFT of two copies of free bosons. The twist operator can join the two separated copies in the untwisted sector into a joined copy in the twisted sector. Starting with a state with any number of quanta in the untwisted sector, the state in the twisted sector obtained by the action of the twist operator can be computed by using the covering map method. We develop a new method to compute the effect of a twist operator by using the Bogoliubov ansatz and conformal symmetry. This may lead to more efficient tools to compute correlation functions involving twist operators., Comment: 33 pages, 2 figures

We analyse the BPS equations of $\mathcal{N} = (1,0)$ supergravity theory in six dimensions coupled to a vector and tensor multiplet. We show how these BPS equations can be reduced to a set of linear differential equations. This system is triangular in that each layer of equations, while linear, is quadratically sourced by the solutions of the previous layers. We examine several explicit examples and discuss the construction of new families of microstate geometries. We expect that the result presented here will open up new branches of superstrata in which the momentum is encoded in a new class of charge carriers., Comment: 48 pages, 1 figure

The relative equilibria of planar Newtonian $N$-body problem become coorbital around a central mass in the limit when all but one of the masses becomes zero. We prove a variety of results about the coorbital relative equilibria, with an emphasis on the relation between symmetries of the configurations and symmetries in the masses, or lack thereof. We prove that in the $N=4$, $N=6$, and $N=8$ Newtonian coorbital problems there exist symmetric relative equilibria with asymmetric positive masses. This result can be generalized to other homogeneous potentials, and we conjecture similar results hold for larger even numbers of infinitesimal masses. We prove that some equalities of the masses in the $1+4$ and $1+5$ coorbital problems imply symmetry of a class of convex relative equilibria. We also prove there is at most one convex central configuration of the symmetric $1+5$ problem.

All known horizonless black-hole microstate geometries correspond to brane sources that acquire a finite size, and hence break the spherical symmetry of the black hole. We construct, for the first time, solutions with zero horizon area that have the same charges as a three-charge F1-NS5-P Type-IIA black hole and preserve this spherical symmetry. The momentum of these solutions is carried by longitudinal D0-D4 density fluctuations inside the NS5-branes. We argue that these solutions should be interpreted as the long-throat limit of a family of smooth, horizonless microstate geometries, called superstrata, where such geometries degenerate. The existence of these geometries indicates that a finite-size horizon does not appear even in the singular corners of the moduli space of three-charge microstate geometries., Comment: 38 pages plus references, 5 figures

We consider the entanglement entropy of an arbitrary subregion in a system of $N$ non-relativistic fermions in $2+1$ dimensions in Lowest Landau Level (LLL) states. Using the connection of these states to those of an auxiliary $1+1$ dimensional fermionic system, we derive an expression for the leading large-$N$ contribution in terms of the expectation value of the phase space density operator in $1+1$ dimensions. For appropriate subregions the latter can replaced by its semiclassical Thomas-Fermi value, yielding expressions in terms of explicit integrals which can be evaluated analytically. We show that the leading term in the entanglement entropy is a perimeter law with a shape independent coefficient. Furthermore, we obtain analytic expressions for additional contributions from sharp corners on the entangling curve. Both the perimeter and the corner pieces are in good agreement with existing calculations for special subregions. Our results are relevant to the integer quantum Hall effect problem, and to the half-BPS sector of $\mathcal N=4$ Yang Mills theory on $S^3$. In this latter context, the entanglement we consider is an entanglement in target space. We comment on possible implications to gauge-gravity duality., Comment: 57 pages, 5 figures

The two-dimensional $\mathcal{N}=4$ superconformal algebra has a free field realization with four bosons and four fermions. There is an automorphism of the algebra called spectral flow. Under spectral flow, the four fermions are transformed together. In this paper, we study {\it partial spectral flow} where only two of the four fermions are transformed. Partial spectral flow is applied to the D1D5 CFT where a marginal deformation moves the CFT away from the free point. The partial spectral flow is broken by the deformation. We show that this effect can be studied due to a transformation of the deformation which is well-defined under partial spectral flow. As a result in the spectrum, we demonstrate how to compute the second-order energy lift of D1D5P states from their partial spectral flowed states. We find that D1D5P states related by partial spectral flow do not have the same lift in general., Comment: 29 pages, 1 figure

There is growing evidence for physical influence between supermassive black holes and their host galaxies. We present a case study of nearby galaxy NGC 7582, for which we find evidence that galactic substructure plays an important role in affecting the collimation of ionized outflows as well as contributing to the heavy active galactic nucleus (AGN) obscuration. This result contrasts with a simple, small-scale AGN torus model, according to which AGN wind collimation may take place inside the torus itself, at subparsec scale. Using 3D spectroscopy with the MUSE instrument, we probe the kinematics of the stellar and ionized gas components as well as the ionization state of the gas from a combination of emission line ratios. We report for the first time a kinematically distinct core (KDC) in NGC 7582, on a scale of ~600pc. This KDC coincides spatially with dust lanes and starbursting complexes previously observed. We interpret it as a circumnuclear ring of stars and dusty, gas-rich material. We obtain a clear view of the outflowing cones over kpc scales, and demonstrate that they are predominantly photoionized by the central engine. We detect the back cone (behind the galaxy), and confirm previous results of a large nuclear obscuration of both the stellar continuum and HII regions. While we tentatively associate the presence of the KDC to a large-scale bar and/or a minor galaxy merger, we stress the importance of gaining a better understanding of the role of galaxy substructure in controlling the fueling, feedback and obscuration of AGN., Comment: 31 pages, 16 figures (including two animations), accepted for publication in the Astrophysical Journal. For a preview of the animations, see https://noirlab.edu/science/videos/rotating-gas-disk-of-the-ngc-7582-galaxy for Figure 15, and see https://noirlab.edu/science/videos/black-hole-wind-from-the-ngc-7582-galaxy for Figure 16

Central configurations and relative equilibria are an important facet of the study of the $N$-body problem, but become very difficult to rigorously analyze for $N>3$. In this paper we focus on a particular but interesting class of configurations of the 5-body problem: the equilateral pentagonal configurations, which have a cycle of five equal edges. We prove a variety of results concerning central configurations with this property, including a computer-assisted proof of the finiteness of such configurations for any positive five masses with a range of rational-exponent homogeneous potentials (including the Newtonian case and point-vortex model), some constraints on their shapes, and we determine some exact solutions for particular N-body potentials.

The Multilevel Monte Carlo (MLMC) method has proven to be an effective variance-reduction statistical method for Uncertainty Quantification (UQ) in Partial Differential Equation (PDE) models, combining model computations at different levels to create an accurate estimate. Still, the computational complexity of the resulting method is extremely high, particularly for 3D models, which requires advanced algorithms for the efficient exploitation of High Performance Computing (HPC). In this article we present a new implementation of the MLMC in massively parallel computer architectures, exploiting parallelism within and between each level of the hierarchy. The numerical approximation of the PDE is performed using the finite element method but the algorithm is quite general and could be applied to other discretization methods as well, although the focus is on parallel sampling. The two key ingredients of an efficient parallel implementation are a good processor partition scheme together with a good scheduling algorithm to assign work to different processors. We introduce a multiple partition of the set of processors that permits the simultaneous execution of different levels and we develop a dynamic scheduling algorithm to exploit it. The problem of finding the optimal scheduling of distributed tasks in a parallel computer is an NP-complete problem. We propose and analyze a new greedy scheduling algorithm to assign samples and we show that it is a 2-approximation, which is the best that may be expected under general assumptions. On top of this result we design a distributed memory implementation using the Message Passing Interface (MPI) standard. Finally we present a set of numerical experiments illustrating its scalability properties., Comment: 21 pages, 13 figures

Unsupervised learning algorithms are beginning to achieve accuracies comparable to their supervised counterparts on benchmark computer vision tasks, but their utility for practical applications has not yet been demonstrated. In this work, we present a novel application of unsupervised learning to the task of auroral image classification. Specifically, we modify and adapt the Simple framework for Contrastive Learning of Representations (SimCLR) algorithm to learn representations of auroral images in a recently released auroral image dataset constructed using image data from Time History of Events and Macroscale Interactions during Substorms (THEMIS) all-sky imagers. We demonstrate that (a) simple linear classifiers fit to the learned representations of the images achieve state-of-the-art classification performance, improving the classification accuracy by almost 10 percentage points over the current benchmark; and (b) the learned representations naturally cluster into more clusters than exist manually assigned categories, suggesting that existing categorizations are overly coarse and may obscure important connections between auroral types, near-earth solar wind conditions, and geomagnetic disturbances at the earth's surface. Moreover, our model is much lighter than the previous benchmark on this dataset, requiring in the area of fewer than 25\% of the number of parameters. Our approach exceeds an established threshold for operational purposes, demonstrating readiness for deployment and utilization., Comment: 6 pages, 5 figures, 1 table

It was demonstrated that a string probe falling radially within a superstratum geometry would experience tidal forces. These tidal forces were shown to excite the string by converting its kinetic energy into stringy excitations. Using the AdS/CFT correspondence we seek to understand this behavior from the perspective of the dual D1D5 CFT. To study this process we turn on an interaction of the theory which is described by a deformation operator. We start with an initial state which is dual to a graviton probe moving within the superstratum geometry. We then use two deformation operators to compute transition amplitudes between this state and a final state that corresponds to stringy excitations. We show that this amplitude grows as $t^2$ with $t$ being the amount of time for which the deformation operators are turned on. We argue that this process in the CFT is suggestive of the tidal effects experienced by the probe propagating within the dual superstratum geometry., Comment: 21 pages, 4 figures

We study a freely falling graviton propagating in AdS in the context of the D1D5 CFT, where we introduce an interaction by turning on a deformation operator. We start with one left and right moving boson in the CFT. After applying two deformation operators, the initial bosons split into three left moving and three right moving bosons. We compute the amplitude for various energies and extrapolate the result to the large energy region. At early times, the amplitude is linear in time. This corresponds to an infalling graviton becoming redshifted in AdS. At late times, the amplitude is periodic, which agrees with the fact that a freely falling graviton will not be thermalized., Comment: 22 pages, 3 figures

Power device reliability is a major concern during operation under extreme environments, as doing so reduces the operational lifetime of any power system or sensing infrastructure. Due to a potential for system failure, devices must be experimentally validated before implementation, which is expensive and time-consuming. In this paper, we have utilized machine learning algorithms to predict device reliability, significantly reducing the need for conducting experiments. To train the models, we have tested 224 power devices from 10 different manufacturers. First, we describe a method to process the data for modeling purposes. Based on the in-house testing data, we implemented various ML models and observed that computational models such as Gradient Boosting and LSTM encoder-decoder networks can predict power device failure with high accuracy., Comment: 11 pages, 8 figures. Submitted to IEEE Transactions on Device and Materials Reliability

We study the dynamics of the three-dimensional quadratic diffeomorphism using a concept first introduced thirty years ago for the Frenkel-Kontorova model of condensed matter physics: the anti-integrable (AI) limit. At the traditional AI limit, orbits of a map degenerate to sequences of symbols and the dynamics is reduced to the shift operator, a pure form of chaos. Under nondegeneracy conditions, a contraction mapping argument can show that infinitely many AI states continue to orbits of the deterministic map. For the 3D quadratic map, the AI limit that we study is a quadratic correspondence whose branches, a pair of one-dimensional maps, introduce symbolic dynamics on two symbols. The AI states, however, are nontrivial orbits of this correspondence. The character of these orbits depends on whether the quadratic takes the form of an ellipse, a hyperbola, or a pair of lines. Using contraction arguments, we find parameter domains for each case such that each symbol sequence corresponds to a unique AI state. In some parameter domains, sufficient conditions are then found for each such AI state to continue away from the limit to become an orbit of the original 3D map. Numerical continuation methods extend these results, allowing computation of bifurcations and obtaining orbits with horseshoe-like structures and intriguing self-similarity. We conjecture that pairs of periodic orbits in saddle-node or period doubling bifurcations have symbol sequences that differ in exactly one position., Comment: 28 pages, 34 figures

Tidal effects in capped geometries computed in previous literature display no dynamics along internal (toroidal) directions. However, the dual CFT picture suggests otherwise. To resolve this tension, we consider a set of infalling null geodesics in a family of black hole microstate geometries with a smooth cap at the bottom of a long BTZ-like throat. Using the Penrose limit, we show that a string following one of these geodesics feels tidal stresses along all spatial directions, including internal toroidal directions. We find that the tidal effects along the internal directions are of the same order of magnitude as those along other, non-internal, directions. Furthermore, these tidal effects oscillate as a function of the distance from the cap -- as a string falls down the throat it alternately experiences compression and stretching. We explain some physical properties of this oscillation and comment on the dual CFT interpretation., Comment: 18 pages, 5 figures. V2: Minor modifications, accepted in JHEP

A ubiquitous challenge in design space exploration or uncertainty quantification of complex engineering problems is the minimization of computational cost. A useful tool to ease the burden of solving such systems is model reduction. This work considers a stochastic model reduction method (SMR), in the context of polynomial chaos (PC) expansions, where low-fidelity (LF) samples are leveraged to form a stochastic reduced basis. The reduced basis enables the construction of a bi-fidelity (BF) estimate of a quantity of interest from a small number of high-fidelity (HF) samples. A successful BF estimate approximates the quantity of interest with accuracy comparable to the HF model and computational expense close to the LF model. We develop new error bounds for the SMR approach and present a procedure to practically utilize these bounds in order to assess the appropriateness of a given pair of LF and HF models for BF estimation. The effectiveness of the SMR approach, and the utility of the error bound are presented in three numerical examples., Comment: 22 pages, 21 figures

Synthetic Aperture Sonar (SAS) surveys produce imagery with large regions of transition between seabed types. Due to these regions, it is difficult to label and segment the imagery and, furthermore, challenging to score the image segmentations appropriately. While there are many approaches to quantify performance in standard crisp segmentation schemes, drawing hard boundaries in remote sensing imagery where gradients and regions of uncertainty exist is inappropriate. These cases warrant weak labels and an associated appropriate scoring approach. In this paper, a labeling approach and associated modified version of the Rand index for weakly-labeled data is introduced to address these issues. Results are evaluated with the new index and compared to traditional segmentation evaluation methods. Experimental results on a SAS data set containing must-link and cannot-link labels show that our Weakly-Labeled Rand index scores segmentations appropriately in reference to qualitative performance and is more suitable than traditional quantitative metrics for scoring weakly-labeled data.

GenoML is a Python package automating machine learning workflows for genomics (genetics and multi-omics) with an open science philosophy. Genomics data require significant domain expertise to clean, pre-process, harmonize and perform quality control of the data. Furthermore, tuning, validation, and interpretation involve taking into account the biology and possibly the limitations of the underlying data collection, protocols, and technology. GenoML's mission is to bring machine learning for genomics and clinical data to non-experts by developing an easy-to-use tool that automates the full development, evaluation, and deployment process. Emphasis is put on open science to make workflows easily accessible, replicable, and transferable within the scientific community. Source code and documentation is available at https://genoml.com.

This paper unapologetically reflects on the critical role that Black feminism can and should play in abolishing algorithmic oppression. Positioning algorithmic oppression in the broader field of feminist science and technology studies, I draw upon feminist philosophical critiques of science and technology and discuss histories and continuities of scientific oppression against historically marginalized people. Moreover, I examine the concepts of invisibility and hypervisibility in oppressive technologies a l\'a the canonical double bind. Furthermore, I discuss what it means to call for diversity as a solution to algorithmic violence, and I critique dialectics of the fairness, accountability, and transparency community. I end by inviting you to envision and imagine the struggle to abolish algorithmic oppression by abolishing oppressive systems and shifting algorithmic development practices, including engaging our communities in scientific processes, centering marginalized communities in design, and consensual data and algorithmic practices., Comment: 12 pages, accepted to ACM Conference on Fairness, Accountability, and Transparency 2021

The multilevel Monte Carlo (MLMC) method has proven to be an effective variance-reduction statistical method for Uncertainty quantification in PDE models. It combines approximations at different levels of accuracy using a hierarchy of meshes in a similar way as multigrid. The generation of body-fitted mesh hierarchies is only possible for simple geometries. On top of that, MLMC for random domains involves the generation of a mesh for every sample. Instead, here we consider the use of embedded methods which make use of simple background meshes of an artificial domain (a bounding-box) for which it is easy to define a mesh hierarchy, thus eliminating the need of body-fitted unstructured meshes, but can produce ill-conditioned discrete problems. To avoid this complication, we consider the recent aggregated finite element method (AgFEM). In particular, we design an embedded MLMC framework for (geometrically and topologically) random domains implicitly defined through a random level-set function, which makes use of a set of hierarchical background meshes and the AgFEM. Performance predictions from existing theory are verified statistically in three numerical experiments, namely the solution of the Poisson equation on a circular domain of random radius, the solution of the Poisson equation on a topologically identical but more complex domain, and the solution of a heat-transfer problem in a domain that has geometric and topological uncertainties. Finally, the use of AgFE is statistically demonstrated to be crucial for complex and uncertain geometries in terms of robustness and computational cost.

It is generally agreed that black hole formation in gravity corresponds to thermalization in the dual CFT. It is sometimes argued that if the CFT evolution shows evidence of large redshift in gravity, then we have seen black hole formation in the CFT. We argue that this is not the case: a clock falling towards the horizon increases its redshift but remains intact as a clock; thus it is not `thermalized'. Instead, thermalization should correspond to a new phase after the phase of large redshift, where the infalling object turns into fuzzballs on reaching within planck distance of the horizon. We compute simple examples of the scattering vertex in the D1D5 CFT which, after many iterations, would lead to thermalization. An initial state made of two left-moving and two right-moving excitations corresponds, in gravity, to two gravitons heading towards each other. The thermalization vertex in the CFT breaks these excitations into multiple excitations on the left and right sides; we compute the amplitudes for several of these processes. We find secular terms that grow as $t^2$ instead of oscillating with $t$; we conjecture that this may be a feature of processes leading to thermalization., Comment: 25 pages

We consider quantum quench in large-N singlet sector quantum mechanics of a single hermitian matrix in the double scaling limit. The time dependent parameter is the self-coupling of the matrix. We find exact classical solutions of the collective field theory of eigenvalue density with abrupt smooth quench profiles which asymptote to constant couplings at early and late times, and the system is initially in its ground state. With adiabatic initial conditions we find that adiabaticity is always broken regardless of the quench speed. In a class of quench profiles the saddle point solution for the collective field diverges at a finite time, and a further time evolution becomes ambiguous. However the underlying matrix model expressed in terms of fermions predict a smooth time evolution across this point. By studying fluctuations around the saddle point solution we interpret the emergent space-times. They generically have spacelike boundaries where the couplings of the fluctuations diverge. Only for very finely tuned quench profiles, the space-time is normal., Comment: 44 pages, 15 figures; version accepted for publication in JHEP

The strongly coupled dynamics of black hole formation in bulk AdS is conjectured to be dual to the thermalization of a weakly interacting CFT on the boundary for low $N$ which, for $N\to\infty$, becomes strongly coupled. We search for this thermalization effect by utilizing the D1D5 CFT to compute effective string interactions for $N=2$. This is done by turning on a marginal deformation of the theory which twists together or untwists effective strings. For a system to thermalize, the initial state, which is far from thermal, must redistribute it's energy via interactions until a thermal state is achieved. In our case, we consider excited states of the effective strings. We compute splitting amplitudes for 1) one excitation going to three excitations and 2) two excitations going to four excitations using two insertions of the deformation. Scenario 1) corresponds to a single particle moving in AdS. Scenario 2) corresponds to two particles moving and colliding in AdS. We find that the `1 to 3' amplitude has terms which oscillate with time, $t$, where $t$ is the duration of the two deformations. We find that the `2 to 4' amplitude has similar oscillatory terms as well as secular terms which grow like $t^2$. For this case the growth implies that for large $t$ the excitations in the initial state, which carry a given energy, prefer to redistribute themselves amongst lower energy modes in the final state. This is a key feature of thermalization. Albeit in a simplified setting, we therefore argue that we have identified the thermalization vertex in the D1D5 CFT, which after repeated applications, should lead to thermalization. This ultimately maps to two particles colliding and forming a black hole in AdS, which in our case, is a fuzzball., Comment: PhD Dissertation, 443 pages, 18 figures

We investigate a class of exactly solvable quantum quench protocols with a finite quench rate in systems of one dimensional non-relativistic fermions in external harmonic oscillator or inverted harmonic oscillator potentials, with time dependent masses and frequencies. These hamiltonians arise, respectively, in harmonic traps, and the $c=1$ Matrix Model description of two dimensional string theory with time dependent string coupling. We show how the dynamics is determined by a single function of time which satisfies a generalized Ermakov-Pinney equation. The quench protocols we consider asymptote to constant masses and frequencies at early times, and cross or approach a gapless potential. In a right side up harmonic oscillator potential we determine the scaling behavior of the one point function and the entanglement entropy of a subregion by obtaining analytic approximations to the exact answers. The results are consistent with Kibble-Zurek scaling for slow quenches and with perturbation calculations for fast quenches. For cis-critical quench protocols the entanglement entropy oscillates at late times around its initial value. For end-critical protocols the entanglement entropy monotonically goes to zero inversely with time, reflecting the spread of fermions over the entire line. For the inverted harmonic oscillator potential, the dual collective field description is a scalar field in a time dependent metric and dilaton background., Comment: 44 pages, 10 figures

Topology optimization under uncertainty (TOuU) often defines objectives and constraints by statistical moments of geometric and physical quantities of interest. Most traditional TOuU methods use gradient-based optimization algorithms and rely on accurate estimates of the statistical moments and their gradients, e.g., via adjoint calculations. When the number of uncertain inputs is large or the quantities of interest exhibit large variability, a large number of adjoint (and/or forward) solves may be required to ensure the accuracy of these gradients. The optimization procedure itself often requires a large number of iterations, which may render TOuU computationally expensive, if not infeasible. To tackle this difficulty, we here propose an optimization approach that generates a stochastic approximation of the objective, constraints, and their gradients via a small number of adjoint (and/or forward) solves, per iteration. A statistically independent (stochastic) approximation of these quantities is generated at each optimization iteration. The total cost of this approach is only a small factor larger than that of the corresponding deterministic TO problem. We incorporate the stochastic approximation of objective, constraints and their design sensitivities into two classes of optimization algorithms. First, we investigate the stochastic gradient descent (SGD) method and a number of its variants, which have been successfully applied to large-scale optimization problems for machine learning. Second, we study the use of the proposed stochastic approximation approach within conventional nonlinear programming methods, focusing on the Globally Convergent Method of Moving Asymptotes (GCMMA). The performance of these algorithms is investigated with structural design optimization problems utilizing a Solid Isotropic Material with Penalization (SIMP), as well as an explicit level set method., Comment: 44 pages, 22 figures

Central configurations give rise to self-similar solutions to the Newtonian $N$-body problem, and play important roles in understanding its complicated dynamics. Even the simple question of whether or not there are finitely many planar central configurations for $N$ positive masses remains unsolved in most cases. Considering central configurations as critical points of a function $f$, we explicity compute the eigenvalues of the Hessian of $f$ for all $N$ for the point vortex potential for the regular polygon with equal masses. For homogeneous potentials including the Newtonian case we compute bounds on the eigenvalues for the regular polygon with equal masses, and give estimates on where bifurcations occur. These eigenvalue computations imply results on the Morse indices of $f$ for the regular polygon. Explicit formulae for the eigenvalues of the Hessian are also given for all central configurations of the equal mass 4-body problem with a homogeneous potential. Classic results on collinear central configurations are also generalized to the homogeneous potential case. Numerical results, conjectures, and suggestions for future work in the context of a homogeneous potential are given.

Over half a million individuals are diagnosed with head and neck cancer each year worldwide. Radiotherapy is an important curative treatment for this disease, but it requires manual time consuming delineation of radio-sensitive organs at risk (OARs). This planning process can delay treatment, while also introducing inter-operator variability with resulting downstream radiation dose differences. While auto-segmentation algorithms offer a potentially time-saving solution, the challenges in defining, quantifying and achieving expert performance remain. Adopting a deep learning approach, we demonstrate a 3D U-Net architecture that achieves expert-level performance in delineating 21 distinct head and neck OARs commonly segmented in clinical practice. The model was trained on a dataset of 663 deidentified computed tomography (CT) scans acquired in routine clinical practice and with both segmentations taken from clinical practice and segmentations created by experienced radiographers as part of this research, all in accordance with consensus OAR definitions. We demonstrate the model's clinical applicability by assessing its performance on a test set of 21 CT scans from clinical practice, each with the 21 OARs segmented by two independent experts. We also introduce surface Dice similarity coefficient (surface DSC), a new metric for the comparison of organ delineation, to quantify deviation between OAR surface contours rather than volumes, better reflecting the clinical task of correcting errors in the automated organ segmentations. The model's generalisability is then demonstrated on two distinct open source datasets, reflecting different centres and countries to model training. With appropriate validation studies and regulatory approvals, this system could improve the efficiency, consistency, and safety of radiotherapy pathways.

Female role models reduce the impact on women of stereotype threat, i.e., of being at risk of conforming to a negative stereotype about one's social, gender, or racial group. This can lead women scientists to underperform or to leave their scientific career because of negative stereotypes such as, not being as talented or as interested in science as men. Sadly, history rarely provides role models for women scientists; instead, it often renders these women invisible. In response to this situation, we present a selection of twelve outstanding women who helped to develop nuclear astrophysics., Comment: 5 pages; to appear in Springer Proceedings in Physics (Proc. of Intl. Conf. "Nuclei in the Cosmos XV", LNGS Assergi, Italy, June 2018)

Although studying outflows in the host galaxies of AGN have become the forefront of extra-galactic astronomy in recent years, estimating the energy associated with these outflows have been a major challenge. Determination of the energy associated with an outflow often involves an assumption of uniform density in the NLR, which span a wide range in literature leading to large systematic uncertainties in energy estimation. In this paper, we present electron density maps for a sample of outflowing and non-outflowing Seyfert galaxies at z<0.02 drawn from the S7 survey and understand the origin and values of the observed density structures to reduce the systematic uncertainties in outflow energy estimation. We use the ratio of the [SII]6716,6731 emission lines to derive spatially resolved electron densities (<50-2000 cm$^{-3}$). Using optical IFU observations, we are able to measure densities across the central 2-5 kpc of the selected AGN host galaxies. We compare the density maps with the positions of the HII regions derived from the narrow H$\alpha$ component, ionization maps from [OIII], and spatially resolved BPT diagrams, to infer the origin of the observed density structures. We also use the electron density maps to construct density profiles as a function of distance from the central AGN. We find a spatial correlation between the sites of high star formation and high electron density for targets without an active ionized outflow. The non-outflowing targets also show an exponential drop in the electron density as a function of distance from the center, with a mean exponential index of ~0.15. The correlation between the star forming sites and electron density ceases for targets with an outflow. The density within the outflowing medium is not uniform and shows both low and high density sites, most likely due to the presence of shocks and highly turbulent medium., Comment: Accepted for publication in A&A, 16 pages, 10 figures, 2 tables

We consider states of the D1-D5 CFT where only the left-moving sector is excited. As we deform away from the orbifold point, some of these states will remain BPS while others can `lift'. We compute this lifting for a particular family of D1-D5-P states, at second order in the deformation off the orbifold point. We note that the maximally twisted sector of the CFT is special: the covering surface appearing in the correlator can only be genus one while for other sectors there is always a genus zero contribution. We use the results to argue that fuzzball configurations should be studied for the full class including both extremal and near-extremal states; many extremal configurations may be best seen as special limits of near extremal configurations., Comment: 51 pages, 6 figures

We present the ~800 star formation rate maps for the SAMI Galaxy Survey based on H{\alpha} emission maps, corrected for dust attenuation via the Balmer decrement, that are included in the SAMI Public Data Release 1. We mask out spaxels contaminated by non-stellar emission using the [O III]/H{\beta}, [N II]/H{\alpha}, [S II]/H{\alpha}, and [O I]/H{\alpha} line ratios. Using these maps, we examine the global and resolved star-forming main sequences of SAMI galaxies as a function of morphology, environmental density, and stellar mass. Galaxies further below the star-forming main sequence are more likely to have flatter star formation profiles. Early-type galaxies split into two populations with similar stellar masses and central stellar mass surface densities. The main sequence population has centrally-concentrated star formation similar to late-type galaxies, while galaxies >3{\sigma} below the main sequence show significantly reduced star formation most strikingly in the nuclear regions. The split populations support a two-step quenching mechanism, wherein halo mass first cuts off the gas supply and remaining gas continues to form stars until the local stellar mass surface density can stabilize the reduced remaining fuel against further star formation. Across all morphologies, galaxies in denser environments show a decreased specific star formation rate from the outside in, supporting an environmental cause for quenching, such as ram-pressure stripping or galaxy interactions., Comment: accepted to MNRAS

In the field of uncertainty quantification, sparse polynomial chaos (PC) expansions are commonly used by researchers for a variety of purposes, such as surrogate modeling. Ideas from compressed sensing may be employed to exploit this sparsity in order to reduce computational costs. A class of greedy compressed sensing algorithms use least squares minimization to approximate PC coefficients. This least squares problem lends itself to the theory of optimal design of experiments (ODE). Our work focuses on selecting an experimental design that improves the accuracy of sparse PC approximations for a fixed computational budget. We propose DSP, a novel sequential design, greedy algorithm for sparse PC approximation. The algorithm sequentially augments an experimental design according to a set of the basis polynomials deemed important by the magnitude of their coefficients, at each iteration. Our algorithm incorporates topics from ODE to estimate the PC coefficients. A variety of numerical simulations are performed on three physical models and manufactured sparse PC expansions to provide a comparative study between our proposed algorithm and other non-adaptive methods. Further, we examine the importance of sampling by comparing different strategies in terms of their ability to generate a candidate pool from which an optimal experimental design is chosen. It is demonstrated that the most accurate PC coefficient approximations, with the least variability, are produced with our design-adaptive greedy algorithm and the use of a studied importance sampling strategy. We provide theoretical and numerical results which show that using an optimal sampling strategy for the candidate pool is key, both in terms of accuracy in the approximation, but also in terms of constructing an optimal design.

In the fuzzball paradigm the information paradox is resolved because the black hole is replaced by an object with no horizon. One may therefore ask if observations can distinguish a traditional hole from a fuzzball. We find: (a) It is very difficult to reflect quanta off the surface of a fuzzball, mainly because geodesics starting near the horizon radius cannot escape to infinity unless their starting direction is very close to radial. (b) If infalling particles interact with the emerging radiation before they are engulfed by the horizon, then we say that we have a `firewall behavior'. We consider several types of interactions, but find no evidence for firewall behavior in any theory that obeys causality. (c) Photons with wavelengths {\it larger} than the black hole radius can be scattered off the emerging radiation, but a very small fraction of the backscattered photons will be able to escape back to infinity., Comment: 52 pages, 4 figures

For practical model-based demands, such as design space exploration and uncertainty quantification (UQ), a high-fidelity model that produces accurate outputs often has high computational cost, while a low-fidelity model with less accurate outputs has low computational cost. It is often possible to construct a bi-fidelity model having accuracy comparable with the high-fidelity model and computational cost comparable with the low-fidelity model. This work presents the construction and analysis of a non-intrusive (i.e., sample-based) bi-fidelity model that relies on the low-rank structure of the map between model parameters/uncertain inputs and the solution of interest, if exists. Specifically, we derive a novel, pragmatic estimate for the error committed by this bi-fidelity model. We show that this error bound can be used to determine if a given pair of low- and high-fidelity models will lead to an accurate bi-fidelity approximation. The cost of this error bound is relatively small and depends on the solution rank. The value of this error estimate is demonstrated using two example problems in the context of UQ, involving linear and non-linear partial differential equations.

We present the second and final data release of the Siding Spring Southern Seyfert Spectroscopic Snapshot Survey (S7). Data are presented for 63 new galaxies not included in the first data release, and we provide 2D emission-line fitting products for the full S7 sample of 131 galaxies. The S7 uses the WiFeS instrument on the ANU 2.3m telescope to obtain spectra with a spectral resolution of R = 7000 in the red (540 - 700 nm) and R = 3000 in the blue (350 - 570 nm), over an integral field of 25x38 arcsec^2 with 1x1 arcsec^2 spatial pixels. The S7 contains both the largest sample of active galaxies and the highest spectral resolution of any comparable integral field survey to date. The emission-line fitting products include line fluxes, velocities and velocity dispersions across the WiFeS field of view, and an artificial neural network has been used to determine the optimal number of Gaussian kinematic components for emission-lines in each spaxel. Broad Balmer lines are subtracted from the spectra of nuclear spatial pixels in Seyfert 1 galaxies before fitting the narrow lines. We bin nuclear spectra and measure reddening-corrected nuclear fluxes of strong narrow lines for each galaxy. The nuclear spectra are classified on optical diagnostic diagrams, where the strength of the coronal line [FeVII]6087 is shown to be correlated with [OIII]/Hbeta. Maps revealing gas excitation and kinematics are included for the entire sample, and we provide notes on the newly observed objects., Comment: 37 pages, 18 figures, 4 tables. Accepted for ApJSS. Data release available at: https://datacentral.aao.gov.au/asvo/surveys/s7/

We present the first major release of data from the SAMI Galaxy Survey. This data release focuses on the emission-line physics of galaxies. Data Release One includes data for 772 galaxies, about 20% of the full survey. Galaxies included have the redshift range 0.004 < z < 0.092, a large mass range (7.6 < log(Mstellar/M$_\odot$) < 11.6), and star-formation rates of 10^-4 to 10^1\ M$_\odot$/yr. For each galaxy, we include two spectral cubes and a set of spatially resolved 2D maps: single- and multi-component emission-line fits (with dust extinction corrections for strong lines), local dust extinction and star-formation rate. Calibration of the fibre throughputs, fluxes and differential-atmospheric-refraction has been improved over the Early Data Release. The data have average spatial resolution of 2.16 arcsec (FWHM) over the 15~arcsec diameter field of view and spectral (kinematic) resolution R=4263 (sigma=30km/s) around Halpha. The relative flux calibration is better than 5\% and absolute flux calibration better than $\pm0.22$~mag, with the latter estimate limited by galaxy photometry. The data are presented online through the Australian Astronomical Observatory's Data Central., Comment: Submitted to MNRAS. SAMI DR1 data products available from http://datacentral.aao.gov.au/asvo/surveys/sami/

The most rapidly evolving regions of galaxies often display complex optical spectra with emission lines excited by massive stars, shocks and accretion onto supermassive black holes. Standard calibrations (such as for the star formation rate) cannot be applied to such mixed spectra. In this paper we isolate the contributions of star formation, shock excitation and active galactic nucleus (AGN) activity to the emission line luminosities of individual spatially resolved regions across the central 3 $\times$ 3 kpc$^2$ region of the active barred spiral galaxy NGC$\sim$613. The star formation rate and AGN luminosity calculated from the decomposed emission line maps are in close agreement with independent estimates from data at other wavelengths. The star formation component traces the B-band stellar continuum emission, and the AGN component forms an ionization cone which is aligned with the nuclear radio jet. The optical line emission associated with shock excitation is cospatial with strong $H_2$ and [Fe II] emission and with regions of high ionized gas velocity dispersion ($\sigma > 100$ km s$^{-1}$). The shock component also traces the outer boundary of the AGN ionization cone and may therefore be produced by outflowing material interacting with the surrounding interstellar medium. Our decomposition method makes it possible to determine the properties of star formation, shock excitation and AGN activity from optical spectra, without contamination from other ionization mechanisms., Comment: 16 pages, 12 figures. Accepted for publication in MNRAS

Foldable paper constructions known as flexagons have been studied since 1939. In this paper we review the construction of hexagonal flexagons (hexaflexagons) and compute the number of distinct hexaflexagons with n faces.