Dual feed progressive cavity pump extrusion system for functionally graded direct ink write 3D printing

Material extrusion Additive Manufacturing (AM), is one of the most widely practiced methods of AM. Fused Filament Fabrication (FFF) is what most associate with AM, as it is relatively inexpensive, and highly accessible, involving feeding plastic filament into a hot-end that melts and extrudes from a nozzle as the toolhead moves along the toolpath. Direct Ink Write (DIW) 3D printing falls into this same category of AM, however is primarily practiced in laboratory settings to construct novel parts from flowable feedstock materials. DIW printers are relatively expensive and often depend on custom software to print a part, limiting user-specificity. There have been recent advancements in multi-material and functionally graded DIW, but the systems are highly custom and the methods used to achieve multi-material prints are openly available to the public. The following article outlines the construction and operation method of a DIW system that is capable of printing that can produce compositionally-graded components using a dual feed progressive cavity pump extruder equipped with a dynamic mixer. The extruder and its capabilities to vary material composition while printing are demonstrated using a Prusa i3 MK3S+ desktop fused filament fabrication printer as the gantry system. This provides users ease of operation, and the capability of further tailoring to specific needs.

Suppression of Richtmyer-Meshkov Instability via Special Pairs of Shocks and Phase Transitions

The classical Richtmyer-Meshkov instability (RMI) is a hydrodynamic instability characterizing the evolution of an interface following shock loading. In contrast to other hydrodynamic instabilities such as Rayleigh-Taylor, it is known for being unconditionally unstable: regardless of the direction of shock passage, any deviations from a flat interface will be amplified. In this article, we show that for negative Atwood numbers, there exist special sequences of shocks which result in a nearly perfectly suppressed instability growth. We demonstrate this principle computationally and experimentally with stepped fliers and phase transition materials. A fascinating immediate corollary is that in specific instances, a phase-transitioning material may self-suppress RMI.

Scaling Law for the Onset of Solidification at Extreme Undercooling

Quasi-isentropic compression enables one to study the solidification of metastable liquid states that are inaccessible through other experimental means. The onset of this nonequilibrium solidification is known to depend on the compression rate and material-specific factors, but this complex interdependence has not been well characterized. In this study, we use a combination of experiments, theory, and computational simulations to derive a general scaling law that quantifies this dependence. One of its applications is a novel means to elucidate melt temperatures at high pressures.

Model for the Solid-Liquid Interfacial Free Energy at High Pressures

The free energy involved in the formation of an interface between two phases (e.g., a solid-liquid interface) is referred to as the interfacial free energy. For the case of solidification, the interfacial free energy dictates the height of the energy barrier required to nucleate stable clusters of the newly forming solid phase and is essential for producing an accurate solidification kinetics model using classical nucleation theory (CNT)-based methods. While various methods have been proposed for modeling the interfacial free energy for solid-liquid interfaces in prior literature, many of these formulations involve making restrictive assumptions or approximations, such as the system being at or near equilibrium (i.e., the system temperature is approximately equal to the melt temperature) or that the system is at pressures close to atmospheric. However, these approximations and assumptions may break down in highly non-equilibrium situations, such as in dynamic-compression experiments where metastable liquids that are undercooled by hundreds of kelvin or overpressurized by several gigapascals or more are formed before eventually solidifying. We derive a solid-liquid interfacial free-energy model for such high-pressure conditions by considering the enthalpies of interactions between pairs of atoms or molecules. We also consider the contribution of interface roughness (disordering) by incorporating a multilayer interface model known as the Temkin n-layer model. Our formulation is applicable to a diverse variety of materials, and we demonstrate it by developing models specifically for two different materials: water and gallium. We apply our interfacial free-energy formulation to CNT-based kinetics simulations of several suites of dynamic-compression experiments that cause liquid water to solidify to the high-pressure solid polymorph ice VII and have found good agreement to the observed kinetics with only minor empirical fitting.

Apparatus for measuring strength in biaxial compression

Most measurements of compressive strength of ductile materials have involved Hopkinson-Kolsky bars or Taylor anvils placing samples in uniaxial compression. In these geometries, strain is limited by the tendency of the sample to petal, in analogy to necking in uniaxial tension. Estimation of strength for any other form of the stress tensor requires assuming a shape of the yield surface; because data exist only for uniaxial compression, these assumptions are untested. In an imploding spherical shell, compression is biaxial, the plastic strain may not be small, and the material behavior may be nonlinear as a result of work hardening and heating by plastic work. We outline a method of measuring the strengths of materials in biaxial compression, both quasistatically and dynamically, using the compression of thin spherical shells. We suggest surrounding the shell with an annulus filled with a mixture of H₂ and Cl₂ gases whose homogeneous ignition is initiated by a flash of blue and near-ultraviolet light. Less promising approaches are described in Appendixes A-C.

Measuring the melting curve of iron at super-Earth core conditions

The discovery of more than 4500 extrasolar planets has created a need for modeling their interior structure and dynamics. Given the prominence of iron in planetary interiors, we require accurate and precise physical properties at extreme pressure and temperature. A first-order property of iron is its melting point, which is still debated for the conditions of Earth’s interior. We used high-energy lasers at the National Ignition Facility and in situ x-ray diffraction to determine the melting point of iron up to 1000 gigapascals, three times the pressure of Earth’s inner core. We used this melting curve to determine the length of dynamo action during core solidification to the hexagonal close-packed (hcp) structure. We find that terrestrial exoplanets with four to six times Earth’s mass have the longest dynamos, which provide important shielding against cosmic radiation.

Metastability of Liquid Water Freezing into Ice VII under Dynamic Compression

The ubiquitous nature and unusual properties of water have motivated many studies on its metastability under temperature- or pressure-induced phase transformations. Here, nanosecond compression by a high-power laser is used to create the nonequilibrium conditions where liquid water persists well into the stable region of ice VII. Through our experiments, as well as a complementary theoretical-computational analysis based on classical nucleation theory, we report that the metastability limit of liquid water under nearly isentropic compression from ambient conditions is at least 8 GPa, higher than the 7 GPa previously reported for lower loading rates.

Minimization of Gibbs Energy in High-Pressure Multiphase, Multicomponent Mixtures through Particle Swarm Optimization

We present a global optimization method to construct phase boundaries in multicomponent mixtures by minimizing the Gibbs energy. The minimization method is, in essence, an extension of the Maxwell construction procedure that is used in single-component systems. For a given temperature, pressure, and overall mixture composition, it reveals the mole fractions of the thermodynamically stable phases and the composition of these phases. Our approach is based on particle swarm optimization (PSO), which is a gradient-free, stochastic method. It is not reliant on good initial guesses for the phase fractions and compositions, which is an important requirement for the high-pressure applications considered in this study because data on phase boundaries at high pressures tend to be extremely limited. One practical use of this method is to create equation-of-state tables needed by continuum-scale, multiphysics codes that are ubiquitous in high-pressure science. Currently, there does not exist a method to generate such tables that rigorously account for changes in phase boundaries due to mixing. We have done extensive testing to demonstrate that PSO can reliably determine the Gibbs energy minimum and can capture nontrivial features like eutectic and peritectic temperatures to produce coherent phase diagrams. As part of our testing, we have developed a PSO-based Helmholtz-energy minimization procedure that we have used to cross-check the results of the Gibbs energy minimization. We conclude with a critique of our approach and provide suggestions for future work, including a PSO-based entropy-maximization method that would enable the aforementioned continuum codes to perform on-the-fly, phase-equilibria calculations of multicomponent mixtures.

Numerical modeling of solid-cluster evolution applied to the nanosecond solidification of water near the metastable limit

Classical nucleation theory (CNT) is a promising way to predictively model the submicrosecond kinetics of phase transitions that occur under dynamic compression, such as the suite of experiments performed over the past two decades on the solidification of liquid water to the high-pressure ice VII phase. Myint et al. [Phys. Rev. Lett. 121, 155701 (2018)] presented the first CNT-based model for these types of rapid phase transitions, but relied on an empirical scaling parameter in their transient induction model to simulate the lag time that occurs prior to the onset of significant formation of ice VII clusters in the system. To build on that study, we model the liquid water-ice VII phase transformation using a numerical discretization scheme to solve the Zel'dovich-Frenkel partial differential equation, which is a fundamental CNT-based kinetic equation that describes the statistical time-dependent behavior of solid cluster formation. The Zel'dovich-Frenkel equation inherently accounts for transience in the nucleation kinetics and eliminates the need for the empirical scaling factor used by Myint et al. One major result of this research is that transience is found to play a relatively small role in the nucleation process for the dynamic-compression time scales of the liquid water-ice VII experiments being simulated. Instead, we show that it is possible to properly model the lag time using steady-state CNT by making small refinements to the interfacial free energy value. We have also developed a new dimensionless parameter that may be applied a priori to predict whether or not transient nucleation will be important in a given dynamic-compression experiment.

Nanosecond Freezing of Water at High Pressures: Nucleation and Growth near the Metastability Limit

The fundamental study of phase transition kinetics has motivated experimental methods toward achieving the largest degree of undercooling possible, more recently culminating in the technique of rapid, quasi-isentropic compression. This approach has been demonstrated to freeze water into the high-pressure ice VII phase on nanosecond timescales, with some experiments undergoing heterogeneous nucleation while others, in apparent contradiction, suggest a homogeneous nucleation mode. In this study, we show through a combination of theory, simulation, and analysis of experiments that these seemingly contradictory results are in agreement when viewed from the perspective of classical nucleation theory. We find that, perhaps surprisingly, classical nucleation theory is capable of accurately predicting the solidification kinetics of ice VII formation under an extremely high driving force (|Δμ/k_{B}T|≈1) but only if amended by two important considerations: (i) transient nucleation and (ii) separate liquid and solid temperatures. This is the first demonstration of a model that is able to reproduce the experimentally observed rapid freezing kinetics.

The thermodynamics of a liquid-solid interface at extreme conditions: A model close-packed system up to 100 GPa

The first experimental insight into the nature of the liquid-solid interface occurred with the pioneering experiments of Turnbull, which simultaneously demonstrated both that metals could be deeply undercooled (and therefore had relatively large barriers to nucleation) and that the inferred interfacial free energy γ was linearly proportional to the enthalpy of fusion [D. Turnbull, J. Appl. Phys. 21, 1022 (1950)]. By an atomistic simulation of a model face-centered cubic system via adiabatic free energy dynamics, we extend Turnbull's result to the realm of high pressure and demonstrate that the interfacial free energy, evaluated along the melting curve, remains linear with the bulk enthalpy of fusion, even up to 100 GPa. This linear dependence of γ on pressure is shown to be a consequence of the entropy dominating the free energy of the interface in conjunction with the fact that the entropy of fusion does not vary greatly along the melting curve for simple monoatomic metals. Based on this observation, it appears that large undercoolings in liquid metals can be achieved even at very high pressure. Therefore, nucleation rates at high pressure are expected to be non-negligible, resulting in observable solidification kinetics.

Rapid freezing of water under dynamic compression

Understanding the behavior of materials at extreme pressures is a central issue in fields like aerodynamics, astronomy, and geology, as well as for advancing technological grand challenges such as inertial confinement fusion. Dynamic compression experiments to probe high-pressure states often encounter rapid phase transitions that may cause the materials to behave in unexpected ways, and understanding the kinetics of these phase transitions remains an area of great interest. In this review, we examine experimental and theoretical/computational efforts to study the freezing kinetics of water to a high-pressure solid phase known as ice VII. We first present a detailed analysis of dynamic compression experiments in which water has been observed to freeze on sub-microsecond time scales to ice VII. This is followed by a discussion of the limitations of currently available molecular and continuum simulation methods in modeling these experiments. We then describe how our phase transition kinetics models, which are based on classical nucleation theory, provide a more physics-based framework that overcomes some of these limitations. Finally, we give suggestions on future experimental and modeling work on the liquid-ice VII transition, including an outline of the development of a predictive multiscale model in which molecular and continuum simulations are intimately coupled.

Extraction of effective solid-liquid interfacial free energies for full 3D solid crystallites from equilibrium MD simulations

Molecular dynamics simulations of an embedded atom copper system in the isobaric-isenthalpic ensemble are used to study the effective solid-liquid interfacial free energy of quasi-spherical solid crystals within a liquid. This is within the larger context of molecular dynamics simulations of this system undergoing solidification, where single individually prepared crystallites of different sizes grow until they reach a thermodynamically stable final state. The resulting equilibrium shapes possess the full structural details expected for solids with weakly anisotropic surface free energies (in these cases, ∼5% radial flattening and rounded [111] octahedral faces). The simplifying assumption of sphericity and perfect isotropy leads to an effective interfacial free energy as appearing in the Gibbs-Thomson equation, which we determine to be ∼177 erg/cm², roughly independent of crystal size for radii in the 50-250 Å range. This quantity may be used in atomistically informed models of solidification kinetics for this system.

Grain-size-independent plastic flow at ultrahigh pressures and strain rates

A basic tenet of material science is that the flow stress of a metal increases as its grain size decreases, an effect described by the Hall-Petch relation. This relation is used extensively in material design to optimize the hardness, durability, survivability, and ductility of structural metals. This Letter reports experimental results in a new regime of high pressures and strain rates that challenge this basic tenet of mechanical metallurgy. We report measurements of the plastic flow of the model body-centered-cubic metal tantalum made under conditions of high pressure (>100  GPa) and strain rate (∼10(7)  s(-1)) achieved by using the Omega laser. Under these unique plastic deformation ("flow") conditions, the effect of grain size is found to be negligible for grain sizes >0.25  μm sizes. A multiscale model of the plastic flow suggests that pressure and strain rate hardening dominate over the grain-size effects. Theoretical estimates, based on grain compatibility and geometrically necessary dislocations, corroborate this conclusion.

Fabrication and application of high impedance graded density impactors in light gas gun experiments

Recent advances in Graded Density Impactor fabrication technique have increased the maximum achievable pressure in gas gun quasi-isentropic experiments to 5 Mbars. In this report, we outline the latest methodologies and applications of Graded Density Impactors in experiments at extreme conditions. These new Graded Density Impactors are essentially metallic discs made of nearly one hundred layers of precisely mixed Mg, Cu, and W. The density gradients in these impactors are specifically designed to generate the desired thermodynamic path required for each experiment. We carried out a number of experiments at various pressures using these Graded Density Impactors. These experimental results and their simulations will be presented here.

A molecular H2 potential for heterogeneous simulations including polarization and many-body van der Waals interactions

A highly accurate aniostropic intermolecular potential for diatomic hydrogen has been developed that is transferable for molecular modeling in heterogeneous systems. The potential surface is designed to be efficacious in modeling mixed sorbates in metal-organic materials that include sorption interactions with charged interfaces and open metal sites. The potential parameters are compatible for mixed simulations but still maintain high accuracy while deriving dispersion parameters from a proven polarizability model. The potential includes essential physical interactions including: short-range repulsions, dispersion, and permanent and induced electrostatics. Many-body polarization is introduced via a point-atomic polarizability model that is also extended to account for many-body van der Waals interactions in a consistent fashion. Permanent electrostatics are incorporated using point partial charges on atomic sites. However, contrary to expectation, the best potentials are obtained by permitting the charges to take on values that do not reproduce the first non-vanishing moment of the electrostatic potential surface, i.e., the quadrupole moment. Potential parameters are fit to match ab initio energies for a representative range of dimer geometries. The resulting potential is shown to be highly effective by comparing to electronic structure calculations for a thermal distribution of trimer geometries, and by reproducing experimental bulk pressure-density isotherms. The surface is shown to be superior to other similarly portable potential choices even in tests on homogeneous systems without strong polarizing fields. The present streamlined approach to developing such potentials allows for a simple adaptation to other molecules amenable to investigation by high-level electronic structure methods.

Characterization of Tunable Radical Metal-Carbenes: Key Intermediates in Catalytic Cyclopropanation

A new class of radical metal-carbene complex has been characterized as having Fischer-like orbital interactions and adjacent π acceptor stabilization. Density Functional Theory (DFT) along with Natural Bond Orbital (NBO) analysis and Charge Decomposition Analysis (CDA) has given insight into the electronics of this catalytic intermediate in an open-shell cobalt-porphyrin, [Co(Por)], system. The complex has a single bond from the metal to the carbene and has radical character with localized spin density on the carbene carbon. In addition, the carbene carbon is found to be nucleophilic and "tunable" through the introduction of different α-carbon substituents. Finally, based on these findings, rational design strategies are proposed which should lead to the enhancement of catalytic activity.

Evidence for substrate preorganization in the peptidylglycine α-amidating monooxygenase reaction describing the contribution of ground state structure to hydrogen tunneling

Peptidylglycine α-amidating monooxygenase (PAM) is a bifunctional enzyme which catalyzes the post-translational modification of inactive C-terminal glycine-extended peptide precursors to the corresponding bioactive α-amidated peptide hormone. This conversion involves two sequential reactions both of which are catalyzed by the separate catalytic domains of PAM. The first step, the copper-, ascorbate-, and O(2)-dependent stereospecific hydroxylation at the α-carbon of the C-terminal glycine, is catalyzed by peptidylglycine α-hydroxylating monooxygenase (PHM). The second step, the zinc-dependent dealkylation of the carbinolamide intermediate, is catalyzed by peptidylglycine amidoglycolate lyase. Quantum mechanical tunneling dominates PHM-dependent C(α)-H bond activation. This study probes the substrate structure dependence of this chemistry using a set of N-acylglycine substrates of varying hydrophobicity. Primary deuterium kinetic isotope effects (KIEs), molecular mechanical docking, alchemical free energy perturbation, and equilibrium molecular dynamics were used to study the role played by ground-state substrate structure on PHM catalysis. Our data show that all Ν-acylglycines bind sequentially to PHM in an equilibrium-ordered fashion. The primary deuterium KIE displays a linear decrease with respect to acyl chain length for straight-chain N-acylglycine substrates. Docking orientation of these substrates displayed increased dissociation energy proportional to hydrophobic pocket interaction. The decrease in KIE with hydrophobicity was attributed to a preorganization event which decreased reorganization energy by decreasing the conformational sampling associated with ground state substrate binding. This is the first example of preorganization in the family of noncoupled copper monooxygenases.

Glucosamine-induced increase in Akt phosphorylation corresponds to increased endoplasmic reticulum stress in astroglial cells

Increased glucose flux through the hexosamine biosynthetic pathway (HBP) is known to affect the activity of a number of signal transduction pathways and lead to insulin resistance. Although widely studied in insulin responsive tissues, the effect of increased HBP activity on largely insulin unresponsive tissues, such as the brain, remains relatively unknown. Herein, we investigate the effects of increased HBP flux on Akt activation in a human astroglial cells line using glucosamine, a compound commonly used to mimic hyperglycemic conditions by increasing HBP flux. Cellular treatment with 8 mM glucosamine resulted in a 96.8% +/- 24.6 increase in Akt phosphorylation after 5 h of treatment that remained elevated throughout the 9-h time course. Glucosamine treatment also resulted in modest increases in global levels of the O-GlcNAc protein modification. Increasing O-GlcNAc levels using the O-GlcNAcase inhibitor streptozotocin (STZ) also increased Akt phosphorylation by 96.8% +/- 11.0 after only 3 h although for a shorter duration than glucosamine; however, the more potent O-GlcNAcase inhibitors O-(2-acetamido-2-deoxy-D-glucopyranosylidene)amino-N-phenylcarbamate (PUGNAc) and 1,2-dideoxy-2'-propyl-alpha-D-glucopyranoso-[2,1-d]-Delta2'-thiazoline (NAGBT) failed to mimic the increases in phospho-Akt indicating that the Akt phosphorylation is not a result of increased O-GlcNAc protein modification. Further analysis indicated that this increased phosphorylation was also not due to increased osmotic stress and was not attenuated by N-acetylcysteine eliminating the potential role of oxidative stress in the observed phospho-Akt increases. Glucosamine treatment, but not STZ treatment, did correlate with a large increase in the expression of the endoplasmic reticulum (ER) stress marker GRP 78. Altogether, these results indicate that increased HBP flux in human astroglial cells results in a rapid, short-term phosphorylation of Akt that is likely a result of increased ER stress. The mechanism by which STZ increases Akt phosphorylation, however, remains unknown.