Biological evolution provides a creative fount of complex and subtle adaptations, often surprising the scientists who discover them. However, because evolution is an algorithmic process that transcends the substrate in which it occurs, evolution's creativity is not limited to nature. Indeed, many researchers in the field of digital evolution have observed their evolving algorithms and organisms subverting their intentions, exposing unrecognized bugs in their code, producing unexpected adaptations, or exhibiting outcomes uncannily convergent with ones in nature. Such stories routinely reveal creativity by evolution in these digital worlds, but they rarely fit into the standard scientific narrative. Instead they are often treated as mere obstacles to be overcome, rather than results that warrant study in their own right. The stories themselves are traded among researchers through oral tradition, but that mode of information transmission is inefficient and prone to error and outright loss. Moreover, the fact that these stories tend to be shared only among practitioners means that many natural scientists do not realize how interesting and lifelike digital organisms are and how natural their evolution can be. To our knowledge, no collection of such anecdotes has been published before. This paper is the crowd-sourced product of researchers in the fields of artificial life and evolutionary computation who have provided first-hand accounts of such cases. It thus serves as a written, fact-checked collection of scientifically important and even entertaining stories. In doing so we also present here substantial evidence that the existence and importance of evolutionary surprises extends beyond the natural world, and may indeed be a universal property of all complex evolving systems.

We perform a search for stellar streams around the Milky Way using the first three years of multi-band optical imaging data from the Dark Energy Survey (DES). We use DES data covering $\sim 5000$ sq. deg. to a depth of g > 23.5 with a relative photometric calibration uncertainty of < 1 \%. This data set yields unprecedented sensitivity to the stellar density field in the southern celestial hemisphere, enabling the detection of faint stellar streams to a heliocentric distance of $\sim 50$ kpc. We search for stellar streams using a matched-filter in color-magnitude space derived from a synthetic isochrone of an old, metal-poor stellar population. Our detection technique recovers four previously known thin stellar streams: Phoenix, ATLAS, Tucana III, and a possible extension of Molonglo. In addition, we report the discovery of eleven new stellar streams. In general, the new streams detected by DES are fainter, more distant, and lower surface brightness than streams detected by similar techniques in previous photometric surveys. As a by-product of our stellar stream search, we find evidence for extra-tidal stellar structure associated with four globular clusters: NGC 288, NGC 1261, NGC 1851, and NGC 1904. The ever-growing sample of stellar streams will provide insight into the formation of the Galactic stellar halo, the Milky Way gravitational potential, as well as the large- and small-scale distribution of dark matter around the Milky Way.

We propose and study a family of complex matrix models computing the protected two- and three-point correlation functions in $\mathcal{N}=4$ SYM. Our description allows us to directly relate the eigenvalue density of the matrix model for ``Huge" operators with $\Delta \sim N^2$ to the shape of droplets in the dual Lin-Lunin-Maldacena (LLM) geometry. We demonstrate how to determine the eigenvalue distribution for various choices of operators such as those of exponential, character, or coherent state type, which then allows us to efficiently compute one-point functions of light chiral primaries in generic LLM backgrounds. In particular, we successfully match the results for light probes with the supergravity calculations of Skenderis and Taylor. We provide a large $N$ formalism for one-point functions of ``Giant" probes, such as operators dual to giant graviton branes in LLM backgrounds, and explicitly apply it for particular backgrounds. We also explicitly compute the correlator of three huge half-BPS operators of exponential type and stacks of determinant operators by reducing them to the known matrix model problems such as the Potts or $O(n)$ model on random planar graphs. Finally, we point out a curious relation between the correlators of $\frac{1}{4}$-BPS and $\frac{1}{8}$-BPS coherent state operators and the Eguchi-Kawai reduction of the Principal Chiral Model in $2D$ and $3D$ correspondingly.

We consider the observation of stellar-mass black holes binaries with the Laser Interferometer Space Antenna (LISA). Preliminary results based on Fisher information matrix analyses have suggested that gravitational waves from those sources could be very sensitive to possible deviations from the theory of general relativity and from the strong equivalence principle during the low-frequency binary inspiral. We perform a full Markov Chain Monte Carlo Bayesian analysis to quantify the sensitivity of these signals to two phenomenological modifications of general relativity, namely a putative gravitational dipole emission and a non-zero mass for the graviton, properly accounting for the detector's response. Moreover, we consider a scenario where those sources could be observed also with Earth-based detectors, which should measure the coalescence time with precision better than $1 \ {\rm ms}$. This constraint on the coalescence time further improves the bounds that we can set on those phenomenological deviations from general relativity. We show that tests of dipole radiation and the graviton's mass should improve respectively by seven and half an order(s) of magnitude over current bounds. Finally, we discuss under which conditions one may claim the detection of a modification to general relativity.

The mechanisms leading to the formation of disks around young stellar objects (YSOs) and to the launching of the associated jets are crucial to the understanding of the earliest stages of star and planet formation. HH 212 is a privileged laboratory to study a pristine jet-disk system. Therefore we investigate the innermost region (<100 AU) around the HH 212-MM1 protostar through ALMA band\,7 observations of methanol. The 8 GHz bandwidth spectrum towards the peak of the continuum emission of the HH 212 system reveals at least 19 transitions of methanol. Several of these lines (among which several vibrationally excited lines in the v$_{\rm t}=1,2$ states) have upper energies above 500 K. They originate from a compact (<135 AU in diameter), hot ($\sim 295$ K) region elongated along the direction of the SiO jet. We performed a fit in the $uv$ plane of various velocity channels of the strongest high-excitation lines. The blue- and red-shifted velocity centroids are shifted roughly symmetrically on either side of the jet axis, indicating that the line-of-sight velocity beyond 0.7 km s$^{-1}$ from systemic is dominated by rotational motions. The velocity increases moving away from the protostar further indicating that the emission of methanol is not associated with a Keplerian disk or rotating-infalling cavity, and it is more likely associated with outflowing gas. We speculate that CH$_3$OH traces a disk wind gas accelerated at the base. The launching region would be at a radius of a few astronomical units from the YSO.

The Eigenstate Thermalization Hypothesis (ETH) implies a form for the matrix elements of local operators between eigenstates of the Hamiltonian, expected to be valid for chaotic systems. Another signal of chaos is a positive Lyapunov exponent, defined on the basis of Loschmidt echo or out-of-time-order correlators. For this exponent to be positive, correlations between matrix elements unrelated by symmetry, usually neglected, have to exist. The same is true for the peak of the dynamic heterogeneity length, relevant for systems with slow dynamics. These correlations, as well as those between elements of different operators, are encompassed in a generalized form of ETH.

We used newly obtained observations from the Mid-Infrared Instrument (MIRI) equipped with the Medium Resolution Spectrometer (MRS) aboard the James Webb Space Telescope (JWST) to extract spectra covering the entire spectral range from 5 to 27~$\mu$m in the CND and in the CC. We used the photoionization code CLOUDY to generate synthetic spectra with the same spectral range and resolution, simulating a wide range of gas phases and abundances. We then determined the contribution of each phase to the spectra. Once the abundances and contribution from each phase of the gas were determined, we identified four dominant phases and performed a spatial analysis to determine their contribution to each spaxel of the datacubes. We find that in both the CND and the CC, the bulk of the emission originates from warm ionized gas with temperatures of between $10^4$ and $10^{4.8}$~K. In the CND, molecular gas contributes significantly to the flux and is spatially structured, while the CC shows minimal molecular gas content, as is expected from these regions. Coronal gas is detected in both regions at the interface between molecular and warm ionized gas. The observed abundance pattern (enhanced CNO and $\alpha$ elements with suppressed Fe) indicates a chemically young environment, recently enriched by core-collapse supernovae and stellar winds, with a limited contribution from older Type Ia supernovae. This favors a scenario of massive, recent star formation rather than cumulative long-term enrichment.

As robots leave the controlled environments of factories to autonomously function in more complex, natural environments, they will have to respond to the inevitable fact that they will become damaged. However, while animals can quickly adapt to a wide variety of injuries, current robots cannot "think outside the box" to find a compensatory behavior when damaged: they are limited to their pre-specified self-sensing abilities, can diagnose only anticipated failure modes, and require a pre-programmed contingency plan for every type of potential damage, an impracticality for complex robots. Here we introduce an intelligent trial and error algorithm that allows robots to adapt to damage in less than two minutes, without requiring self-diagnosis or pre-specified contingency plans. Before deployment, a robot exploits a novel algorithm to create a detailed map of the space of high-performing behaviors: This map represents the robot's intuitions about what behaviors it can perform and their value. If the robot is damaged, it uses these intuitions to guide a trial-and-error learning algorithm that conducts intelligent experiments to rapidly discover a compensatory behavior that works in spite of the damage. Experiments reveal successful adaptations for a legged robot injured in five different ways, including damaged, broken, and missing legs, and for a robotic arm with joints broken in 14 different ways. This new technique will enable more robust, effective, autonomous robots, and suggests principles that animals may use to adapt to injury.

Very faint X-ray binaries appear to be transient in many cases with peak luminosities much fainter than that of usual soft X-ray transients, but their nature still remains elusive. We investigate the possibility that this transient behaviour is due to the same thermal/viscous instability which is responsible for outbursts of bright soft X-ray transients, occurring in ultracompact binaries for adequately low mass-transfer rates. More generally, we investigate the observational consequences of this instability when it occurs in ultracompact binaries. We use our code for modelling the thermal-viscous instability of the accretion disc, assumed here to be hydrogen poor. We also take into account the effects of disc X-ray irradiation, and consider the impact of the mass-transfer rate on the outburst brightness. We find that one can reproduce the observed properties of both the very faint and the brighter short transients (peak luminosity, duration, recurrence times), provided that the viscosity parameter in quiescence is slightly smaller (typically a factor of between two and four) than in bright soft X-ray transients and normal dwarf nova outbursts, the viscosity in outburst being unchanged. This possibly reflects the impact of chemical composition on non-ideal MHD effects affecting magnetically driven turbulence in poorly ionized discs.

We present significant evidence that the powerful property of Yangian invariance extends to a new large class of conformally invariant Feynman integrals. Our results apply to planar Feynman diagrams in any spacetime dimension dual to an arbitrary network of intersecting straight lines on a plane (Baxter lattice), with propagator powers determined by the geometry. We formulate Yangian symmetry in terms of a chain of Lax operators acting on the fixed coordinates around the graph, and we also extend this construction to the case of infinite-dimensional auxiliary space. Yangian invariance leads to new differential and integral equations for individual, highly nontrivial, Feynman graphs, and we present them explicitly for several examples. The graphs we consider determine correlators in the recently proposed loom fishnet CFTs. We also describe a generalization to the case with interaction vertices inside open faces of the diagram. Our construction unifies and greatly extends the known special cases of Yangian invariance to likely the most general family of integrable scalar planar graphs.

We compute the phase diagram of the simplest holographic bottom-up model of conformal interfaces. The model consists of a thin domain wall between three-dimensional Anti-de Sitter (AdS) vacua, anchored on a boundary circle. We distinguish five phases depending on the existence of a black hole, the intersection of its horizon with the wall, and the fate of inertial observers. We show that, like the Hawking-Page phase transition, the capture of the wall by the horizon is also a first order transition and comment on its field-theory interpretation. The static solutions of the domain-wall equations include gravitational avatars of the Faraday cage, black holes with negative specific heat, and an intriguing phenomenon of suspended vacuum bubbles corresponding to an exotic interface/anti-interface fusion. Part of our analysis overlaps with recent work by Simidzija and Van Raamsdonk but the interpretation is different.

We compute exactly various 4-point correlation functions of shortest scalar operators in bi-scalar planar four-dimensional "fishnet" CFT. We apply the OPE to extract from these functions the exact expressions for the scaling dimensions and the structure constants of all exchanged operators with an arbitrary Lorentz spin. In particular, we determine the conformal data of the simplest unprotected two-magnon operator analogous to the Konishi operator, as well as of the one-magnon operator. We show that at weak coupling 4-point correlation functions can be systematically expanded in terms of harmonic polylogarithm functions and verify our results by explicit calculation of Feynman graphs at a few orders in the coupling. At strong coupling we obtain that the correlation functions exhibit the scaling behaviour typical for semiclassical description hinting at the existence of the holographic dual.

We study the second-order phase transition in the $d$-dimensional Ising model with long-range interactions decreasing as a power of the distance $1/r^{d+s}$. For $s$ below some known value $s_*$, the transition is described by a conformal field theory without a local stress tensor operator, with critical exponents varying continuously as functions of $s$. At $s=s_*$, the phase transition crosses over to the short-range universality class. While the location $s_*$ of this crossover has been known for 40 years, its physics has not been fully understood, the main difficulty being that the standard description of the long-range critical point is strongly coupled at the crossover. In this paper we propose another field-theoretic description which, on the contrary, is weakly coupled near the crossover. We use this description to clarify the nature of the crossover and make predictions about the critical exponents. That the same long-range critical point can be reached from two different UV descriptions provides a new example of infrared duality.

Herschel/HIFI observations toward the compact HII region W51 has revealed the presence of a cold dense core along its line of sight in a high-velocity stream located just in front of W51. This detection has been made possible through absorption measurements of low-energy transitions of HDO, NH3, and C3 against the bright background emitted by the star-forming region. We present a follow-up study of this core using the high sensitivity and high spectral resolution provided by the IRAM 30-meter telescope. We report new detections of this core in absorption for DCO+ (2-1, 3-2), H13CO+ (1-0), DNC (3-2), HN13C (1-0), p-H2CO (202-101, 303-202), and in emission for o-NH2D. We also report interferometric observation of this last species using the IRAM/NOEMA telescope, revealing the fragmented nature of the source through the detection of two cores, separated by 0.19 - 0.24 pc, with average sizes of less than 0.16 - 0.19 pc. From a non-LTE analysis, we are able to estimate the density (~2.5e4 cm-3) and temperature (~10 K) of this component, typical of what is found in dark clouds. This component (called W51-core) has the same DCO+/HCO+ ratio (0.02) as TMC-1 and a high DNC/HNC ratio (0.14). Detection of these deuterated species indicates that W51-core is similar to an early-phase low-mass star-forming region, formed from the interaction between the W51 giant molecular cloud and the high-velocity stream in front of it. The W51 complex being at about 5 kpc, these findings lead to what is the first detection of the earliest phase of low-mass star-forming region at such a large distance.

Observations reveal a `bulk flow' in the local Universe which is faster and extends to much larger scales than is expected around a typical observer in the standard $\Lambda$CDM cosmology. This is expected to result in a scale-dependent dipolar modulation of the acceleration of the expansion rate inferred from observations of objects within the bulk flow. From a maximum-likelihood analysis of the Joint Lightcurve Analysis (JLA) catalogue of Type Ia supernovae we find that the deceleration parameter, in addition to a small monopole, indeed has a much bigger dipole component aligned with the CMB dipole which falls exponentially with redshift $z$: $q_0 = q_\mathrm{m} + \vec{q}_\mathrm{d}.\hat{n}\exp(-z/S)$. The best fit to data yields $q_\mathrm{d} = -8.03$ and $S = 0.0262~(\Rightarrow d \sim 100~\mathrm{Mpc})$, rejecting isotropy ($q_\mathrm{d} = 0$) with $3.9\sigma$ statistical significance, while $q_\mathrm{m} = -0.157$ and consistent with no acceleration ($q_\mathrm{m} = 0$) at $1.4\sigma$. Thus the cosmic acceleration deduced from supernovae may be an artefact of our being non-Copernican observers, rather than evidence for a dominant component of `dark energy' in the Universe.

Based on the rate of resolved stellar origin black hole and neutron star mergers measured by LIGO and Virgo, it is expected that these detectors will also observe an unresolved Stochastic Gravitational Wave Background (SGWB) by the time they reach design sensitivity. Since the same binaries observed by LIGO and Virgo pass through the LISA mHz frequency band at an earlier stage of their orbital evolution,it is foreseen that their SGWB will also be observable by LISA with Signal to Noise Ratio (SNR) $\sim 53$. Unlike the stochastic signal from Galactic white dwarf binaries, for which a subtraction is expected to be possible by exploiting its yearly modulation (induced by the motion of the LISA constellation), the background from unresolved stellar origin black hole and neutron star binaries acts as a foreground for other stochastic signals of cosmological or astrophysical origin, which may also be present in the LISA band. Here, we employ a principal component analysis to model and extract an additional hypothetical SGWB in the LISA band, without making any a priori assumptions on its spectral shape. At the same time, we account for the presence of the foreground from stellar origin black holes and neutron stars, as well as for possible uncertainties in the LISA noise calibration. We find that our technique leads to a linear problem and is therefore suitable for fast and reliable extraction of SGWBs with SNR up to ten times weaker than the LIGO/Virgo foreground, quite independently of their spectral shape.

Aiming at increasing the yield strength of transformation and twinning induced plasticity (TRIP and TWIP) titanium alloys, a dual-phase $\alpha$/$\beta$ alloy is designed and studied. The composition Ti 7Cr 1.5Sn (wt.%) is proposed, based on an approach coupling Calphad calculations and classical Bo-Md design tool used in Ti-alloys. Its microstructure is made of 20% of $\alpha$ precipitates in a $\beta$ matrix, the matrix having optimal Bo and Md parameters for deformation twinning and martensitic transformation. The alloy indeed displays a yield strength of 760 MPa, about 200 MPa above that of a Ti 8.5Cr 1.5Sn (wt.%) single beta phase TRIP/TWIP alloy, combined with good ductility and work-hardening. In situ synchrotron X ray diffraction and post-mortem electron back-scattered analyses are performed to characterize the deformation mechanisms. They evidence that the TRIP and TWIP mechanisms are successfully obtained in the material, validating the design strategy. The interaction of the precipitates with the {332}<113> $\beta$ twins is analyzed, evidencing that the precipitates are sheared when hit by a twin, and therefore do not hinder the propagation of the twins. The detailed nature of the interaction is discussed, as well as the impact of the precipitates on the mechanical properties.

We evaluate the impact of imaging systematics on the clustering of luminous red galaxies (LRG), emission-line galaxies (ELG) and quasars (QSO) targeted for the upcoming Dark Energy Spectroscopic Instrument (DESI) survey. Using Data Release 7 of the DECam Legacy Survey, we study the effects of astrophysical foregrounds, stellar contamination, differences between north galactic cap and south galactic cap measurements, and variations in imaging depth, stellar density, galactic extinction, seeing, airmass, sky brightness, and exposure time before presenting survey masks and weights to mitigate these effects. With our sanitized samples in hand, we conduct a preliminary analysis of the clustering amplitude and evolution of the DESI main targets. From measurements of the angular correlation functions, we determine power law fits $r_0 = 7.78 \pm 0.26$$h^{-1}$Mpc,$\gamma = 1.98 \pm 0.02$for LRGs and$r_0 = 5.45 \pm 0.1$$h^{-1}$Mpc,$\gamma = 1.54 \pm 0.01$ for ELGs. Additionally, from the angular power spectra, we measure the linear biases and model the scale dependent biases in the weakly nonlinear regime. Both sets of clustering measurements show good agreement with survey requirements for LRGs and ELGs, attesting that these samples will enable DESI to achieve precise cosmological constraints. We also present clustering as a function of magnitude, use cross-correlations with external spectroscopy to infer $dN/dz$ and measure clustering as a function of luminosity, and probe higher order clustering statistics through counts-in-cells moments.