Search

Dualities play a central role in the study of quantum spin chains, providing insight into the structure of quantum phase diagrams and phase transitions. In this work we study categorical dualities, which are defined as bounded-spread isomorphisms between algebras of symmetry-respecting local operators on a spin chain. We consider generalized global symmetries that correspond to unitary fusion categories, which are represented by matrix-product operator algebras. A fundamental question about dualities is whether they can be extended to quantum cellular automata on the larger algebra generated by all local operators that respect the unit matrix-product operator. For conventional global symmetries, which are on-site representations of finite groups, this larger algebra is simply the tensor product of algebras associated to individual spins in the chain. We present a solution to the extension problem using the machinery of Doplicher-Haag-Roberts bimodules. Our solution provides a crisp categorical criterion for when an extension of a duality exists. We show that the set of possible extensions form a torsor over the invertible objects in the relevant symmetry category. Our approach recovers existing results from the literature when applied to quantum cellular automata that respect a conventional global symmetry., Comment: Comments welcome!

Significant success has been achieved in fabricating angstrom-scale artificial solid ionic channels aiming to replicate the biological ion channels (BICs).Besides high selectivity, BICs also exhibit sophisticated ion gating and interplay. However, such behavior and functionality are seldomly recreated in the artificial counterparts due to the insufficient understanding of the molecular origin. Here we report cooperative and inhibitory ion transport in angstrom-scale acetate functionalized MoS2 two dimensional channels., Comment: 34 pages for main text and 40 pages for SI

This white paper is the result of a collaboration by those that attended a workshop at the Facility for Rare Isotope Beams (FRIB), organized by the FRIB Theory Alliance (FRIB-TA), on Theoretical Justifications and Motivations for Early High-Profile FRIB Experiments. It covers a wide range of topics related to the science that will be explored at FRIB. After a brief introduction, the sections address: (II) Overview of theoretical methods, (III) Experimental capabilities, (IV) Structure, (V) Near-threshold Physics, (VI) Reaction mechanisms, (VII) Nuclear equations of state, (VIII) Nuclear astrophysics, (IX) Fundamental symmetries, and (X) Experimental design and uncertainty quantification., Comment: 227 pages, 24 figures

We report how the unique temporal and spectral features of pulsed entangled photons from a parametric downconversion source can be utilized for manipulating electronic excitations through the optimization of their spectral phase. A new comprehensive optimization protocol based on Bayesian optimization has been developed in this work to selectively excite electronic states accessible by two-photon absorption. Using our optimization method, the entangled two-photon absorption probability for a thiophene dendrimer can be enhanced by up to a factor of 20 while classical light turns out to be nonoptimizable. Moreover, the optimization involving photon entanglement enables selective excitation that would not be possible otherwise. In addition to optimization, we have explored entangled two-photon absorption in the small entanglement time limit showing that entangled light can excite molecular electronic states that are vanishingly small for classical light. We demonstrate these opportunities with an application to a thiophene dendrimer.

Rare gynecological tumors (RGTs) present major clinical challenges due to their low incidence and heterogeneity. The lack of clear guidelines leads to suboptimal management and poor prognosis. Molecular tumor boards accelerate access to effective therapies by tailoring treatment based on biomarkers, beyond cancer type. Unstructured data that requires manual curation hinders efficient use of biomarker profiling for therapy matching. This study explores the use of large language models (LLMs) to construct digital twins for precision medicine in RGTs. Our proof-of-concept digital twin system integrates clinical and biomarker data from institutional and published cases (n=21) and literature-derived data (n=655 publications with n=404,265 patients) to create tailored treatment plans for metastatic uterine carcinosarcoma, identifying options potentially missed by traditional, single-source analysis. LLM-enabled digital twins efficiently model individual patient trajectories. Shifting to a biology-based rather than organ-based tumor definition enables personalized care that could advance RGT management and thus enhance patient outcomes., Comment: 20 pages, 2 figures, 3 tables, supplements, original article

At the low temperature and high density conditions of a neutron star crust neutrons are degenerate. In this work, we study the effect of this degeneracy on the capture rates of neutrons on neutron rich nuclei in accreted crusts. We use a statistical Hauser-Feshbach model to calculate neutron capture rates and find that neutron degeneracy can increase rates significantly. Changes increase from a factor of a few to many orders of magnitude near the neutron drip line. We also quantify uncertainties due to model inputs for masses, $\gamma$-strength functions, and level densities. We find that uncertainties increase dramatically away from stability and that degeneracy tends to increase these uncertainties further, except for cases near the neutron drip line where degeneracy leads to more robustness. As in the case of capture of classically distributed neutrons, variations in the mass model have the strongest impact. Corresponding variations in the reaction rates can be as high as 3 to 4 orders of magnitude, and be more than 5 times larger than under classical conditions. To ease the incorporation of neutron degeneracy in nucleosynthesis networks, we provide tabulated results of capture rates as well as analytical expressions as function of temperature and neutron chemical potential, for proton numbers between $3 \le Z \le 85$, derived from fits to our numerical results. Fits are based on a new parametrization that complements previously employed power law approximations with additional Lorentzian terms that account for low energy resonances, significantly improving accuracy., Comment: 16 pages, 11 figures. Tabulated reaction rates can be found at https://github.com/Nyrb1001/Degenerate-Neutron-Captures

We elaborate an algebraic framework for describing internal topological symmetries of gapped boundaries of (2+1)D topological orders. We present a categorical obstruction to the coherence of bulk group symmetry and boundary symmetries in terms of liftings of categorical actions on the bulk theory to a certain 2-group of boundary symmetries., Comment: 18 pages, 11 figures. Comments welcome

As a first step towards a complete computational model of speech learning involving perception-production loops, we investigate the forward mapping between pseudo-motor commands and articulatory trajectories. Two phonological feature sets, based respectively on generative and articulatory phonology, are used to encode a phonetic target sequence. Different interpolation techniques are compared to generate smooth trajectories in these feature spaces, with a potential optimisation of the target value and timing to capture co-articulation effects. We report the Pearson correlation between a linear projection of the generated trajectories and articulatory data derived from a multi-speaker dataset of electromagnetic articulography (EMA) recordings. A correlation of 0.67 is obtained with an extended feature set based on generative phonology and a linear interpolation technique. We discuss the implications of our results for our understanding of the dynamics of biological motion., Comment: accepted at Interspeech 2024

The high contrast and spatial resolution requirements for directly imaging exoplanets requires effective coordination of wavefront control, coronagraphy, observation techniques, and post-processing algorithms. However, even with this suite of tools, identifying and retrieving exoplanet signals embedded in resolved scattered light regions can be extremely challenging due to the increased noise from scattered light off the circumstellar disk and the potential misinterpretation of the true nature of the detected signal. This issue pertains not only to imaging terrestrial planets in habitable zones within zodiacal and exozodiacal emission but also to young planets embedded in circumstellar, transitional, and debris disks. This is particularly true for H{\alpha} detection of exoplanets in transitional disks. This work delves into recent H{\alpha} observations of three transitional disks systems with MagAO-X, an extreme adaptive optics system for the 6.5-meter Magellan Clay telescope. We employed angular differential imaging (ADI) and simultaneous spectral differential imaging (SSDI) in combination with KLIP, a PCA algorithm in post-processing, for optimal starlight suppression and quasi-static noise removal. We discuss the challenges in protoplanet identification with MagAO-X in environments rich with scattered and reflected light from disk structures and explore a potential solution for removing noise contributions from real astronomical objects with current observation and post-processing techniques.

The next generation of extreme adaptive optics (AO) must be calibrated exceptionally well to achieve the desired contrast for ground-based direct imaging exoplanet targets. Current wavefront sensing and control system responses deviate from lab calibration throughout the night due to non linearities in the wavefront sensor (WFS) and signal loss. One cause of these changes is the optical gain (OG) effect, which shows that the difference between actual and reconstructed wavefronts is sensitive to residual wavefront errors from partially corrected turbulence. This work details on-sky measurement of optical gain on MagAO-X, an extreme AO system on the Magellan Clay 6.5m. We ultimately plan on using a method of high-temporal frequency probes on our deformable mirror to track optical gain on the Pyramid WFS. The high-temporal frequency probes, used to create PSF copies at 10-22 lambda /D, are already routinely used by our system for coronagraph centering and post-observation calibration. This method is supported by the OG measurements from the modal response, measured simultaneously by sequenced pokes of each mode. When tracked with DIMM measurements, optical gain calibrations show a clear dependence on Strehl Ratio, and this relationship is discussed. This more accurate method of calibration is a crucial next step in enabling higher fidelity correction and post processing techniques for direct imaging ground based systems., Comment: 7 pages, 6 figures, proceedings of Adaptive Optics Systems IX at SPIE Astronomical Telescopes + Instrumentation 2024

High-contrast imaging data analysis depends on removing residual starlight from the host star to reveal planets and disks. Most observers do this with principal components analysis (i.e. KLIP) using modes computed from the science images themselves. These modes may not be orthogonal to planet and disk signals, leading to over-subtraction. The wavefront sensor data recorded during the observation provide an independent signal with which to predict the instrument point-spread function (PSF). MagAO-X is an extreme adaptive optics (ExAO) system for the 6.5-meter Magellan Clay telescope and a technology pathfinder for ExAO with GMagAO-X on the upcoming Giant Magellan Telescope. MagAO-X is designed to save all sensor information, including kHz-speed wavefront measurements. Our software and compressed data formats were designed to record the millions of training samples required for machine learning with high throughput. The large volume of image and sensor data lets us learn a PSF model incorporating all the information available. This will eventually allow us to probe smaller star-planet separations at greater sensitivities, which will be needed for rocky planet imaging., Comment: 9 pages, 5 figures, proceedings of Adaptive Optics Systems IX at SPIE Astronomical Telescopes + Instrumentation 2024

In this study, we investigate second- and third-harmonic generation processes in Au nanorod systems using the real-time time-dependent density functional tight binding method. Our study focuses on the computation of nonlinear signals based on the time dependent dipole response induced by linearly polarized laser pulses interacting with nanoparticles. We systematically explore the influence of various laser parameters, including pump intensity, duration, frequency, and polarization directions, on harmonic generation. We demonstrate all the results using Au nanorod dimer systems arranged in end-to-end configurations, and disrupting the spatial symmetry of regular single nanorod systems is crucial for second-harmonic generation processes. Furthermore, we study the impact of nanorod lengths, which lead to variable plasmon energies, on harmonic generation, and estimates of polarizabilities and hyper-polarizabilities are provided.

Excited atomic nitrogen atoms play an important role in plasma formation in hypersonic shock-waves, as happens during spacecraft reentry and other high velocity vehicle applications. In this study, we have thoroughly studied collision induced excitation (CIE) associated with two colliding nitrogen atoms in the N(4S), N(2D) and N(2P) states at collisions energies up to 6 eV, using time-independent scattering calculations to determine cross sections and temperature-dependent rate coefficients. The calculations are based on potential curves and couplings determined in earlier MRCI calculations with large basis sets, and the results are in good agreement with experiment where comparisons are possible. To properly consider the spin-orbit coupling matrix, we have developed a scaling method for treating transitions between different fine-structure components with calculations that only require calculations with two coupled states, and with this we define accurate degeneracy factors for determining cross sections and rate coefficients that include all states. The results indicate that both spin-orbit and derivative coupling effects can play important roles in collisional excitation and quenching, and that although derivative coupling is always much stronger than spin-orbit, there are many transitions where only spin-orbit can contribute. As part of this, we identify two distinct pathways associated with N(2D) relaxation, including one Auger-like mechanism leading to 2N(2D) that could be important at high temperature.

The wide array of currently available genomes display a wonderful diversity in size, composition and structure with many more to come thanks to several global biodiversity genomics initiatives starting in recent years. However, sequencing of genomes, even with all the recent advances, can still be challenging for both technical (e.g. small physical size, contaminated samples, or access to appropriate sequencing platforms) and biological reasons (e.g. germline restricted DNA, variable ploidy levels, sex chromosomes, or very large genomes). In recent years, k-mer-based techniques have become popular to overcome some of these challenges. They are based on the simple process of dividing the analysed sequences (e.g. raw reads or genomes) into a set of sub-sequences of length k, called k-mers. Despite this apparent simplicity, k-mer-based analysis allows for a rapid and intuitive assessment of complex sequencing datasets. Here, we provide the first comprehensive review to the theoretical properties and practical applications of k-mers in biodiversity genomics, serving as a reference manual for this powerful approach., Comment: Main text is 25 pages, 4 figures; With supplement it is 44

In-situ environmental transmission electron microscopy (E-TEM) provides sub-Angstrom spatiotemporal information on nanoscale chemical transformations. However, the high-energy electron probe can significantly affect the surrounding medium's chemistry and induce side reactions within the area of interest. Radiolysis, or the dissociation of molecules under electron beam irradiation, has become an increasingly important consideration in in-situ nanomaterials characterization. While numerous studies have explored this phenomenon in liquid-cell TEM, there are no reported extensions to gas phase TEM thus far. Herein, we present a generalized computational model to elucidate radiation chemistry in both gas and liquid E-TEM. Our model reveals that gas phase E-TEM produces radiolytic species with lower reactivity than in the liquid phase. However, these radiolysis species can accumulate to levels that interfere with desired reactions, especially at higher pressures. We verify our model by studying the radiation-promoted oxidation of aluminum nanocubes and the disproportionation of carbon monoxide, where we find increasing the electron beam dose rate increases the rate of AlOx production and deposited carbon, respectively. Additionally, we propose guidelines to control the production of radiolysis species within these closed-cell in-situ nanoreactors. Taken together, these findings should pave the way for nanoscale control of chemical reactions in E-TEM.

This study examines the Raman scattering and charge transfer properties of molecules adsorbed on the surface of a tetrahedral Au$_{120}$ nanoparticle based on the time-dependent density functional tight-binding (TD-DFTB) method. We study Raman scattering (SERS) enhancements for pyridine where the molecule is adsorbed either on the tip (V complex) or surface (S complex) of the nanoparticle. The scattering intensity is enhanced by a factor of 3-15 due to chemical effects while significantly larger enhancements (in the order of 10$^2$-10$^4$) are observed for plasmon resonance excitation at an energy of 2.5 eV depending on the adsorption site. Furthermore, we demonstrate charge transfer between the nanoparticle and a fullerene-based molecule after pulsed excitation of the plasmon resonance which shows how plasmon excitation can lead to negative molecular ion formation. All of these results are consistent with earlier studies using either TD-DFT theory or experimental measurements.

Almost all current and future high-contrast imaging instruments will use a Pyramid wavefront sensor (PWFS) as a primary or secondary wavefront sensor. The main issue with the PWFS is its nonlinear response to large phase aberrations, especially under strong atmospheric turbulence. Most instruments try to increase its linearity range by using dynamic modulation, but this leads to decreased sensitivity, most prominently for low-order modes, and makes it blind to petal-piston modes. In the push toward high-contrast imaging of fainter stars and deeper contrasts, there is a strong interest in using the PWFS in its unmodulated form. Here, we present closed-loop lab results of a nonlinear reconstructor for the unmodulated PWFS of the Magellan Adaptive Optics eXtreme (MagAO-X) system based on convolutional neural networks (CNNs). We show that our nonlinear reconstructor has a dynamic range of >600 nm root-mean-square (RMS), significantly outperforming the linear reconstructor that only has a 50 nm RMS dynamic range. The reconstructor behaves well in closed loop and can obtain >80% Strehl at 875 nm under a large variety of conditions and reaches higher Strehl ratios than the linear reconstructor under all simulated conditions. The CNN reconstructor also achieves the theoretical sensitivity limit of a PWFS, showing that it does not lose its sensitivity in exchange for dynamic range. The current CNN's computational time is 690 microseconds, which enables loop speeds of >1 kHz. On-sky tests are foreseen soon and will be important for pushing future high-contrast imaging instruments toward their limits., Comment: Accepted for publication in A&A

In this study, we numerically investigate hysteretic behaviors of secondary bifurcations in the Taylor-Couette flow of a semidilute, neutrally buoyant, and noncolloidal suspension. We consider a suspension with a bulk particle volume fraction of 0.1, a radius ratio of 0.877, and a particle size of 60. The suspension balance model (SBM) is introduced for numerical simulations to model the dynamics of particles undergoing shear-induced particle migration with rheological constitutive laws. The suspension flow transitions from circular Couette flow (CCF) via ribbons (RIB), spiral vortex flow (SVF), and wavy spiral vortex flow (WSVF) to wavy vortex flow (WVF) with the increase of suspension Reynolds number. The primary transition arises very slowly with an oscillatory critical mode and appears through a supercritical (or non-hysteretic) bifurcation. However, we find hysteretic behaviors in subsequent bifurcations (spiral vortex flow (SVF) to wavy spiral vortex flow (WSVF) and WSVF to wavy vortex flow (WVF)) during increasing-Re and decreasing-Re procedures with a rapid-step change near the transition boundaries. The WSVF and WVF states are more sustained below the transition boundaries when the Reynolds number is rapidly decreased in stages. However, the SVF and WSVF become WSVF and WVF more sharply with increasing Re, respectively. To conclude this study, we also examine in detail, a standing wave of weak counter-rotating vortices that occurs as the primary instability and analyze the wave that oscillates in time but is stationary in space in the RIB state.

The elemental abundances between strontium and silver ($Z = 38-47$) observed in the atmospheres of very metal-poor stars (VMP) in the Galaxy may contain the fingerprint of the weak $r$-process and $\nu p$-process occurring in early core-collapse supernovae explosions. In this work, we combine various astrophysical conditions based on a steady-state model to cover the richness of the supernova ejecta in terms of entropy, expansion timescale, and electron fraction. The calculated abundances based on different combinations of conditions are compared with stellar observations with the aim of constraining supernova ejecta conditions. We find that some conditions of the neutrino-driven outflows consistently reproduce the observed abundances of our sample. In addition, from the successful combinations, the neutron-rich trajectories better reproduce the observed abundances of Sr-Zr ($Z= 38-40$), while the proton-rich ones, Mo-Pd ($Z= 42-47$)., Comment: 18 pages, 11 figures, submitted to ApJ

The future quantum internet will leverage existing communication infrastructures, including deployed optical fibre networks, to enable novel applications that outperform current information technology. In this scenario, we perform a feasibility study of quantum communications over an industrial 224 km submarine optical fibre link deployed between Southport in the United Kingdom (UK) and Portrane in the Republic of Ireland (IE). With a characterisation of phase drift, polarisation stability and arrival time of entangled photons, we demonstrate the suitability of the link to enable international UK-IE quantum communications for the first time., Comment: 11 pages, 6 figures

Solution self-assembly of amphiphilic block copolymers (BCs) is typically performed by a solvent-to-water exchange. However, BC assemblies are often trapped in metastable states depending on the mixing conditions such as the magnitude and rate of water addition. BC self-assembly can be performed under near thermodynamic control by dialysis, which accounts for a slow and gradual water addition. In this Letter we report the use of a specifically designed dialysis cell to continuously monitor by dynamic light scattering and small-angle neutron scattering the morphological changes of PDMS-b-PEG BCs self-assemblies during THF-to-water exchange. The complete phase diagrams of near-equilibrium structures can then be established. Spherical micelles first form before evolving to rod-like micelles and vesicles, decreasing the total developed interfacial area of self-assembled structures in response to increasing interfacial energy as the water content increases. The dialysis kinetics can be tailored to the time scale of BC self-assembly by modifying the membrane pore size, which is of interest to study the interplay between thermodynamics and kinetics in self-assembly pathways.

Absolute cross sections measured using electromagnetic devices to separate and detect heavy recoiling ions need to be corrected for charge state fractions. Accurate prediction of charge state distributions using theoretical models is not always a possibility, especially in energy and mass regions where data is sparse. As such, it is often necessary to measure charge state fractions directly. In this paper we present a novel method of using a scintillation screen along with a CMOS camera to image the charge dispersed beam after a set of magnetic dipoles. A measurement of the charge state distribution for 88Sr passing through a natural carbon foil is performed. Using a Bayesian model to extract statistically meaningful uncertainties from these images, we find agreement between the new method and a more traditional method using Faraday cups. Future work is need to better understand systematic uncertainties. Our technique offers a viable method to measure charge state distributions.

The detection of emission lines associated with accretion processes is a direct method for studying how and where gas giant planets form, how young planets interact with their natal protoplanetary disk and how volatile delivery to their atmosphere takes place. H$\alpha$ ($\lambda=0.656\,\mu$m) is expected to be the strongest accretion line observable from the ground with adaptive optics systems, and is therefore the target of specific high-contrast imaging campaigns. We present MagAO-X and HST data obtained to search for H$\alpha$ emission from the previously detected protoplanet candidate orbiting AS209, identified through ALMA observations. No signal was detected at the location of the candidate, and we provide limits on its accretion. Our data would have detected an H$\alpha$ emission with $F_\mathrm{H\alpha}>2.5\pm0.3 \times10^{-16}$ erg s$^{-1}$ cm$^{-2}$, a factor 6.5 lower than the HST flux measured for PDS70b (Zhou et al., 2021). The flux limit indicates that if the protoplanet is currently accreting it is likely that local extinction from circumstellar and circumplanetary material strongly attenuates its emission at optical wavelengths. In addition, the data reveal the first image of the jet north of the star as expected from previous detections of forbidden lines. Finally, this work demonstrates that current ground-based observations with extreme adaptive optics systems can be more sensitive than space-based observations, paving the way to the hunt for small planets in reflected light with extremely large telescopes., Comment: 13 pages, 5 figures, accepted for publication in AJ

Rare, energetic (long) thermonuclear (Type I) X-ray bursts are classified either as intermediate-duration or superbursts, based on their duration. Intermediate-duration bursts lasting a few to tens of minutes are thought to arise from the thermonuclear runaway of a relatively thick (10^10 g/cm2) helium layer, while superbursts lasting hours are attributed to the detonation of an underlying carbon layer. We present a catalogue of 84 long thermonuclear bursts from 40 low-mass X-ray binaries, and defined from a new set of criteria distinguishing them from the more frequent short bursts. The three criteria are: (1) a total energy release larger than 10^40 erg, (2) a photospheric radius expansion phase longer than 10 s, and (3) a burst time-scale longer than 70 s. This work is based on a comprehensive systematic analysis of 70 bursts found with INTEGRAL, RXTE, Swift, BeppoSAX, MAXI, and NICER, as well as 14 long bursts from the literature that were detected with earlier generations of X-ray instruments. For each burst, we measure its peak flux and fluence, which eventually allows us to confirm the distinction between intermediate-duration bursts and superbursts. Additionally, we list 18 bursts that only partially meet the above inclusion criteria, possibly bridging the gap between normal and intermediate-duration bursts. With this catalogue, we significantly increase the number of long-duration bursts included in the MINBAR and thereby provide a substantial sample of these rare X-ray bursts for further study., Comment: Published in MNRAS, 20 pages, 7 figures, 4 tables

Achieving reliable and quantifiable performance in large-area surface-enhanced Raman spectroscopy (SERS) substrates has long been a formidable challenge. It requires substantial signal enhancement while maintaining a reproducible and uniform response. Conventional SERS substrates are typically made of inhomogeneous materials with random resonator geometries and distributions. As a result, they exhibit several or broadened plasmonic resonances, undesired absorptive losses, and inhomogeneous field enhancement. These limitations diminish the signal strength and hamper reproducibility, making it difficult to conduct comparative studies with high sensitivity. In this study, we propose an approach that utilizes monocrystalline gold flakes to fabricate well-defined gratings composed of plasmonic double-wire resonators, which are fabricated through focused ion-beam lithography. The geometry of the double wire grating substrate (DWGS) was evolutionary optimized to achieve efficient enhancement for both excitation and emission processes. The use of monocrystalline material minimizes absorption losses while enhancing the shape fidelity during the nanofabrication process. The DWGS shows notable reproducibility (RSD=6.6%), repeatability (RSD=5.6%), and large-area homogeneity over areas $>10^4\,\mu$m$^2$. Moreover, it provides a SERS enhancement factor of $\approx 10^6$ for 4-Aminothiophenol (4-ATP) analyte and detection capability for sub-monolayer coverage. The DWGS demonstrates reusability, as well as long-term stability on the shelf. Experimental validation with various analytes, in different states of matter, including biological macromolecules, confirms the sensitive and reproducible nature of DWGSs, thereby establishing them as a promising SERS substrate design for future sensing applications., Comment: 22 pages, 6 Figures, 12 pages Supplementary

The observation of X-rays during quiescence from transiently accreting neutron stars provides unique clues about the nature of dense matter. This, however, requires extensive modeling of the crusts and matching the results to observations. The pycnonuclear fusion reaction rates implemented in these models are theoretically calculated by extending phenomenological expressions and have large uncertainties spanning many orders of magnitude. We present the first sensitivity studies of these pycnonuclear fusion reactions in realistic network calculations. We also couple the reaction network with the thermal evolution code dStar to further study their impact on the neutron star cooling curves in quiescence. Varying the pycnonuclear fusion reaction rates alters the depth at which nuclear heat is deposited although the total heating remains constant. The enhancement of the pycnonuclear fusion reaction rates leads to an overall shallower deposition of nuclear heat. The impurity factors are also altered depending on the type of ashes deposited on the crust. These total changes correspond to a variation of up to 9 eV in the modeled cooling curves. While this is not sufficient to explain the shallow heat source, it is comparable to the observational uncertainties and can still be important for modeling the neutron star crust., Comment: AASTeX63, 11 pages with 9 figures