title: Kinematics of Simulated Galaxies I: Connecting Dynamical and Morphological Properties of Early-Type Galaxies at Different Redshifts

abstract: State-of-the-art integral field surveys like $\mathrm{ATLAS^{3D}}$, SLUGGS, CALIFA, SAMI, and MaNGA provide large data sets of kinematical observations of early-type galaxies (ETGs), yielding constraints on the formation of ETGs. Using the cosmological hydrodynamical \textit{Magneticum Pathfinder} simulations, we investigate the paradigm of fast and slow rotating ETGs in a fully cosmological context. We show that the ETGs within the \textit{Magneticum} simulation are in remarkable agreement with the observations, revealing fast and slow rotators quantified by the angular momentum proxy $\lambda_{\mathrm{R}}$ and the flattening $\epsilon$ with the observed prevalence. Taking full advantage of the three-dimensional data, we demonstrate that the dichotomy between fast and slow rotating galaxies gets enhanced, showing an upper and lower population separated by an underpopulated region in the edge-on $\lambda_{\mathrm{R}}$-$\epsilon$ plane. Following the evolution of the $\lambda_{\mathrm{R}}$-$\epsilon$ plane through cosmic time, we find that, while the upper population is already in place at $z=2$, the lower population gets statistically significant below $z=1$ with a gradual increase. At least $50\%$ of the galaxies transition from fast to slow rotators on a short timescale, in most cases associated to a significant merger event. Furthermore, we connect the $M_{*}$-$j_{*}$ plane, quantified by the $b$-value, with the $\lambda_{\mathrm{R}}$-$\epsilon$ plane, revealing a strong correlation between the position of a galaxy in the $\lambda_{\mathrm{R}}$-$\epsilon$ plane and the $b$-value. Going one step further, we classify our sample based on features in their velocity map, finding all five common kinematic groups, also including the recently observed group of prolate rotators, populating distinct regions in the $\lambda_{\mathrm{R}}$-$b$ plane.


title: Measuring the value of the Hubble constant "\`a la Refsdal"

abstract: Realizing Refsdal's original idea from 1964, we present estimates of the Hubble constant that are complementary to and potentially competitive with those of other cosmological probes. We use the observed positions of 89 multiple images, with extensive spectroscopic information, from 28 background sources and the measured time delays between the images S1-S4 and SX of supernova "Refsdal" ($z = 1.489$), which were obtained thanks to Hubble Space Telescope (HST) deep imaging and Multi Unit Spectroscopic Explorer (MUSE) data. We extend the strong lensing modeling of the Hubble Frontier Fields (HFF) galaxy cluster MACS J1149.5$+$2223 ($z = 0.542$), published by Grillo et al. (2016), and explore different $\Lambda$CDM models. Taking advantage of the lensing information associated to the presence of very close pairs of multiple images at various redshifts and to the extended surface brightness distribution of the SN Refsdal host, we can reconstruct the total mass density profile of the cluster very precisely. The combined dependence of the multiple image positions and time delays on the cosmological parameters allows us to infer the values of $H_{0}$ and $\Omega_{\rm m}$ with relative (1$\sigma$) statistical errors of, respectively, 6% (7%) and 31% (26%) in flat (general) cosmological models, assuming a conservative 3% uncertainty on the final time delay of image SX and, remarkably, no priors from other cosmological experiments. Our best estimate of $H_{0}$, based on the model described in this work, will be presented when the final time-delay measurement becomes available. Our results show that it is possible to utilize time delays in lens galaxy clusters as an important alternative tool for measuring the expansion rate and the geometry of the Universe.


title: Positive gravitational subsystem energies from CFT cone relative entropies

abstract: The positivity of relative entropy for spatial subsystems in a holographic CFT implies the positivity of certain quantities in the dual gravitational theory. In this note, we consider CFT subsystems whose boundaries lie on the lightcone of a point $p$. We show that the positive gravitational quantity which corresponds to the relative entropy for such a subsystem $A$ is a novel notion of energy associated with a gravitational subsystem bounded by the minimal area extremal surface $\tilde{A}$ associated with $A$ and by the AdS boundary region $\hat{A}$ corresponding to the part of the lightcone from $p$ bounded by $\partial A$. This generalizes the results of arXiv:1605.01075 for ball-shaped regions by making use of the recent results in arXiv:1703.10656 for the vacuum modular Hamiltonian of regions bounded on lightcones. As part of our analysis, we give an analytic expression for the extremal surface in pure AdS associated with any such region $A$. We note that its form immediately implies the Markov property of the CFT vacuum (saturation of strong subadditivity) for regions bounded on the same lightcone. This gives a holographic proof of the result proven for general CFTs in arXiv:1703.10656. A similar holographic proof shows the Markov property for regions bounded on a lightsheet for non-conformal holographic theories defined by relevant perturbations of a CFT.


title: A Census of the LyC Photons that Form the UV Background During Reionization

abstract: We present a new, on-the-fly photon flux and absorption tracer algorithm designed to directly measure the contribution of different source populations to the metagalactic UV background and to the ionisation fraction of gas in the Universe. We use a suite of multifrequency radiation hydrodynamics simulations that are carefully calibrated to reproduce a realistic reionization history and galaxy properties at $z \ge 6$, to disentangle the contribution of photons emitted by different mass haloes and by stars with different metallicities and ages to the UV background during reionization. While at very early cosmic times low mass, metal poor haloes provide most of the LyC photons, their contribution decreases steadily with time. At $z = 6$ it is the photons emitted by massive systems (${\rm M_{halo}}/{\rm M_\odot} > 10^{10} \, {\rm h ^{-1}}$) and by the metal enriched stars ($10^{-3} < Z/Z_{\rm \odot} < 10^{-1.5}$) that provide the largest contribution to the ionising UV background. We demonstrate that there are large variations in the escape fraction depending on the source, with the escape fraction being highest ($\sim 45-60\%$) for photons emitted by the oldest stars that penetrate into the IGM via low opacity channels carved by the ionising photons and supernova from younger stars. Before HII regions begin to overlap, the photoionisation rate strongly fluctuates between different, isolated HII bubbles, depending on the embedded ionising source, which we suggest may result in spatial variations in the properties of dwarf galaxies.


title: Resilience of scrambling measurements

abstract: Most experimental protocols for measuring scrambling require time evolution with a Hamiltonian and with the Hamiltonian's negative counterpart (backwards time evolution). Engineering controllable quantum many-body systems for which such forward and backward evolution is possible is a significant experimental challenge. Furthermore, if the system of interest is quantum-chaotic, one might worry that any small errors in the time reversal will be rapidly amplified, obscuring the physics of scrambling. This paper undermines this expectation: We exhibit a renormalization protocol that extracts nearly ideal out-of-time-ordered-correlator measurements from imperfect experimental measurements. We analytically and numerically demonstrate the protocol's effectiveness, up to the scrambling time, in a variety of models and for sizable imperfections. The scheme extends to errors from decoherence by an environment.

categories: quant-ph cond-mat.dis-nn cond-mat.stat-mech hep-th


title: The First Detection of Neutral Hydrogen in Emission in a Strong Spiral Lens

abstract: We report HI observations of eight spiral galaxies that are strongly lensing background sources. Our targets were selected from the Sloan WFC (Wide Field Camera) Edge-on Late-type Lens Survey (SWELLS) using the Arecibo, Karl G. Jansky Very Large Array, and Green Bank telescopes. We securely detect J1703+2451 at z=0.063 with a signal-to-noise of 6.7 and W50=79+/-13 km/s, obtaining the first detection of HI emission in a strong spiral lens. We measure a mass of M(HI)= 1.77+/-0.06(+0.35/-0.75) x 10^9 M_(sol) for this source. We find that this lens is a normal spiral, with observable properties that are fairly typical of spiral galaxies. For three other sources we did not secure a detection; however, we are able to place strong constraints on the HI masses of those galaxies. The observations for four of our sources were rendered unusable due to strong radio frequency interference.


title: PISCO: The Pmas/ppak Integral-field Supernova hosts COmpilation

abstract: We present the Pmas/ppak Integral-field Supernova hosts COmpilation (PISCO) which comprises Integral Field Spectroscopy (IFS) of 232 supernova (SN) host galaxies, that hosted 272 SNe, observed over several semesters with the 3.5m telescope at the Calar Alto Observatory (CAHA). PISCO is the largest collection of SN host galaxies observed with wide-field IFS, totaling 466,347 individual spectra covering a typical spatial resolution of $\sim$380 pc. While focused studies regarding specific SN Ia- related topics will be published elsewhere, this paper aims to present the properties of the SN environments with stellar population (SP) synthesis and the gas-phase ISM, providing additional results separating stripped-envelope SNe into their subtypes. With 11,270 HII regions detected in all galaxies, we present for the first time an HII region statistical analysis, that puts HII regions that have hosted SNe in context with all other SF clumps within their galaxies. SNe Ic are associated to more metal-rich, higher EW(H{\alpha}) and higher SF rate environments within their host galaxies than the mean of all HII regions detected within each host, on contrary SNe IIb occur at the most different environments compared to other CC SNe types. We find two clear components of young and old SP at SNe IIn locations. We find that SNe II fast-decliners (IIL) tend to explode at locations where {\Sigma}SFR is more intense. Finally, we outline how a future dedicated IFS survey of galaxies in parallel to an untargeted SN search would overcome the biases in current environmental studies.


title: Computational inverse method for constructing spaces of quantum models from wave functions

abstract: Traditional computational methods for studying quantum many-body systems are "forward methods," which take quantum models, i.e., Hamiltonians, as input and produce ground states as output. However, such forward methods often limit one's perspective to a small fraction of the space of possible Hamiltonians. We introduce an alternative computational "inverse method," the Eigenstate-to-Hamiltonian Construction (EHC), that allows us to better understand the vast space of quantum models describing strongly correlated systems. EHC takes as input a wave function $|\psi_T\rangle$ and produces as output Hamiltonians for which $|\psi_T\rangle$ is an eigenstate. This is accomplished by computing the quantum covariance matrix, a quantum mechanical generalization of a classical covariance matrix. EHC is widely applicable to a number of models and in this work we consider seven different examples. Using the EHC method, we construct a parent Hamiltonian with a new type of antiferromagnetic ground state, a parent Hamiltonian with two different targeted degenerate ground states, and large classes of parent Hamiltonians with the same ground states as well-known quantum models, such as the Majumdar-Ghosh model, the XX chain, the Heisenberg chain, the Kitaev chain, and a 2D BdG model.


title: Clockwork without supersymmetry

abstract: The clockwork mechanism would be completely spoiled by the presence of a cosmological constant in the bulk of its 5D construction, or by analogous terms in the discrete case. It is believed that supersymmetry is required to forbid the appearance of these fatal terms, thus apparently tying the fate of the clockwork to the presence of supersymmetry. In this letter we argue that a robust clockwork can be obtained also without supersymmetry, by providing a model for the clockwork/linear dilaton EFT where this protection is instead related to diffeomorphism invariance in a higher number of extra dimensions. We show that the clockwork mechanism is a quite general setup that emerges, rather minimally, from pure gravity in the presence of $D-5$ additional flat extra dimensions. The linear dilaton model is obtained asymptotically in the limit of infinitely many extra dimensions. We then study the finite-$D$ theory, which is a deformation of the linear dilaton one, but still a clockwork theory for $D$ large enough.


title: Deconstructing S-Duality

abstract: We use the technique of deconstruction to lift dualities from 2+1 to 3+1 dimensions. In this work we demonstrate the basic idea by deriving S-duality of maximally supersymmetric electromagnetism in 3+1 dimensions from mirror symmetry in 2+1. We also study the deconstruction of a non-supersymmetric duality in 3+1 dimensions using Abelian bosonization in 2+1 dimensions.


title: The enhancement of rapidly quenched galaxies in distant clusters at 0.5<z<1.0

abstract: We investigate the relationship between environment and galaxy evolution in the redshift range $0.5 < z < 1.0$. Galaxy overdensities are selected using a Friends-of-Friends algorithm, applied to deep photometric data in the Ultra-Deep Survey (UDS) field. A study of the resulting stellar mass functions reveals clear differences between cluster and field environments, with a strong excess of low-mass rapidly quenched galaxies in cluster environments compared to the field. Cluster environments also show a corresponding deficit of young, low-mass star-forming galaxies, which show a sharp radial decline towards cluster centres. By comparing mass functions and radial distributions, we conclude that young star-forming galaxies are rapidly quenched as they enter overdense environments, becoming post-starburst galaxies before joining the red sequence. Our results also point to the existence of two environmental quenching pathways operating in galaxy clusters, operating on different timescales. Fast quenching acts on galaxies with high specific star-formation rates, operating on timescales shorter than the cluster dynamical time ($ < 1$ Gyr). In contrast, slow quenching affects galaxies with moderate specific star-formation rates, regardless of their stellar mass, and acts on longer timescales ($\gtrsim 1$ Gyr). Of the cluster galaxies in the stellar mass range $9.0 < \log(M_{*}/M_{\odot}) < 10.5$ quenched during this epoch, we find that 73% were transformed through fast quenching, while the remaining 27% followed the slow quenching route.


title: Magnetized advective accretion flows: formation of magnetic barriers in Magnetically Arrested Discs

abstract: We discuss the importance of large scale strong magnetic field in the removal of angular momentum outward, as well as the possible origin of different kinds of magnetic barrier in advective, geometrically thick, sub-Keplerian accretion flows around black holes. The origin of this large scale strong magnetic field near the event horizon is due to the advection of the magnetic flux by the accreting gas from the environment, say, the interstellar medium or a companion star, because of flux freezing. In this simplest vertically averaged, 1.5-dimensional disc model, we choose the maximum upper limit of the magnetic field, which the disc around a black hole can sustain. In this so called magnetically arrested disc (MAD) model, the accreting gas either decelerates or faces the magnetic barrier near the event horizon by the accumulated magnetic field depending on the geometry. The magnetic barrier may knock the matter to infinity. We suggest that these types of flow are the building block to produce jets and outflows in the accreting system. We also find that in some cases, when matter is trying to go back to infinity after knocking the barrier, matter is prevented being escaped by the cumulative action of strong gravity and the magnetic tension, hence by another barrier. In this way, magnetic field can lock the matter in between these two barriers and it might be a possible explanation for the formation of episodic jet.


title: Search for the Neutron Decay n$\rightarrow$ X+$\gamma$ where X is a dark matter particle

abstract: In a recent paper submitted to Physical Review Letters, Fornal and Grinstein have suggested that the discrepancy between two different methods of neutron lifetime measurements, the beam and bottle methods can be explained by a previously unobserved dark matter decay mode, n$\rightarrow$ X+$\gamma$ where X is a dark matter particle. We have performed a search for this decay mode over the allowed range of energies of the monoenergetic gamma ray for X to be a dark matter particle. We exclude the possibility of a sufficiently strong branch to explain the lifetime discrepancy with greater than 4 sigma confidence.


title: The mass-size relation of LRGs from BOSS and DECaLS

abstract: We use the DECaLS DR3 survey photometry matched to the SDSS-III/BOSS DR12 spectroscopic catalog to investigate the morphology and stellar mass-size relation of luminous red galaxies (LRGs) within the CMASS and LOWZ galaxy samples in the redshift range $0.2<z<0.7$. The large majority of both samples is composed of early-type galaxies with De Vaucouleurs profiles, while only less than 20% are late-type exponentials. We calibrate DECaLS effective radii using the higher resolution CFHT/MegaCam observations and optimise the correction for each morphological type. By cross-matching the photometric properties of the early-type population with the Portsmouth stellar mass catalog, we are able to explore the high-mass end of the distribution using a large sample of 313,026 galaxies over 4380 deg$^{2}$. We find a clear correlation between the sizes and the stellar masses of these galaxies, which appears flatter than previous estimates at lower masses. The sizes of these early-type galaxies do not exhibit significant evolution within the BOSS redshift range, but a slightly declining redshift trend is found when these results are combined with $z\sim0.1$ SDSS measurements at the high-mass end. The synergy between BOSS and DECaLS has important applications in other fields, including galaxy clustering and weak lensing.


title: Mapping the core of the Tarantula Nebula with VLT-MUSE: I. Spectral and nebular content around R136

abstract: We introduce VLT-MUSE observations of the central 2$'\times2'$ (30$\times$30 pc) of the Tarantula Nebula in the Large Magellanic Cloud. The observations provide an unprecedented spectroscopic census of the massive stars and ionised gas in the vicinity of R136, the young, dense star cluster located in NGC 2070, at the heart of the richest star-forming region in the Local Group. Spectrophotometry and radial-velocity estimates of the nebular gas (superimposed on the stellar spectra) are provided for 2255 point sources extracted from the MUSE datacubes, and we present estimates of stellar radial velocities for 270 early-type stars (finding an average systemic velocity of 271$\pm$41 km/s). We present an extinction map constructed from the nebular Balmer lines, with electron densities and temperatures estimated from intensity ratios of the [SII], [NII], and [SIII] lines. The interstellar medium, as traced by H$\alpha$ and [NII] $\lambda$6583, provides new insights in regions where stars are probably forming. The gas kinematics are complex, but with a clear bi-modal, blue- and red-shifted distribution compared to the systemic velocity of the gas centred on R136. Interesting point-like sources are also seen in the eastern cavity, western shell, and around R136; these might be related to phenomena such as runaway stars, jets, formation of new stars, or the interaction of the gas with the population of Wolf--Rayet stars. Closer inspection of the core reveals red-shifted material surrounding the strongest X-ray sources, although we are unable to investigate the kinematics in detail as the stars are spatially unresolved in the MUSE data. Further papers in this series will discuss the detailed stellar content of NGC 2070 and its integrated stellar and nebular properties.


title: Scalar one-point functions and matrix product states of AdS/dCFT

abstract: We determine in a closed form all scalar one-point functions of the defect CFT dual to the D3-D5 probe brane system with k units of flux which amounts to calculating the overlap between a Bethe eigenstate of the integrable SO(6) spin chain and a certain matrix product state of bond dimension k. In particular, we show that the matrix product state is annihilated by all the parity odd charges of the spin chain which has recently been suggested as the criterion for such a state to correspond to an integrable initial state. Finally, we discuss the properties of the analogous matrix product state for the SO(5) symmetric D3-D7 probe brane set-up.


title: Cooper-Pair Spin Current in a Strontium Ruthenate Heterostructure

abstract: It has been recognized that the condensation of spin-triplet Cooper pairs requires not only the broken gauge symmetry but also the spin ordering as well. One consequence of this is the possibility of the Cooper-pair spin current analogous to the magnon spin current in magnetic insulators, the analogy also extending to the existence of the Gilbert damping of the collective spin-triplet dynamics. The recently fabricated heterostructure of the thin film of the itinerant ferromagnet SrRuO3 on the bulk Sr2RuO4, the best-known candidate material for the spin-triplet superconductor, offers a promising platform for generating such spin current. We will show how such heterostructure allows us to not only realize the long-range spin valve but also electrically drive the collective spin mode of the spin-triplet order parameter. Our proposal represents both a new realization of the spin superfluidity and a transport signature of the spin-triplet superconductivity.


title: Models of Tidally Induced Gas Filaments in the Magellanic Stream

abstract: The Magellanic Stream and Leading Arm of HI that stretches from the Large and Small Magellanic Clouds (LMC and SMC) and over 200 degrees of the Southern sky is thought to be formed from multiple encounters between the LMC and SMC. In this scenario, most of the gas in the Stream and Leading Arm is stripped from the SMC, yet recent observations have shown a bifurcation of the Trailing Arm that reveals LMC origins for some of the gas. Absorption measurements in the Stream also reveal an order of magnitude more gas than in current tidal models. We present hydrodynamical simulations of the multiple encounters between the LMC and SMC at their first pass around the Milky Way, assuming that the Clouds were more extended and gas rich in the past. Our models create filamentary structures of gas in the Trailing Stream from both the LMC and SMC. While the SMC trailing filament matches the observed Stream location, the LMC filament is offset. In addition, the total observed mass of the Stream in these models is underestimated of a factor of four when the ionized component is accounted for. Our results suggest that there should also be gas stripped from both the LMC and SMC in the Leading Arm, mirroring the bifurcation in the Trailing Stream. This prediction is consistent with recent measurements of spatial variation in chemical abundances in the Leading Arm, which show that gas from multiple sources is present, although the nature is still uncertain.


title: Quantum Fisher information for general spatial deformations of quantum emitters

abstract: We present a framework for the detection and estimation of deformations applied to a grid of sources. Our formalism uses the Hamiltonian formulation of the quantum Fisher information matrix (\textsc{qfim}) as the figure of merit to quantify the amount of information we have on the deformation matrix. Quantum metrology for grid deformations provides an ideal testbed to examine multi-parameter estimations for arbitrarily parameterised channel evolutions with generally non-commuting Hermitian generators. We generalise the local generator of translations for deformation parameters to multi-parameter estimations and use it to explore how well different deformations can be detected and corrected for. This approach holds for any deformation. We explore the application of our theory to the set of affine geometry maps. Both the configuration of the grid and the properties of the sources help to maximise the sensitivity of the \textsc{qfim} to changes in the deformation parameters. For the non-multiplicative Hamiltonian parameterisations resulting from grid rotations about any chosen axis, oscillatory dependence of the \textsc{qfi} surfaces for a specific interplay between mutual source separation distances and grid configurations.


title: Gauge-invariant observables, gravitational dressings, and holography in AdS

abstract: This paper explores construction of gauge (diffeomorphism)-invariant observables in anti de Sitter (AdS) space and the related question of how to find a "holographic map" providing a quantum equivalence to a boundary theory. Observables are constructed perturbatively to leading order in the gravitational coupling by gravitationally dressing local field theory operators in order to solve the gravitational constraints. Many such dressings are allowed and two are explicitly examined, corresponding to a gravitational line and to a Coulomb field; these also reveal an apparent role for more general boundary conditions than considered previously. The observables obey a nonlocal algebra, and we derive explicit expressions for the boundary generators of the SO(D-1,2) AdS isometries that act on them. We examine arguments that gravity {\it explains} holography through the role of such a boundary Hamiltonian. Our leading-order gravitational construction reveals some questions regarding how these arguments work, and indeed construction of such a holographic map appears to require solution of the non-perturbative generalization of the bulk constraint equations.


title: Evading no-hair theorems: hairy black holes in a Minkowski box

abstract: We find hairy black holes of Einstein-Maxwell theory with a complex scalar field that is confined inside a box in a Minkowski background. These regular hairy black holes are asymptotically flat and thus the presence of the box or mirror allows to evade well-known no-hair theorems. We also find the Israel surface stress tensor that the confining box must have to obey the energy conditions. In the zero horizon radius limit, these hairy black holes reduce to a regular asymptotically flat hairy soliton. We find our solutions using perturbation theory. At leading order, a hairy black hole can be seen as a Reissner-Nordstrom black hole placed on top of a hairy soliton with the same chemical potential (so that the system is in thermodynamic equilibrium). The hairy black holes merge with the Reissner-Nordstrom black hole family at the onset of the superradiant instability. When they co-exist, for a given energy and electric charge, hairy black holes have higher entropy than caged Reissner-Nordstrom black holes. Therefore, our hairy black holes are the natural candidates for the endpoint of charged superradiance in the Reissner-Nordstrom black hole mirror system.


title: Learning from Richer Human Guidance: Augmenting Comparison-Based Learning with Feature Queries

abstract: We focus on learning the desired objective function for a robot. Although trajectory demonstrations can be very informative of the desired objective, they can also be difficult for users to provide. Answers to comparison queries, asking which of two trajectories is preferable, are much easier for users, and have emerged as an effective alternative. Unfortunately, comparisons are far less informative. We propose that there is much richer information that users can easily provide and that robots ought to leverage. We focus on augmenting comparisons with feature queries, and introduce a unified formalism for treating all answers as observations about the true desired reward. We derive an active query selection algorithm, and test these queries in simulation and on real users. We find that richer, feature-augmented queries can extract more information faster, leading to robots that better match user preferences in their behavior.


title: Dimension-scalable recurrence threshold estimation

abstract: The appropriate selection of recurrence thresholds is a key problem in applications of recurrence quantification analysis (RQA) and related methods across disciplines. Here, we discuss the distribution of pairwise distances between state vectors in the studied system's state space reconstructed by means of time-delay embedding as the key characteristic that should guide the corresponding choice for obtaining an adequate resolution of a recurrence plot. Specifically, we present an empirical description of the distance distribution, focusing on characteristic changes of its shape with increasing embedding dimension. Based on our results, we recommend selecting the recurrence threshold adaptively according to a fixed quantile of this distribution. We highlight the advantages of this strategy over other previously suggested approaches by discussing the performance of selected RQA measures in detecting chaos--chaos transitions in some prototypical model system.


title: Scaling of polarization amplitude in quantum many-body systems in one dimension

abstract: Resta proposed a definition of the electric polarization in one-dimensional systems in terms of the ground-state expectation value of the large gauge transformation operator. Vanishing of the expectation value in the thermodynamic limit implies that the system is a conductor. We study Resta's polarization amplitude (expectation value) in the $S=1/2$ XXZ chain and its several generalizations, in the gapless conducting Tomonaga-Luttinger Liquid phase. We obtain an analytical expression in the lowest-order perturbation theory about the free fermion point (XY chain), and an exact result for the Haldane-Shastry model with long-range interactions. We also obtain numerical results, mostly using the exact diagonalization method. We find that the amplitude exhibits a power-law scaling in the system size (chain length) and vanishes in the thermodynamic limit. On the other hand, the exponent depends on the model even when the low-energy limit is described by the Tomonaga-Luttinger Liquid with the same Luttinger parameter. We find that a change in the exponent occurs when the Umklapp term(s) are eliminated, suggesting the importance of the Umklapp terms.


title: $hh+\text{jet}$ production at 100 TeV

abstract: Higgs pair production is a crucial phenomenological process in deciphering the nature of the TeV scale and the mechanism underlying electroweak symmetry breaking. At the Large Hadron Collider, this process is statistically limited. Pushing the energy frontier beyond the LHC's reach will create new opportunities to exploit the rich phenomenology at higher centre-of-mass energies and luminosities. In this work, we perform a comparative analysis of the $hh+\text{jet}$ channel at a future 100 TeV hadron collider. We focus on the $hh\to b\bar b b\bar b$ and $hh \to b\bar b \tau^+\tau^-$ channels and employ a range of analysis techniques to estimate the sensitivity potential that can be gained by including this jet-associated Higgs pair production to the list of sensitive collider processes in such an environment. In particular, we observe that $hh \to b\bar b \tau^+\tau^-$ in the boosted regime exhibits a large sensitivity to the Higgs boson self-coupling and the Higgs self-coupling could be constrained at the 8\% level in this channel alone.


title: On The Circular Altitude of Graphs

abstract: Peter Cameron introduced the concept of the circular altitude of graphs; a parameter which was shown by Bamberg et al. that provides a lower bound on the circular chromatic number. In this note, we investigate this parameter and show that the circular altitude of a graph is equal to the maximum of circular altitudes of its blocks. Also, we show that homomorphically equivalent graphs have the same circular altitudes. Finally, we prove that the circular altitude of the Cartesian product of two graphs is equal to the maximum of circular altitudes of its factors.


title: The KMOS Cluster Survey (KCS) II - The Effect of Environment on the Structural Properties of Massive Cluster Galaxies at Redshift $1.39 < z <1.61$

abstract: We present results on the structural properties of massive passive galaxies in three clusters at $1.39<z<1.61$ from the KMOS Cluster Survey. We measure light-weighted and mass-weighted sizes from optical and near-infrared Hubble Space Telescope imaging and spatially resolved stellar mass maps. The rest-frame $R$-band sizes of these galaxies are a factor of $\sim2-3$ smaller than their local counterparts. The slopes of the relation between the stellar mass and the light-weighted size are consistent with recent studies in clusters and the field. Their mass-weighted sizes are smaller than the rest frame $R$-band sizes, with an average mass-weighted to light-weighted size ratio that varies between $\sim0.45$ and $0.8$ among the clusters. We find that the median light-weighted size of the passive galaxies in the two more evolved clusters is $\sim24\%$ larger than for field galaxies, independent of the use of circularized effective radii or semi-major axes. These two clusters also show a smaller size ratio than the less evolved cluster, which we investigate using color gradients to probe the underlying $M_{*}/L_{H_{160}}$ gradients. The median color gradients are $\nabla{z-H} \sim-0.4$ mag dex$^{-1}$, twice the local value. Using stellar populations models, these gradients are best reproduced by a combination of age and metallicity gradients. Our results favor the minor merger scenario as the dominant process responsible for the observed galaxy properties and the environmental differences at this redshift. The environmental differences support that clusters experience accelerated structural evolution compared to the field, likely via an epoch of enhanced minor merger activity during cluster assembly.


title: Fast and accurate approximation of the full conditional for gamma shape parameters

abstract: The gamma distribution arises frequently in Bayesian models, but there is not an easy-to-use conjugate prior for the shape parameter of a gamma. This inconvenience is usually dealt with by using either Metropolis-Hastings moves, rejection sampling methods, or numerical integration. However, in models with a large number of shape parameters, these existing methods are slower or more complicated than one would like, making them burdensome in practice. It turns out that the full conditional distribution of the gamma shape parameter is well approximated by a gamma distribution, even for small sample sizes, when the prior on the shape parameter is also a gamma distribution. This article introduces a quick and easy algorithm for finding a gamma distribution that approximates the full conditional distribution of the shape parameter. We empirically demonstrate the speed and accuracy of the approximation across a wide range of conditions. If exactness is required, the approximation can be used as a proposal distribution for Metropolis-Hastings.


title: A Sterile Neutrino Origin for the Upward Directed Cosmic Ray Showers Detected by ANITA

abstract: The ANITA balloon experiment has recently observed several ~EeV cascade events at an angle below the horizon that renders any Standard Model (SM) interpretation unlikely as the Earth is significantly opaque to all SM particles at such energies. In this paper, we study a sterile neutrino interpretation of these events, calculating the angular acceptance of cascades and the relative sensitivities of several experiments to a cascade initiated by an EeV sterile neutrino. We find that ANITA is uniquely sensitive to this type of upward directed cascade signal over a wide portion of the sky and from the direction of the two observed events has a transient acceptance roughly equivalent to that of the IceCube experiment.


title: Sequence-dependent Three Interaction Site (TIS) Model for Single and Double-stranded DNA

abstract: We develop a robust coarse-grained model for single and double stranded DNA by representing each nucleotide by three interaction sites (TIS) located at the centers of mass of sugar, phosphate, and base. The resulting TIS model includes base-stacking, hydrogen bond, and electrostatic interactions as well as bond-stretching and bond angle potentials that account for the polymeric nature of DNA. The choices of force constants for stretching and the bending potentials were guided by a Boltzmann inversion procedure using a large representative set of DNA structures extracted from the Protein Data Bank. Some of the parameters in the stacking interactions were calculated using a learning procedure, which ensured that the experimentally measured melting temperatures of dimers are faithfully reproduced. Without any further adjustments, the calculations based on the TIS model reproduces the experimentally measured salt and sequence dependence of the size of single stranded DNA (ssDNA), as well as the persistence lengths of poly(dA) and poly(dT) chains. Interestingly, upon application of mechanical force the extension of poly(dA) exhibits a plateau, which we trace to the formation of stacked helical domains. In contrast, the force-extension curve (FEC) of poly(dT) is entropic in origin, and could be described by a standard polymer model. We also show that the persistence length of double stranded DNA is consistent with the prediction based on the worm-like chain. The persistence length, which decreases with increasing salt concentration, is in accord with the Odijk-Skolnick-Fixman theory intended for stiff polyelectrolyte chains near the rod limit. The range of applications, which did not require adjusting any parameter after the initial construction based solely on PDB structures and melting profiles of dimers, attests to the transferability and robustness of the TIS model for ssDNA and dsDNA.


title: The Inhomogeneous Reionization Times of Present-day Galaxies

abstract: Today's galaxies experienced cosmic reionization at different times in different locations. For the first time, reionization ($50\%$ ionized) redshifts, $z_R$, at the location of their progenitors are derived from new, fully-coupled radiation-hydrodynamics simulation of galaxy formation and reionization at $z > 6$, matched to N-body simulation to z = 0. Constrained initial conditions were chosen to form the well-known structures of the local universe, including the Local Group and Virgo, in a (91 Mpc)$^3$ volume large enough to model both global and local reionization. Reionization simulation CoDa I-AMR, by CPU-GPU code EMMA, used (2048)$^3$ particles and (2048)$^3$ initial cells, adaptively-refined, while N-body simulation CoDa I-DM2048, by Gadget2, used (2048)$^3$ particles, to find reionization times for all galaxies at z = 0 with masses $M(z=0)\ge 10^8 M_\odot$. Galaxies with $M(z=0) \gtrsim 10^{11} M_\odot$ reionized earlier than the universe as a whole, by up to $\sim$ 500 Myrs, with significant scatter. For Milky-Way-like galaxies, $z_R$ ranged from 8 to 15. Galaxies with $M(z=0) \lesssim 10^{11} M_\odot$ typically reionized as late or later than globally-averaged $50\%$ reionization at $\langle z_R\rangle =7.8$, in neighborhoods where reionization was completed by external radiation. The spread of reionization times within galaxies was sometimes as large as the galaxy-to-galaxy scatter. The Milky Way and M31 reionized earlier than global reionization but later than typical for their mass, neither dominated by external radiation. Their most massive progenitors at $z>6$ had $z_R$ = 9.8 (MW) and 11 (M31), while their total masses had $z_R$ = 8.2 (both).


title: ComPAS: Community Preserving Sampling for Streaming Graphs

abstract: In the era of big data, graph sampling is indispensable in many settings. Existing sampling methods are mostly designed for static graphs, and aim to preserve basic structural properties of the original graph (such as degree distribution, clustering coefficient etc.) in the sample. We argue that for any sampling method it is impossible to produce an universal representative sample which can preserve all the properties of the original graph; rather sampling should be application specific (such as preserving hubs - needed for information diffusion). Here we consider community detection as an application scenario. We propose ComPAS, a novel sampling strategy that unlike previous methods, is not only designed for streaming graphs (which is a more realistic representation of a real-world scenario) but also preserves the community structure of the original graph in the sample. Empirical results on both synthetic and different real-world graphs show that ComPAS is the best to preserve the underlying community structure with average performance reaching 73.2% of the most informed algorithm for static graphs.


title: Association schemes perspective of microbubble cluster in ultrasonic fields

abstract: Dynamics of a cluster of chaotic oscillators on a network are studied using coupled maps. By introducing the association schemes, we obtain coupling strength in the adjacency matrices form, which satisfies Markov matrices property. We remark that in general, the stability region of the cluster of oscillators at the synchronization state is characterized by Lyapunov exponent which can be defined based on the N-coupled map. As a detailed physical example, dynamics of microbubble cluster in an ultrasonic field are studied using coupled maps. Microbubble cluster dynamics have an indicative highly active nonlinear phenomenon, were not easy to be explained. In this paper, a cluster of microbubbles with a thin elastic shell based on the modified Keller-Herring equation in an ultrasonic field is demonstrated in the framework of the globally coupled map. On the other hand, a relation between the microbubble elements is replaced by a relation between the vertices. Based on this method, the stability region of microbubbles pulsations at complete synchronization state has been obtained analytically. In this way, distances between microbubbles as coupling strength play the crucial role. In the stability region, we thus observe that the problem of the study of dynamics of N-microbubble oscillators reduce to that of a single microbubble. Therefore, the important parameters of the isolated microbubble such as applied pressure, driving frequency and the initial radius have effective behavior on the synchronization state.


title: Re-Weighted Learning for Sparsifying Deep Neural Networks

abstract: This paper addresses the topic of sparsifying deep neural networks (DNN's). While DNN's are powerful models that achieve state-of-the-art performance on a large number of tasks, the large number of model parameters poses serious storage and computational challenges. To combat these difficulties, a growing line of work focuses on pruning network weights without sacrificing performance. We propose a general affine scaling transformation (AST) algorithm to sparsify DNN's. Our approach follows in the footsteps of popular sparse recovery techniques, which have yet to be explored in the context of DNN's. We describe a principled framework for transforming densely connected DNN's into sparsely connected ones without sacrificing network performance. Unlike existing methods, our approach is able to learn sparse connections at each layer simultaneously, and achieves comparable pruning results on the architecture tested.


title: Predictive Second Order Sliding Control of Constrained Linear Systems with Application to Automotive Control Systems

abstract: This paper presents a new predictive second order sliding controller (PSSC) formulation for setpoint tracking of constrained linear systems. The PSSC scheme is developed by combining the concepts of model predictive control (MPC) and second order discrete sliding mode control. In order to guarantee the feasibility of the PSSC during setpoint changes, a virtual reference variable is added to the PSSC cost function to calculate the closest admissible set point. The states of the system are then driven asymptotically to this admissible setpoint by the control action of the PSSC. The performance of the proposed PSSC is evaluated for an advanced automotive engine case study, where a high fidelity physics-based model of a reactivity controlled compression ignition (RCCI) engine is utilized to serve as the virtual test-bed for the simulations. Considering the hard physical constraints on the RCCI engine states and control inputs, simultaneous tracking of engine load and optimal combustion phasing is a challenging objective to achieve. The simulation results of testing the proposed PSSC on the high fidelity RCCI model show that the developed predictive controller is able to track desired engine load and combustion phasing setpoints, with minimum steady state error, and no overshoot. Moreover, the simulation results confirm the robust tracking performance of the PSSC during transient operations, in the presence of engine cyclic variability.


title: Unified Analysis and Optimization of D2D Communications in Cellular Networks Over Fading Channels

abstract: This paper develops an innovative approach to the modeling and analysis of downlink cellular networks with device-to-device (D$2$D) transmissions. The analytical embodiment of the signal-to-noise and-interference ratio (SINR) analysis in general fading channels is unified due to the H-transform theory, a taxonomy never considered before in stochastic geometry-based cellular network modeling and analysis. The proposed framework has the potential, due to versatility of the Fox's H functions, of significantly simplifying the cumbersome analysis procedure and representation of D$2$D and cellular coverage, while subsuming those previously derived for all the known simple and composite fading models. By harnessing its tractability, the developed statistical machinery is employed to launch an investigation into the optimal design of coexisting D$2$D and cellular communications. We propose novel coverage-aware power control combined with opportunistic access control to maximize the area spectral efficiency (ASE) of D$2$D communications. Simulation results substantiate performance gains achieved by the proposed optimization framework in terms of cellular communication coverage probability, average D$2$D transmit power, and the ASE of D$2$D communications under different fading models and link- and network-level dynamics.


title: Convergence of Maximum Bisection Ratio of Sparse Random Graphs

abstract: We consider sequences of large sparse random graphs whose degree distribution approaches a limit with finite mean. This model includes both the random regular graphs and the Erd\"os-Renyi graphs of constant average degree. We prove that the maximum bisection ratio of such a graph sequence converges almost surely to a deterministic limit. We extend this result to so-called 2-spin spin glasses in the paramagnetic to ferromagnetic regime. Our work generalizes the graph interpolation method to some non-additive graph parameters.


title: KLT-type relations for QCD and bicolor amplitudes from color-factor symmetry

abstract: Color-factor symmetry is used to derive a KLT-type relation for tree-level QCD amplitudes containing gluons and an arbitrary number of massive or massless quark-antiquark pairs, generalizing the expression for Yang-Mills amplitudes originally postulated by Bern, De Freitas, and Wong. An explicit expression is given for all amplitudes with two or fewer quark-antiquark pairs in terms of the (modified) momentum kernel. We also introduce the bicolor scalar theory, the "zeroth copy" of QCD, containing massless biadjoint scalars and massive bifundamental scalars, generalizing the biadjoint scalar theory of Cachazo, He, and Yuan. We derive KLT-type relations for tree-level amplitudes of biadjoint and bicolor theories using the color-factor symmetry possessed by these theories.


title: Un-unzippable Convex Caps

abstract: An unzipping of a polyhedron P is a cut-path through its vertices that unfolds P to a non-overlapping shape in the plane. It is an open problem to decide if every convex P has an unzipping. Here we show that there are nearly flat convex caps that have no unzipping. A convex cap is a "top" portion of a convex polyhedron; it has a boundary, i.e., it is not closed by a base.


title: A Matrix Baker-Akhiezer Function Associated with the Maxwell-Bloch Equations and their Finite-Gap Solutions

abstract: The Baker-Akhiezer (BA) function theory was successfully developed in the mid 1970s. This theory brought very interesting and important results in the spectral theory of almost periodic operators and theory of completely integrable nonlinear equations such as Korteweg-de Vries equation, nonlinear Schr\"odinger equation, sine-Gordon equation, Kadomtsev-Petviashvili equation. Subsequently the theory was reproduced for the Ablowitz-Kaup-Newell-Segur (AKNS) hierarchies. However, extensions of the Baker-Akhiezer function for the Maxwell-Bloch (MB) system or for the Karpman-Kaup equations, which contain prescribed weight functions characterizing inhomogeneous broadening of the main frequency, are unknown. The main goal of the paper is to give a such of extension associated with the Maxwell-Bloch equations. Using different Riemann-Hilbert problems posed on the complex plane with a finite number of cuts we propose such a matrix function that has unit determinant and takes an explicit form through Cauchy integrals, hyperelliptic integrals and theta functions. The matrix BA function solves the AKNS equations (the Lax pair for MB system) and generates a quasi-periodic finite-gap solution to the Maxwell-Bloch equations. The suggested function will be useful in the study of the long time asymptotic behavior of solutions of different initial-boundary value problems for the MB equations using the Deift-Zhou method of steepest descent and for an investigation of rogue waves of the Maxwell-Bloch equations.


title: Cost-optimal design of a simplified highly renewable Chinese electricity network

abstract: Rapid economic growth in China has lead to an increasing energy demand in the country. In combination with China's emission control and clean air initiatives, it has resulted in large-scale expansion of the leading renewable energy technologies, wind and solar power. Their intermittent nature and uneven geographic distribution, however, raises the question of how to best exploit them in a future sustainable electricity system, where their combined production may very well exceed that of all other technologies. It is well known that interconnecting distant regions provides more favorable production patterns from wind and solar. On the other hand, long-distance connections challenge traditional local energy autonomy. In this paper, the advantage of interconnecting the contiguous provinces of China is quantified. To this end, two different methodologies are introduced. The first aims at gradually increasing heterogeneity, that is non-local wind and solar power production, to minimize production costs without regard to the match between production and demand. The second method optimizes the trade-off between low cost production and high utility value of the energy. In both cases, the study of a 100% renewable Chinese electricity network is based on 8 years of high-resolution hourly time series of wind and solar power generation and electricity demand for each of the provinces. From the study we conclude that compared to a baseline design of homogeneously distributed renewable capacities, a heterogeneous network not only lowers capital investments but also reduces backup dispatches from thermal units. Installing more capacity in provinces like Inner Mongolia, Jiangsu, Hainan and north-western regions, heterogeneous layouts may lower the levelized cost of electricity (LCOE) by up to 27%, and reduce backup needs by up to 64%.


title: Irredundant Triangular Decomposition

abstract: Triangular decomposition is a classic, widely used and well-developed way to represent algebraic varieties with many applications. In particular, there exist sharp degree bounds for a single triangular set in terms of intrinsic data of the variety it represents, and powerful randomized algorithms for computing triangular decompositions using Hensel lifting in the zero-dimensional case and for irreducible varieties. However, in the general case, most of the algorithms computing triangular decompositions produce embedded components, which makes it impossible to directly apply the intrinsic degree bounds. This, in turn, is an obstacle for efficiently applying Hensel lifting due to the higher degrees of the output polynomials and the lower probability of success. In this paper, we give an algorithm to compute an irredundant triangular decomposition of an arbitrary algebraic set $W$ defined by a set of polynomials in C[x_1, x_2, ..., x_n]. Using this irredundant triangular decomposition, we were able to give intrinsic degree bounds for the polynomials appearing in the triangular sets and apply Hensel lifting techniques. Our decomposition algorithm is randomized, and we analyze the probability of success.


title: A posteriori error estimates for the Laplace-Beltrami operator on parametric $C^2$ surfaces

abstract: We prove new a posteriori error estimates for surface finite element methods (SFEM). Surface FEM approximate solutions to PDE posed on surfaces. Prototypical examples are elliptic PDE involving the Laplace-Beltrami operator. Typically the surface is approximated by a polyhedral or higher-order polynomial approximation. The resulting FEM exhibits both a geometric consistency error due to the surface approximation and a standard Galerkin error. A posteriori estimates for SFEM require practical access to geometric information about the surface in order to computably bound the geometric error. It is thus advantageous to allow for maximum flexibility in representing surfaces in practical codes when proving a posteriori error estimates for SFEM. However, previous a posteriori estimates using general parametric surface representations are suboptimal by one order on $C^2$ surfaces. Proofs of error estimates optimally reflecting the geometric error instead employ the closest point projection, which is defined using the signed distance function. Because the closest point projection is often unavailable or inconvenient to use computationally, a posteriori estimates using the signed distance function have notable practical limitations. We merge these two perspectives by assuming {\it practical} access only to a general parametric representation of the surface, but using the distance function as a {\it theoretical} tool. This allows us to derive sharper geometric estimators which exhibit improved experimentally observed decay rates when implemented in adaptive surface finite element algorithms.


title: Frobenius Heisenberg categorification

abstract: We associate a graded monoidal supercategory $\mathcal{H}\mathit{eis}_{F,k}$ to every graded Frobenius superalgebra $F$ and integer $k$. These categories, which categorify a broad range of lattice Heisenberg algebras, recover many previously defined Heisenberg categories as special cases. In this way, the categories $\mathcal{H}\mathit{eis}_{F,k}$ serve as a unifying and generalizing framework for Heisenberg categorification. Even in the case of previously defined Heisenberg categories, we obtain new, more efficient, presentations of these categories, based on an approach of Brundan. When $k=0$, our construction yields new versions of the affine oriented Brauer category depending on a graded Frobenius superalgebra.


title: Sensitivity of asymmetric rate-dependent critical systems to initial conditions: insights into cellular decision making

abstract: The work reported here aims to address the effects of time-dependent parameters and stochasticity on decision-making in biological systems. We achieve this by extending previous studies that resorted to simple normal forms. Yet, we focus primarily on the issue of the system's sensitivity to initial conditions in the presence of different noise distributions. In addition, we assess the impact of two-way sweeping through the critical region of a canonical Pitchfork bifurcation with a constant external asymmetry. The parallel with decision-making in bio-circuits is performed on this simple system since it is equivalent in its available states and dynamics to more complex genetic circuits. Overall, we verify that rate-dependent effects are specific to particular initial conditions. Information processing for each starting state is affected by the balance between sweeping speed through critical regions, and the type of fluctuations added. For a heavy-tail noise, forward-reverse dynamic bifurcations are more efficient in processing the information contained in external signals, when compared to the system relying on escape dynamics, if it starts at an attractor not favoured by the asymmetry and, in conjunction, if the sweeping amplitude is large.


title: Proposed Spreadsheet Transparency Definition and Measures

abstract: Auditors demand financial models be transparent yet no consensus exists on what that means precisely. Without a clear modeling transparency definition we cannot know when our models are "transparent". The financial modeling community debates which methods are more or less transparent as though transparency is a quantifiable entity yet no measures exist. Without a transparency measure modelers cannot objectively evaluate methods and know which improves model transparency. This paper proposes a definition for spreadsheet modeling transparency that is specific enough to create measures and automation tools for auditors to determine if a model meets transparency requirements. The definition also provides modelers the ability to objectively compare spreadsheet modeling methods to select which best meets their goals.


title: The Cohomology of Unramified Rapoport-Zink Spaces of EL-type and Harris's Conjecture

abstract: We study the $l$-adic cohomology of unramified Rapoport-Zink spaces of EL-type. These spaces were used in Harris and Taylor's proof of the local Langlands correspondence for $\mathrm{GL_n}$ and to show local-global compatibilities of the Langlands correspondence. In this paper we consider certain morphisms, $\mathrm{Mant}_{b, \mu}$, of Grothendieck groups of representations constructed from the cohomology of the above spaces, as studied by Harris and Taylor, Mantovan, Fargues, Shin, and others. Due to earlier work of Fargues and Shin we have a description of $\mathrm{Mant}_{b, \mu}(\rho)$ for $\rho$ a supercuspidal representation. In this paper, we give a conjectural formula for $\mathrm{Mant}_{b, \mu}(\rho)$ for all $\rho$ and prove it when $\rho$ is essentially square integrable. Our proof works for general $\rho$ conditionally on a conjecture appearing in Shin's work. We show that our description agrees with a conjecture of Harris in the case of parabolic inductions of supercuspidal representations of a Levi subgroup.


title: Estimation of Viterbi path in Bayesian hidden Markov models

abstract: The article studies different methods for estimating the Viterbi path in the Bayesian framework. The Viterbi path is an estimate of the underlying state path in hidden Markov models (HMMs), which has a maximum posterior probability (MAP). For an HMM with given parameters, the Viterbi path can be easily found with the Viterbi algorithm. In the Bayesian framework the Viterbi algorithm is not applicable and several iterative methods can be used instead. We introduce a new EM-type algorithm for finding the MAP path and compare it with various other methods for finding the MAP path, including the variational Bayes approach and MCMC methods. Examples with simulated data are used to compare the performance of the methods. The main focus is on non-stochastic iterative methods and our results show that the best of those methods work as well or better than the best MCMC methods. Our results demonstrate that when the primary goal is segmentation, then it is more reasonable to perform segmentation directly by considering the transition and emission parameters as nuisance parameters.


title: A Discrete Spherical Harmonics Method for Radiative Transfer Analysis in Inhomogeneous Polarized Planar Atmosphere

abstract: A discrete spherical harmonics method is developed for the radiative transfer problem in inhomogeneous polarized planar atmosphere illuminated at the top by a collimated sunlight while the bottom reflects the radiation. The method expands both the Stokes vector and the phase matrix in a finite series of generalized spherical functions and the resulting vector radiative transfer equation is expressed in a set of polar directions. Hence, the polarized characteristics of the radiance within the atmosphere at any polar direction and azimuthal angle can be determined without linearization and/or interpolations. The spatial dependent of the problem is solved using the spectral Chebyshev method. The emergent and transmitted radiative intensity and the degree of polarization are predicted for both Rayleigh and Mie scattering. The discrete spherical harmonics method predictions using 36 streams are found in good agreement with benchmark literature results. The maximum deviation between the proposed method and literature results and for absolue value of the polar directions greater than 0.1 is less than 0.5% and 0.9% for the Rayleigh and Mie scattering, respectively. These deviations for directions close to zero are about 3% and 10 % for Rayleigh and Mie scattering, respectively.


title: Electrical-current-induced magnetic hysteresis in self-assembled vertically aligned La_{2/3}Sr_{1/3}MnO_3:ZnO-nanopillar composites

abstract: Magnetoresistive random-access memory (MRAM) is poised to become a next-generation information storage device. Yet, many materials challenges remain unsolved before it can become a widely used memory storage solution. Among them, an urgent need is to identify a material system that is suitable for downscaling and is compatible with low-power logic applications. Self-assembled, vertically-aligned La_{2/3}Sr_{1/3}MnO_3:ZnO nanocomposites, in which La_{2/3}Sr_{1/3}MnO_3 (LSMO) matrix and ZnO nanopillars form an intertwined structure with coincident-site-matched growth occurring between the LSMO and ZnO vertical interfaces, may offer new MRAM applications by combining their superior electric, magnetic (B), and optical properties. In this paper, we show the results of electrical current induced magnetic hysteresis in magneto-resistance measurements in these nano-pillar composites. We observe that when the current level is low, for example, 1 uA, the magneto-resistance displays a linear, negative, non-hysteretic B field dependence. Surprisingly, when a large current is used, I > 10 uA, a hysteretic behavior is observed when the B field is swept in the up and down directions. This hysteresis weakens as the sample temperature is increased. A possible spin-valve mechanism related to this electrical current induced magnetic hysteresis is proposed and discussed.


title: Topological phases, edge modes, and the Hofstadter butterfly in coupled Su-Schrieffer-Heeger systems

abstract: Motivated by recent experimental realizations of topological edge states in Su-Schrieffer-Heeger (SSH) chains, we theoretically study a ladder system whose legs are comprised of two such chains. We show that the ladder hosts a rich phase diagram and related edge mode structure dictated by choice of inter-chain and intra-chain couplings. Namely, we exhibit three distinct physical regimes: a topological hosting localized zero energy edge modes, a topologically trivial phase having no edge mode structure, and a regime reminiscent of a weak topological insulator having unprotected edge modes resembling a "twin-SSH" construction. In the topological phase, the SSH ladder system acts as an analog of the Kitaev chain, which is known to support localized Majorana fermion end modes, with the difference that bound states of the SSH ladder having the same spatial wavefunction profiles correspond to Dirac fermion modes. Further, inhomogeneity in the couplings can have a drastic effect on the topological phase diagram of the ladder system. In particular for quasiperiodic variations of the inter-chain coupling, the phase diagram reproduces Hofstadter's butterfly pattern. We thus identify the SSH ladder system as a potential candidate for experimental observation of such fractal structure.


title: Quantum anomaly and thermodynamics of one-dimensional fermions with three-body interactions

abstract: We show that a system of three species of one-dimensional fermions, with an attractive three-body contact interaction, features a scale anomaly directly related to the anomaly of two-dimensional fermions with two-body forces. We show, furthermore, that those two cases (and their multi species generalizations) are the only non-relativistic systems with contact interactions that display a scale anomaly. While the two-dimensional case is well-known and has been under study both experimentally and theoretically for years, the one-dimensional case presented here has remained unexplored. For the latter, we calculate the impact of the anomaly on the equation of state, which appears through the generalization of Tan's contact for three-body forces, and determine the pressure at finite temperature. In addition, we show that the third-order virial coefficient is proportional to the second-order coefficient of the two-dimensional two-body case.

categories: cond-mat.quant-gas hep-th nucl-th quant-ph


title: Imaging a non-singular rotating black hole at the center of the Galaxy

abstract: We show that the rotating generalization of Hayward's non-singular black hole previously studied in the literature is geodesically incomplete, and that its straightforward extension leads to a singular spacetime. We present another extension, which is devoid of any curvature singularity. The obtained metric depends on three parameters and, depending on their values, yields an event horizon or not. These two regimes, named respectively regular rotating Hayward black hole and naked rotating wormhole, are studied both numerically and analytically. In preparation for the upcoming results of the Event Horizon Telescope, the images of an accretion torus around Sgr A*, the supermassive object at the center of the Galaxy, are computed. These images contain, even in the absence of a horizon, a central faint region which bears a resemblance to the shadow of Kerr black holes and emphasizes the difficulty of claiming the existence of an event horizon from the analysis of strong-field images. The frequencies of the co- and contra-rotating orbits at the innermost stable circular orbit (ISCO) in this geometry are also computed, in the hope that quasi-periodic oscillations may permit to compare this model with Kerr's black hole on observational grounds.


title: Do You Want Your Autonomous Car To Drive Like You?

abstract: With progress in enabling autonomous cars to drive safely on the road, it is time to start asking how they should be driving. A common answer is that they should be adopting their users' driving style. This makes the assumption that users want their autonomous cars to drive like they drive - aggressive drivers want aggressive cars, defensive drivers want defensive cars. In this paper, we put that assumption to the test. We find that users tend to prefer a significantly more defensive driving style than their own. Interestingly, they prefer the style they think is their own, even though their actual driving style tends to be more aggressive. We also find that preferences do depend on the specific driving scenario, opening the door for new ways of learning driving style preference.


title: On residue maps for affine curves

abstract: We establish several compatibility results between residue maps in \'etale and Galois cohomology that arise naturally in the analysis of smooth affine algebraic curves having good reduction over discretely valued fields. These results are needed, and in fact have already been used, for the study of finiteness properties of the unramified cohomology of function fields of affine curves over number fields.


title: Quantum many-body dynamics of the Einstein-de Haas effect

abstract: In 1915, Einstein and de Haas and Barnett demonstrated that changing the magnetization of a magnetic material results in mechanical rotation, and vice versa. At the microscopic level, this effect governs the transfer between electron spin and orbital angular momentum, and lattice degrees of freedom, understanding which is key for molecular magnets, nano-magneto-mechanics, spintronics, and ultrafast magnetism. Until now, the timescales of electron-to-lattice angular momentum transfer remain unclear, since modeling this process on a microscopic level requires addition of an infinite amount of quantum angular momenta. We show that this problem can be solved by reformulating it in terms of the recently discovered angulon quasiparticles, which results in a rotationally invariant quantum many-body theory. In particular, we demonstrate that non-perturbative effects take place even if the electron--phonon coupling is weak and give rise to angular momentum transfer on femtosecond timescales.


title: A Data as a Service (DaaS) Model for GPU-based Data Analytics

abstract: Cloud-based services with resources to be provisioned for consumers are increasingly the norm, especially with respect to Big data, spatiotemporal data mining and application services that impose a user's agreed Quality of Service (QoS) rules or Service Level Agreement (SLA). Considering the pervasive nature of data centers and cloud system, there is a need for a real-time analytics of the systems considering cost, utility and energy. This work presents an overlay model of GPU system for Data As A Service (DaaS) to give a real-time data analysis of network data, customers, investors and users' data from the datacenters or cloud system. Using a modeled layer to define a learning protocol and system, we give a custom, profitable system for DaaS on GPU. The GPU-enabled pre-processing and initial operations of the clustering model analysis is promising as shown in the results. We examine the model on real-world data sets to model a big data set or spatiotemporal data mining services. We also produce results of our model with clustering, neural networks' Self-organizing feature maps (SOFM or SOM) to produce a distribution of the clustering for DaaS model. The experimental results thus far show a promising model that could enhance SLA and or QoS based DaaS.


title: Mitigating Spreadsheet Risk in Complex Multi-Dimensional Models in Excel

abstract: Microsoft Excel is the most ubiquitous analytical tool ever built. Companies around the world leverage it for its power, flexibility and ease of use. However, spreadsheets are manually intensive and prone to error, making it difficult for companies to control spreadsheet risk. The following solution is designed to mitigate spreadsheet risk for a set of problems commonly addressed in a spreadsheet defined as "complex multi-dimensional models". "Complex" referring to certain types of applications that require functionality such as sophisticated algorithms, challenging hierarchies and database write-back (i.e. planning, forecasting, etc.) and "multi-dimensional" referring to providing capabilities such as reporting, data input forms and ad hoc analysis on the different attributes associated with the resulting model. The solution is defined as a "PivotModel" because it works similarly to a PivotTable but is designed to leverage the robust capabilities of the Microsoft Excel platform.


title: Volatility options in rough volatility models

abstract: We discuss the pricing and hedging of volatility options in some rough volatility models. First, we develop efficient Monte Carlo methods and asymptotic approximations for computing option prices and hedge ratios in models where log-volatility follows a Gaussian Volterra process. While providing a good fit for European options, these models are unable to reproduce the VIX option smile observed in the market, and are thus not suitable for VIX products. To accommodate these, we introduce the class of modulated Volterra processes, and show that they successfully capture the VIX smile.


title: Modal Noise Mitigation through Fiber Agitation for Fiber-fed Radial Velocity Spectrographs

abstract: Optical fiber modal noise is a limiting factor for high precision spectroscopy signal-to-noise in the near-infrared and visible. Unabated, especially when using highly coherent light sources for wavelength calibration, modal noise can induce radial velocity (RV) errors that hinder the discovery of low-mass (and potentially Earth-like) planets. Previous research in this field has found sufficient modal noise mitigation through the use of an integrating sphere, but this requires extremely bright light sources, a luxury not necessarily afforded by the next generation of high-resolution optical spectrographs. Otherwise, mechanical agitation, which "mixes" the fiber's modal patterns and allows the noise to be averaged over minutes-long exposures, provides some noise reduction but the exact mechanism behind improvement in signal-to-noise and RV drift has not been fully explored or optimized by the community. Therefore, we have filled out the parameter space of modal noise agitation techniques in order to better understand agitation's contribution to mitigating modal noise and to discover a better method for agitating fibers. We find that modal noise is best suppressed by the quasi-chaotic motion of two high-amplitude agitators oscillating with varying phase for fibers with large core diameters and low azimuthal symmetry. This work has subsequently influenced the design of a fiber agitator, to be installed with the EXtreme PREcision Spectrograph, that we estimate will reduce modal-noise-induced RV error to less than 3.2 cm/s.


title: $C^{1,\alpha}$ regularity for fully nonlinear elliptic equations with superlinear growth in the gradient

abstract: We extend the Caffarelli-\'Swiech-Winter $C^{1,\alpha}$ regularity estimates to $L^p$-viscosity solutions of fully nonlinear uniformly elliptic equations in nondivergence form with superlinear growth in the gradient and unbounded coefficients. As an application, in addition to the usual $W^{2,p}$ results, we prove the existence of positive eigenvalues for proper operators with nonnegative unbounded weight, in particular for Pucci's operators with unbounded coefficients.


title: Nonlinear effects for three-terminal heat engine and refrigerator

abstract: The three-terminal heat device consisting of a cavity and coupled to a heat bath is established. By tuning the temperatures of the electrodes and the phonon bath, the device can function as a heat engine or a refrigerator. We study the characteristic performance in the linear and nonlinear regime for both setups. It is our focus here to analyze how the efficiency of the heat engine and coefficient of performance of the refrigerator are affected by the nonlinear transport. With such considerations, the maximum efficiency and power are then optimized for various energy levels, temperatures and other parameters.


title: Semi-Transparent Solar Cell enabled by Frequency Selective Light Trapping

abstract: We propose a frequency selective light trapping scheme that enables the creation of more visually-transparent and yet simultaneously more efficient semitransparent solar cells. A nanoparticle scattering layer and photonic stack back reflector create a selective trapping effect by total internal reflection within a medium, increasing absorption of IR light. We propose a strong frequency selective scattering layer using spherical TiO2 nanoparticles with radius of 255 nm and area density of 1.1% in a medium with index of refraction of 1.5. Using detailed numerical simulations for this configuration, we find that it is possible to create a semitransparent silicon solar cell that has a Shockley Queisser efficiency of 12.0%\pm0.4% with a visible transparency of 60.2%\pm1.3%, 13.3%\pm1.3 more visibly-transparent than a bare silicon cell at the same efficiency.


title: Anagram-Free Chromatic Number is not Pathwidth-Bounded

abstract: The anagram-free chromatic number is a new graph parameter introduced independently Kam\v{c}ev, {\L}uczak, and Sudakov (2017) and Wilson and Wood (2017). In this note, we show that there are planar graphs of pathwidth 3 with arbitrarily large anagram-free chromatic number. More specifically, we describe $2n$-vertex planar graphs of pathwidth 3 with anagram-free chromatic number $\Omega(\log n)$. We also describe $kn$ vertex graphs with pathwidth $2k-1$ having anagram-free chromatic number in $\Omega(k\log n)$.


title: Strain-mediated spin-orbit torque switching for magnetic memory

abstract: Spin-orbit torque (SOT) represents an energy efficient method to control magnetization in magnetic memory devices. However, deterministically switching perpendicular memory bits usually requires the application of an additional bias field for breaking lateral symmetry. Here we present a new approach of field-free deterministic perpendicular switching using a strain-mediated SOT switching method. The strain-induced magnetoelastic anisotropy breaks the lateral symmetry, and the resulting symmetry-breaking is controllable. A finite element model and a macrospin model are used to numerically simulate the strain-mediated SOT switching mechanism. The results show that a relatively small voltage (${\pm}0.5$ V) along with a modest current ($3.5 \times 10^{7} A/cm^{2}$) can produce a 180{\deg} perpendicular magnetization reversal. The switching direction (up or down) is dictated by the voltage polarity (positive or negative) applied to the piezoelectric layer in the magnetoelastic/heavy metal/piezoelectric heterostructure. The switching speed can be as fast as 10 GHz. More importantly, this control mechanism can be potentially implemented in a magnetic random-access memory system with small footprint, high endurance and high tunnel magnetoresistance (TMR) readout ratio.


title: Surface Tension Prediction for Liquid Mixtures

abstract: Recently Escobedo and Mansoori (AIChE J. 42(5): 1425, 1996) proposed a new expression, originally derived from statistical mechanics, for the surface tension prediction of pure fluids. In this report, the same theory is extended to the case of mixtures of organic liquids. Using the proposed equation surface tensions of 55 binary mixtures are predicted within an overall 0.50 AAD% which is better than the other available prediction methods. Using a corresponding states expression the surface tension of all the 55 binary mixtures are predicted within an overall 2.10 AAD%. The proposed model is also applicable to multicomponent liquid mixtures.


title: Real-Time Rejection and Mitigation of Time Synchronization Attacks on the Global Positioning System

abstract: This paper introduces the Time Synchronization Attack Rejection and Mitigation (TSARM) technique for Time Synchronization Attacks (TSAs) over the Global Positioning System (GPS). The technique estimates the clock bias and drift of the GPS receiver along with the possible attack contrary to previous approaches. Having estimated the time instants of the attack, the clock bias and drift of the receiver are corrected. The proposed technique is computationally efficient and can be easily implemented in real time, in a fashion complementary to standard algorithms for position, velocity, and time estimation in off-the-shelf receivers. The performance of this technique is evaluated on a set of collected data from a real GPS receiver. Our method renders excellent time recovery consistent with the application requirements. The numerical results demonstrate that the TSARM technique outperforms competing approaches in the literature.


title: Centimeter-scale superfine three-dimensional printing with femtosecond laser two-photon polymerization

abstract: Nowadays three-dimensional (3D) printing has been widely used for producing geometrically complex 3D structures from a broad range of materials such as ceramics, metals, polymers, semiconductors, etc. Although it has been demonstrated that a fabrication resolution as high as ~100 nm can be achieved in 3D printing based on two photon polymerization (TPP), the end product size of TPP is typically on millimeter scale limited by the short working distance of high-numerical-aperture focal lens. Here we present a method based on simultaneous spatiotemporal focusing (SSTF) of the femtosecond laser pulses that enables to fabricate centimeter-scale 3D structures of fine features with TPP. We also demonstrate an isotropic spatial resolution which can be continuously tuned in the range of ~10 um and ~40 um by only varying the power of femtosecond laser, making this technique extremely flexible and easy to implement. We fabricate several Chinese guardian lions of a maximum height of 0.6 cm and a Terra Cotta Warrior of a height of 1.3 cm using this method.


title: High sensitivity refractive index sensor based on the semicircular bent fiber

abstract: A refractive index sensor based on a semicircular bent fiber is presented. The interference occurs between the cladding mode excited in the bending region and the core mode. Both the theoretical and experimental results show that the resonant dip wavelength decreases linearly with the increase of the refractive index of the surrounding environment. A high sensitivity of 1031 nm per refractive index unit is obtained over the refractive index range of 1.3324 to 1.3435 by using a bent fiber with a bending radius of 500 {\mu}m.


title: Revisiting the charging-capacitor problem

abstract: The charging capacitor is used as a standard paradigm for illustrating the concept of the Maxwell "displacement current". A certain aspect of the problem, however, is often overlooked. It concerns the conditions for satisfaction of the Faraday-Henry law both in the interior and in the exterior of the capacitor. In this article the situation is analyzed and a recursive process is described for obtaining (at least approximate) solutions of Maxwell's equations inside and outside the capacitor.


title: Levy walk with multiple internal states

abstract: Levy walk is a fundamental model with applications ranging from quantum physics to paths of animal foraging. Taking animal foraging as an example, a natural idea that comes to one's mind is to introduce the multiple internal states for dealing with the dependence of the PDF of waiting time on the energy of the animal and richness of the food at a particular location, etc; the framework can also be used to model the moving trajectories of smart animals without returning to the directions or locations which they come from immediately. After building the Levy walk model with multiple internal states and deriving the governing equation of the distribution of the positions of the particles, some applications are discussed with specific transition matrices. The type of diffusion for non-immediately-repeating L\'{e}vy walk is uncovered, and the distribution and average of first passage time are numerically simulated.


title: Quantum mechanical analysis of nonlinear optical response of interacting graphene nanoflakes

abstract: We propose a distant-neighbor quantum-mechanical (DNQM) approach to study the linear and nonlinear optical properties of graphene nanoflakes (GNFs). In contrast to the widely used tight-binding description of the electronic states that considers only the nearest-neighbor coupling between the atoms, our approach is more accurate and general, as it captures the electron-core interactions between all atoms in the structure. Therefore, as we demonstrate, the DNQM approach enables the investigation of the optical coupling between two closely separated, but chemically unbound GNFs. We also find that the optical response of GNFs depends crucially on their shape, size, and symmetry properties. Specifically, increasing the size of nanoflakes is found to shift their accommodated quantum plasmon oscillations to lower frequency. Importantly, we show that by embedding a cavity into GNFs, one can change their symmetry properties, tune their optical properties, or enable otherwise forbidden second-harmonic generation processes.


title: 3D Printing of Flexible Room-Temperature Liquid Metal Battery

abstract: The entirely soft and transformable room-temperature liquid metal battery based on 3D hybrid printing was proposed and experimentally demonstrated. Liquid metal gallium and conductive gel were employed as the negative and positive electrode material, respectively. The discharge process of batteries with different structure, single-cell and multi-cell, were clarified. The results indicated that the obtained battery possessed rather stable and persisted discharge property, but the multi-cell battery presented stair-step decrease overtime. SEM and EDS images were obtained to evaluate the morphological change of the electrode during the discharge process, and the liquid metal was found to be easily exposed and active according to the solubility of the alkaline electrolyte to oxidized gallium. The successful 3D hybrid printing allowed different polymers, including PLA, ABS and TPU, to be precisely patterned in any required battery body structures, with liquid electrodes in appropriate specified positions. The flexibility and practicability of this new generation battery was illustrated through working devices. The present soft battery offers important power supply for future wearable and epidermal electronics.


title: Innovative thermo-solar air heater

abstract: In the present work we elaborate the innovative design of the solar air heater and justify it by a Computational Fluid Dynamics (CFD) simulation, implementing and experimentally testing a sample. We propose to use this device for maintenance of constant ambient conditions for thermal comfort and low energy consumption for indoor environments, inside greenhouses, passive houses, and to protect buildings against temperature fluctuations. We tested the functionality of our sample of the solar air heater for 50 weeks and obtained an agreement between the results of the numerical simulation, implemented using OpenFOAM (an open source numerical CFD software) and the experimental results.


title: Hybrid Nanophotonics

abstract: Advances in the field of plasmonics, that is, nanophotonics based on optical properties of metal nanostructures, paved the way for the development of ultrasensitive biological sensors and other devices whose operating principles are based on localization of an electromagnetic field at the nanometer scale. However, high dissipative losses of metal nanostructures limit their performance in many modern areas, including metasurfaces, metamaterials, and optical interconnections, which required the development of new devices that combine them with high refractive index dielectric nanoparticles. Resulting metal-dielectric (hybrid) nanostructures demonstrated many superior properties from the point of view of practical application, including moderate dissipative losses, resonant optical magnetic response, strong nonlinear optical properties, which made the development in this field the vanguard of the modern light science. This review is devoted to the current state of theoretical and experimental studies of hybrid metal-dielectric nanoantennas and nanostructures based on them, capable of selective scattering light waves, amplifying and transmitting optical signals in the desired direction, controlling the propagation of such signals, and generating optical harmonics.


title: Almost Hyperbolic Groups with Almost Finitely Presented Subgroups

abstract: We construct new examples of CAT(0) groups containing non finitely presented subgroups that are of type $FP_2$, these CAT(0) groups do not contain copies of $\mathbb{Z}^3$. We also give a construction of groups which are of type $F_n$ but not $F_{n+1}$ with no free abelian subgroups of rank greater than $\lceil\frac{n}{3}\rceil$.


title: Transmutation of Minor Actinides and Power Flattening in PWR MOX Fuel

abstract: In order to transmute long-lived MAs and flatten power distribution in a PWR MOX fuel assembly, the authors have proposed to replace some high power fuel rods by MAs loaded fuel rods. The present work proves the high efficiency of long-lived MAs transmutation in a PWR MOX assembly with 92 fuel rods with 3% MAs loaded. The power and the burnup distribution have been flatted by using the MAs loading method proposed by authors. The 237Np loading method is expected for the transmutation by comparing with natural decay time to achieve the same reduction. The mixed MAs loading has better behaviors on power flattening and negative contribution of reactivity. In addition, the mixed MAs loading can largely reduce the quantity of 241Am, while the addition of other MAs has no influence on the transmutation efficiency of 237Np.


title: Significant efficiency enhancement in thin film solar cells using laser beam-induced graphene transparent conductive electrodes

abstract: Thin film solar cells have been attractive for decades in advanced green technology platforms due to its possibilities to be integrated with buildings and on-chip applications. However, the bottleneck issues involved to consider the current solar cells as a major electricity source includes the lower efficiencies and cost-effectiveness. We numerically demonstrate the concept of the absorption enhancement in thin-film amorphous silicon solar cells using the laser beam-induced graphene material based on the insensitive polarization space-filling fractal design as transparent conductive electrodes. With the optimization of parameters such as thickness, width, and period of fractals, an enhancement of photocurrent generation of solar cells by a factor of 24.5% is achieved compared to reference solar cell with a traditional ITO.


title: A priori bounds and multiplicity for fully nonlinear equations with quadratic growth in the gradient

abstract: We consider fully nonlinear uniformly elliptic equations with quadratic growth in the gradient, such as $$ -F(x,u,Du,D^2u) =\lambda c(x)u+\langle M(x)D u, D u \rangle +h(x) $$ in a bounded domain with a Dirichlet boundary condition, here $\lambda \in\mathbb{R}$, $c,\, h \in L^p(\Omega)$, $p>n\geq 1$, $c\gneqq 0$ and the matrix $M$ satisfies $0<\mu_1 I\leq M\leq \mu_2 I$. Recently this problem was studied in the "coercive" case $\lambda c\le0$, where uniqueness of solutions can be expected, and it was conjectured that the solution set is more complex for noncoercive equations. This conjecture was verified in 2015 by Arcoya, de Coster, Jeanjean and Tanaka for equations in divergence form, by exploiting the integral formulation of the problem. Here we show that similar phenomena occur for general, even fully nonlinear, equations in nondivergence form. We use different techniques based on the maximum principle. We develop a new method to obtain the crucial uniform a priori bounds, which permit to us to use degree theory. This method is based on basic regularity estimates such as half-Harnack inequalities, and on a V\'azquez type strong maximum principle for our kind of equations.


title: Characterizing Germanium Junction Transistors

abstract: Transistors have provided the basis of modern electronics. Being relatively intricate devices, and often exhibiting intense parameter variation, this type of electronic devices has motivated much research interest especially regarding their characterization and modeling. In this work, we apply a recently reported modeling methodology, based on the Early effect, for characterizing new old stock NPN and PNP small signal germanium junction transistors and comparing them to more modern silicon bipolar junction devices. The Early approach is special in the sense that its two parameters, namely the Early voltage $V_a$ and a proportionality parameter $s$, are fixed and independent of transistor operation. Remarkable results are obtained, including the fact that the four considered groups, namely NPN and PNP germanium and silicon devices, occupy mostly non-overlapping regions in the Early parameter space, with PNP devices presenting larger parameter variability in both cases. Surprisingly, the considered germanium devices exhibited smaller parameter variability than observed for the silicon counterparts. When mapped into the more traditional space defined by the current gain and output resistance parameters, the four transistor groups exhibited much larger overlaps and yielded a clustering structure much less organized than allowed by the Early mapping. This result suggest that the Early representation of transistors is more compatible and inherently related to the structure of amplifying devices such as those considered in this work. In addition, it was verified that the center of mass of each of the NPN-PNP pairs of germanium and silicon devices are crossed by respective $\beta$-isolines for gains of 130 and 250, respectively. Germanium devices were also characterized as having smaller output resistance and smaller magnitudes of Early voltage.


title: Super Scatterers Based on Artificial Localized Magnon Resonances

abstract: The interaction between electromagnetic waves and objects is strongly affected by the shape and material composition of the latter. Artificially created materials, formed by a subwavelength structuring of their unit cells, namely metamaterials, can exhibit peculiar responses to electromagnetic radiation and provide additional powerful degrees of freedom to the scatterer design. In particular, negative material susceptibilities give rise to strong resonant interactions with deeply subwavelength particles. While the negative electrical permittivity of natural noble metals manifests itself in localized plasmon resonant oscillations, negative magnetic permeability is virtually non-existent in nature. Here the concept of artificial magnon resonance in subwavelength objects with effective negative permeability, designed based on the metamaterial approach, is demonstrated. Strong localized oscillations of the magnetic fields within an array of split ring resonators, forming a sphere, hybridize in a collective mode of the structure. As a result, extremely high scattering cross section, exceeding that of a steel sphere with the same radius by four orders of magnitude, was demonstrated. Furthermore, the scattering cross section of subwavelength metamaterial-based sphere was shown to be comparable to the low frequency (MHz) radar signature of a big military aircraft. Super scatterers, based on tunable resonances within artificially created materials, can find use in a broad range of electromagnetic applications, including wireless communications, radars, RFID, internet of things hardware and many others.


title: Bistability of buoyancy-driven exchange flows in vertical tubes

abstract: Buoyancy-driven exchange flows are common to a variety of natural and engineering systems ranging from persistently active volcanoes to counterflows in oceanic straits. Experiments of exchange flows in closed vertical tubes have been used as surrogates to elucidate the basic features of such flows. The resulting data have historically been analyzed and interpreted through core-annular flow solutions, the most common flow configuration at finite viscosity contrasts. These models have been successful in fitting experimental data, but less effective at explaining the variability observed in natural systems. In this paper, we formulate a core-annular solution to the classical problem of buoyancy-driven exchange flows in vertical tubes. The model posits the existence of two mathematically valid solutions, i.e. thin- and thick-core solutions. The theoretical existence of two solutions, however, does not necessarily imply that the system is bistable in the sense that flow switching may occur. Using direct numerical simulations, we test the hypothesis that core-annular flow in vertical tubes is bistable, which implies that the realized flow field is not uniquely defined by the material parameters of the flow. Our numerical experiments, which fully predict experimental data without fitting parameters, demonstrate that buoyancy-driven exchange flows are indeed inherently bistable systems. This finding is consistent with previous experimental data, but in contrast to the underlying hypothesis of previous analytical models that the solution is unique and can be identified by maximizing the flux or extremizing the dissipation in the system. These results have important implications for data interpretation by analytical models, and may also have relevant ramifications for understanding volcanic degassing.


title: Particle Production in Ultra-Strong Coupling Waveguide QED

abstract: Understanding large-scale interacting quantum matter requires dealing with the huge number of quanta that are produced by scattering even a few particles against a complex quantum object. Prominent examples are found from high energy cosmic ray showers to the optical or electrical driving of degenerate Fermi gases. We tackle this challenge in the context of many-body quantum optics, as motivated by the recent developments of circuit quantum electrodynamics at ultrastrong coupling. The issue of particle production is addressed quantitatively with a simple yet powerful concept rooted in the quantum superposition principle. This key idea is illustrated by the study of multi-photon emission from a single two-level artificial atom coupled to a high impedance waveguide. We find surprisingly that the off-resonant inelastic emission lineshape is dominated by broadband particle production, due to the large phase space associated with contributions that do not conserve the number of excitations. Such frequency conversion processes produce striking signatures in time correlation measurements, which can be tested experimentally in quantum waveguides. These ideas open new directions for the simulation of a variety of physical systems, from polaron dynamics in solids to complex superconducting quantum architectures.