SEARCH · Engineering Papers

Results for “molecular functionalization”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 145 records · Page 8

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Bridging material models across scales: An integrated approach to equation of state and molecular dynamics modeling of copper

New uncertainty-aware equation of state (EOS) and electrical conductivity (EC) models for copper have been developed. The multiphase EOS/EC models are fit to experimental solid/liquid EC isobar measurements as well as density-functional theory molecular dynamics (DFT-MD) EC calculations in both expanded and compressed regimes (0.1–16 g/ cm 3 ⁠). The liquid and solid EOS phases were fit to available experimental data along with additional DFT-MD data over the same range as the EC. Leveraging the DFT-MD data, a corresponding machine-learned interatomic potential (MLIAP) for copper was trained using genetic-algorithm optimization. The copper MLIAP was constrained by EOS shock points at high compressions. The final EOS bounded MLIAP proves to be stable over a large density range (approximately 0.1–20 g/ cm 3 ) with good agreement to an isothermal compression curve, shock Hugoniot, and liquid speed of sound measurements at high pressures (100s of GPa).

Acoustic measurements and instrumentation↗

DOE OSTI2025

Quantum Ornstein-Zernike theory for two-temperature two-component plasmas

Laboratory plasma production almost always preferentially heats either the ions or electrons, leading to a two-temperature state. In this state, density functional theory molecular dynamic simulation is the state of the art for modeling bulk material properties. We construct a statistical mechanics model for the two temperature limit that is theoretically consistent with the molecular dynamics method. We proceed to derive the electron-ion multi-temperature quantum Ornstein-Zernike equations for the first time. This allows the construction of a two-temperature two-component plasma model using the average atom from which we can compute bulk material properties at a fraction of the computation time of the two-temperature density functional theory simulation. The accuracy of the model is benchmarked against ion pair correlation and self-diffusion results from ab initio simulation. Here, we proceed to compute the viscosity and ion thermal conductivity as a function of both ion and electron temperature.

Ab initio molecular dynamics↗

DOE OSTI2024

Experimental study of mechanistic factors influencing solvent-driven fractional crystallization of calcium sulfate

To advance dimethyl ether-driven fractional crystallization (DME-FC), a more sustainable method of water treatment and mineral recovery, a range of chemical equilibria were measured. These include varying concentrations of miscible organic solvents (MOS) used to experimentally measure the solvent-induced solid-liquid equilibrium (SLE) of calcium sulfate (CaSO 4 ) in water. Seven MOS, including dimethyl ether (DME), acetonitrile (MeCN), 1,4-dioxane, tetrahydrofuran (THF), acetone, ethanol, and diethylamine, were screened to establish trends associated with molecular volume, functional groups, and physical properties. The effect of MOS on CaSO 4 removal differed at concentrations <0.15 mol fraction MOS; MOS with greater molecular volume (THF, 1,4-dioxane, and diethylamine) induced greater CaSO 4 precipitation on a per mole basis. The solvent-induced SLE for all MOS converged between 0.15 and 0.2 mol fraction MOS, reaching a CaSO 4 concentration consistent with a water to MOS hydration ratio of 5:1 to 6:1, which may correspond to the solvent generating a solution-based pseudo-clathrate structure with continuity within the solution. Finally, solution pseudo-clathrate structures provide a mechanistic basis for DME-FC.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Calculation of self-diffusion coefficients in supercritical carbon dioxide using mean force kinetic theory

This paper presents an application of mean force kinetic theory (MFT) to the calculation of the self-diffusivity of CO 2 in the supercritical fluid regime. Two modifications to the typical application of MFT are employed to allow its application to a system of molecular species. Therefore, the first is the assumption that the inter-particle potential of mean force can be obtained from the molecule center-of-mass pair correlation function, which in the case of CO 2 is the C–C pair correlation function. The second is a new definition of the Enskog factor that describes the effect of correlations at the surface of the collision volume. The new definition retains the physical picture that this quantity represents a local density increase, resulting from particle correlations, relative to that in the zero density homogeneous fluid limit. These calculations are facilitated by the calculation of pair correlation functions from molecular dynamics (MD) simulations using the FEPM2 molecular CO2 model. The self-diffusivity calculated from theory is in good agreement with that from MD simulations up to and slightly beyond the density at the location of the Frenkel line. The calculation is compared with and is found to perform similarly well to other commonly used models but has a greater potential for application to systems of mixed species and to systems of particles with long range interatomic potentials due to electrostatic interactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Analytical derivatives of the individual state energies in ensemble density functional theory. II. Implementation on graphical processing units (GPUs)

Conical intersections control excited state reactivity, and thus, elucidating and predicting their geometric and energetic characteristics are crucial for understanding photochemistry. Locating these intersections requires accurate and efficient electronic structure methods. Unfortunately, the most accurate methods (e.g., multireference perturbation theories such as XMS-CASPT2) are computationally challenging for large molecules. The state-interaction state-averaged restricted ensemble referenced Kohn–Sham (SI-SA-REKS) method is a computationally efficient alternative. The application of SI-SA-REKS to photochemistry was previously hampered by a lack of analytical nuclear gradients and nonadiabatic coupling matrix elements. We have recently derived analytical energy derivatives for the SI-SA-REKS method and implemented the method effectively on graphical processing units. We demonstrate that our implementation gives the correct conical intersection topography and energetics for several examples. Furthermore, our implementation of SI-SA-REKS is computationally efficient, with observed sub-quadratic scaling as a function of molecular size. This demonstrates the promise of SI-SA-REKS for excited state dynamics of large molecular systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Coordinate-based simulation of pair distance distribution functions for small and large molecular assemblies: implementation and applications

X-ray scattering has become a major tool in the structural characterization of nanoscale materials. Thanks to the widely available experimental and computational atomic models, coordinate-based X-ray scattering simulation has played a crucial role in data interpretation in the past two decades. However, simulation of real-space pair distance distribution functions (PDDFs) from small- and wide-angle X-ray scattering, SAXS/WAXS, has been relatively less exploited. This study presents a comparison of PDDF simulation methods, which are applied to molecular structures that range in size from β-cyclodextrin [1 kDa molecular weight (MW), 66 non-hydrogen atoms] to the satellite tobacco mosaic virus capsid (1.1 MDa MW, 81 960 non-hydrogen atoms). The results demonstrate the power of interpretation of experimental SAXS/WAXS from the real-space view, particularly by providing a more intuitive method for understanding of partial structure contributions. Furthermore, the computational efficiency of PDDF simulation algorithms makes them attractive as approaches for the analysis of large nanoscale materials and biological assemblies. The simulation methods demonstrated in this article have been implemented in stand-alone software, SolX 3.0 , which is available to download from https://12idb.xray.aps.anl.gov/solx.html.

Zuo, Xiaobing (ORCID:0000000201344804)↗

DOE OSTI2024

Critical Role of Framework Flexibility and Disorder in Driving High Ionic Conductivity in LiNbOCl 4

Understanding Li-ion transport is key for the rational design of superionic solid electrolytes with exceptional ionic conductivities. LiNbOCl 4 is reported to be one of the most highly conducting materials in the recently realized new class of soft oxyhalide solid electrolytes, exhibiting an ionic conductivity of ~11 mS·cm -1 . Here, we apply X-ray/neutron diffraction and pair distribution function analysis - coupled with density functional theory/ab-initio molecular dynamics - to determine a structural model that provides a rationale for the high conductivity that we observe experimentally in this nanocrystalline solid. We show that it arises from unusually high framework flexibility at room temperature. This owes to isolated 1-D [NbOCl 4 ] - anionic chains which exhibit energetically favorable orientational disorder that is - in turn - correlated to multiple, disordered and equi-energetic Li + sites in the lattice. As the Li-ions sample the 3-D energy landscape with a fast predicted diffusion coefficient of 5.1 x 10 -7 cm 2 /s at room temperature (σ i calc = 17.4 mS·cm -1 ), the inorganic polymer chains can reorient or vice versa. The activation energy barrier for Li migration through the frustrated energy landscape is especially reduced by the elastic nature of the NbO 2 Cl 4 octahedra evident from very widely dispersed Cl-Nb-Cl bond angles in AIMD snapshots at 300 K. The phonon spectra are predominantly influenced by Cl vibrations in the low energy range, and there is strong overlap between the framework (Cl, Nb) and Li partial pDOS in the region between 1.2 - 4.0 THz. The framework flexibility is also reflected in a relatively low bulk modulus of 22 GPa. In conclusion, our findings pave the way for investigation of future “flex-ion” inorganic solids and open up a new direction for the design of high conductivity, soft solid electrolytes for all-solid-state batteries.

AIMD↗

DOE OSTI2020

Equation of state of hot, dense magnesium derived with first-principles computer simulations

Using two first-principles computer simulation techniques, path integral Monte Carlo and density functional theory molecular dynamics, we derive the equation of state of magnesium in the regime of warm dense matter, with densities ranging from 0.43 to 86.11 g cm –3 and temperatures from 20,000 K to 5×10 8 K. These conditions are relevant for the interiors of giant planets and stars as well as for shock compression measurements and inertial confinement fusion experiments. Here, we study ionization mechanisms and the electronic structure of magnesium as a function of density and temperature. We show that the L shell electrons, 2s and 2p energy bands, merge at high densities. This results in gradual ionization of the L-shell with increasing density and temperature. In this regard, Mg differs from MgO, which is also reflected in the shape of its principal shock Hugoniot curve. For Mg, we predict a single broad pressure-temperature region, where the shock compression ratio is approximately 4.9. Mg thus differs from Si and Al plasmas that exhibit two well-separated compression maxima on the Hugoniot curve for L and K shell ionizations. Finally, we study multiple shocks and effects of preheat and precompression.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Temperature and density dependent pair potential for deuterium under shock

Large-scale classical molecular dynamics (CMD) simulations naturally include the microscopic physics necessary for atomistic modeling of shock release at the ablator-fuel interface in an inertial confinement fusion (ICF) capsule. Here, the multi-megabar shocks utilized in ICF experiments can drive the deuterium fuel from ambient to electron volt temperatures (T) and multi-fold compression. Modeling interatomic interactions over such an extreme range of conditions is challenging for empirical bond order potentials. We generate a pair potential for deuterium with explicit temperature and mass density dependence from ab initio density functional theory molecular dynamics using the iterative Boltzmann inversion method. This potential accurately reproduces the radial distribution functions and pressures from DFT in CMD equilibrium simulations across a wide range of thermodynamic conditions, yet fails to return the expected Hugoniot relations when used in direct CMD shock simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Unified understanding of the impact of semiflexibility, concentration, and molecular weight on macromolecular-scale ring diffusion

Conformationally fluctuating, globally compact macromolecules such as polymeric rings, single-chain nanoparticles, microgels, and many-arm stars display complex dynamic behaviors due to their rich topological structure and intermolecular organization. Synthetic rings are hybrid objects with conformations that display both ideal random walk and compact globular features, which can serve as models of genomic DNA. To date, emphasis has been placed on the effect of ring molecular weight on their unusual behaviors. Here, we combine simulations and a microscopic force-level theory to build a unified understanding for how key aspects of ring dynamics depend on different tunable molecular properties including backbone rigidity, monomer concentration, degree of traditional entanglement, and molecular weight. Our large-scale molecular dynamics simulations of ring melts with very different backbone stiffnesses reveal unanticipated behaviors which agree well with our generalized theory. This includes a universal master curve for center-of-mass diffusion constants as a function of molecular weight scaled by a chemistry and thermodynamic state-dependent critical molecular weight that generalizes the concept of an entanglement cross-over for linear chains. The key physics is how backbone rigidity and monomer concentration induced changes of the entanglement length, interring packing, degree of interpenetration, and liquid compressibility slow down space-time dynamic-force correlations on macromolecular scales. A power law decay of the center-of-mass diffusion constant with inverse molecular weight squared is the first consequence, followed by an ultraslow activated hopping transport regime. Our results set the stage to address slow dynamics and kinetic arrest in different families of compact synthetic and biological polymeric systems.

Science & Technology - Other Topics↗

DOE OSTI2023

Covalent Triazine Framework-Derived Membranes: Engineered Sol–Gel Construction and Gas Separation Application

Covalent triazine frameworks (CTFs) represent one of the most extensively studied organic networks characterized by graphitic π-conjugated structures linked by aza-fused rings, possessing unique features such as compositions of light elements (e.g., C, H, and N), porous architectures abundant heteroatom involvement, and extensively conjugated structures. In addition, the textural and chemical structures of CTFs could be engineered via synthesis control to accommodate diverse applications. CTF materials with notable characteristics, including plentiful (ultra-)micropores, high surface areas, and the presence of CO 2 -philic functional groups involving nitrogen (N), oxygen (O), and fluorine (F), hold great promise as potential candidates for anthropogenic CO 2 capture and sequestration (CCS) applications. However, the conventional high-temperature involved ionothermal procedures and the solution-based coupling pathway only afforded CTF materials in powder form, which is difficult to be processed toward membrane formation. Successful fabrication of CTF-derived membranes will rely on the development of alternative polymerization approaches as well as structural engineering to afford membrane architectures with controllable porosity distribution and active interaction sites with CO 2 benefiting the CO 2 separation procedure. In this Account, a demonstration of the latest progress in the development of CTF-derived membranes was provided. The CTF membranes were mainly synthesized via a superacid (e.g., CF 3 SO 3 H)-promoted sol–gel approach involving the polymerization of aromatic nitrile monomers. The formation of the triazine unit through the trimerization of cyano groups served as the cross-linkers, resulting in the creation of π-conjugated networks alongside the arenes present in the starting materials. The aromatic nitrile monomers with rigid and sterically hindered structures were required to afford CTF membranes with nanoporous architectures. The acidity of the superacid and reactivity of the aromatic monomers played critical roles in the polymerization efficiency. The monomer diversity and synthesis tunability endowed the introduction of CO 2 -philic functionalities (e.g., pyrazole and fluorine) within the CTF skeletons, and integration of ionic moieties was achieved by adopting FSO 3 H with stronger acidity as the catalyst and aromatic nitrile monomers with pyrazine structures. To ensure the successful construction of fluorinated CTF membranes, it is important to avoid any fluorines on the ortho-position of the cyano groups on the benzene ring. Through control over the monomers and reaction conditions, flexible, transparent, and insoluble CTF membranes could be fabricated. The sol–gel method could be further expanded to membrane fabrication through acetyl-to-benzene transformation through synthesis control. The mild oxidation-exfoliation-filtration method was also demonstrated to fabricate substrate-supported CTF membranes. The as-afforded membranes are well characterized to determine the structural features and provide information to study the structure-performance relationship. Here, the application of CTF membranes in CO 2 separation was summarized, focusing on the approaches being developed to enhance CO 2 uptake and separation performance. In addition to utilizing the pristine CTF membranes for gas separation, functionalized carbon molecular sieve membranes could be obtained from the pyrolysis of thermally stable CTF membrane precursors toward efficient CO 2 separation, benefiting from the abundant ultramicropores being created during the pyrolysis/decomposition procedure and involvement of CO 2 -philic functionalities such as fluorine and nitrogen-containing moieties. Based on these achievements, unsolved issues in CTF membrane-related fabrication and applications, including the potential solution approaches, have been proposed to advance the application of CTF membranes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Localized Phonon Densities of States at Grain Boundaries in Silicon

Since it is now possible to record vibrational spectra at nanometer scales in the electron microscope, it is of interest to explore whether extended defects in crystals such as dislocations or grain boundaries will result in measurable changes of the phonon densities of states (dos) that are reflected in the spectra. Phonon densities of states were calculated for a set of high angle grain boundaries in silicon. The boundaries are modeled by supercells with up to 160 atoms, and the vibrational densities of states were calculated by taking the Fourier transform of the velocity–velocity autocorrelation function from molecular dynamics simulations with larger supercells doubled in all three directions. In selected cases, the results were checked on the original supercells by comparison with the densities of states obtained by diagonalizing the dynamical matrix calculated using density functional theory. Near the core of the grain boundary, the height of the optic phonon peak in the dos at 60 meV was suppressed relative to features due to acoustic phonons that are largely unchanged relative to their bulk values. This can be attributed to the variation in the strength of bonds in grain boundary core regions where there is a range of bond lengths.

Materials Science↗

DOE OSTI2023

Geospatial Information, Metadata, and Maps for Global River Corridor Science Focus Area Sites (v5)

This dataset provides geospatial information, metadata, and maps for the Pacific Northwest National Laboratory (PNNL) River Corridor Science Focus Area (RC-SFA; https://www.pnnl.gov/projects/river-corridor) sites. The RC-SFA works to transform understanding of spatial and temporal dynamics in river corridor hydrobiogeochemical functions from molecular reaction to watershed and basin scales. The knowledge we gain is used to formulate and test hypotheses and to improve mechanistic representation of river corridor processes and their response to disturbances in multiscale models of integrated hydrobiogeochemical function. The data provided includes Site ID, latitude, longitude, stream name, and common ID (COMID) for sites used across the RC-SFA. The COMID can be used to find and download data from NHDPlus (https://www.epa.gov/waterdata/nhdplus-national-hydrography-dataset-plus) and other platforms. The sites included are non-exhaustive. Sites (including past sites) will be added to this data package in the future. Data generated from the RC SFA can be accessed at https://data.ess-dive.lbl.gov/portals/PNNLRiverCorridorSFA. This data package was originally published in April 2023. It was updated in June 2023 (v2; modified files), December 2023 (v3; modified files), January 2025 (v4; modified files), and December 2025 (v5; modified files). See the change history section in the readme for more details. This dataset is comprised of one main data folder. The data folder consists of (1) file-level metadata; (2) data dictionary; (3) readme; (4) methods codes; (5) geospatial information for all RC SFA sites including International Generic Sample Number (IGSN); (6) maps of all sites and sites in Washington State, USA; and (7) a subfolder with the shapefile of all sites. All files are .csv, .pdf, .shp, .cpg, .dbf, .prj, .qmd, or .shx. We thank the Confederated Tribes and Bands of the Yakama Nation for access to field locations where some data were collected in Washington state. We also thank the Yakama Nation Tribal Council and Yakama Nation Fisheries for working with us to facilitate sample collection and optimization of data usage according to their values and worldview.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Real-time equation-of-motion CC cumulant and CC Green’s function simulations of photoemission spectra of water and water dimer

Newly developed coupled-cluster (CC) methods enable simulations of ionization potentials and spectral functions of molecular systems in a wide range of energy scales ranging from core-binding to valence. This paper discusses the results obtained with the real-time equation-of-motion CC cumulant (RT-EOM-CC) approach and CC Green’s function (CCGF) approaches in applications to the water and water dimer molecules. Here, we compare the ionization potentials obtained with these methods for the valence region with the results obtained with the coupled-cluster with singles, doubles, and perturbative triples formulation as a difference of energies for N and N – 1 electron systems. All methods show good agreement with each other. They also agree well with the experiment with errors usually below 0.1 eV for the ionization potentials. We also analyze unique features of the spectral functions, associated with the position of satellite peaks, obtained with the RT-EOM-CC and CCGF methods employing single and double excitations, as a function of the monomer OH bond length and the proton transfer coordinate in the dimer. Finally, we analyze the impact of the basis set effects on the quality of calculated ionization potentials and find that the basis set effects are less pronounced for the augmented-type sets.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Molecular-Weight-Dependent Infiltration and Adsorption of Polymers into Nanochannels

The mesoporous silica shell coating hydrogenolysis nano-catalysts alters the molecular weight distributions of cleaved polymer chains compared to catalysts without a shell. The shell, composed of radially aligned narrow cylindrical nanopores, reduces the formation of low-valued gaseous products and increases the median molecular weight of the product, thus enhancing the value of the products for polymer upcycling. Here, to understand the role of the mesoporous shell, we have studied the spatial distribution of polystyrene chains, used as a model polymer, in the nanochannels in both the melt phase and solution phase. In the melt, we observed from small-angle X-ray scattering experiments that the infiltration rate of the polymer into the nanochannels is inversely proportional to the molecular weight, which is consistent with theory. In theta solution experiments using UV–vis spectroscopy, we found that the shell significantly enhances polymer adsorption compared to nanoparticles without pores. In addition, the degree of polymer adsorption is not a monotonic function of molecular weight but initially increases with the molecular weight before eventually decreasing. The molecular weight for the peak adsorption increases with the pore diameter. This adsorption behavior is rationalized as resulting from a balance between the mixing entropy gain by surface adsorption and the conformational entropy penalty incurred by chains confined in the nanochannels. The spatial distribution of polymer chains in the nanochannels is visualized by energy-dispersive X-ray spectroscopy (EDX), and inverse Abel-transformed data reveals a less uniform polymer distribution along the primary pore axis for longer chains.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Intramolecular hydrogen bonding analysis

The quasi-atomic orbital (QUAO) bonding analysis is used to study intramolecular hydrogen bonding (IMHB) in salicylic acid and an intermediate that is crucial to the synthesis of aspirin. The bonding analysis rigorously explores IMHB through directly accessing information that is intrinsic to the molecular wave function, thereby bypassing the need for intrinsically biased methods. The variables that affect the strength of IMHB are determined using kinetic bond orders, QUAO populations, and QUAO hybridizations. Important properties include both the interatomic distance between hydrogen and oxygen participating in the IMHB and the hybridization on the oxygen. The bonding analysis further shows that each intramolecular hydrogen bond is a four-electron three-center bond. Here, the bonding analysis is used to understand how aromatic reactivity changes due to the effect of functional groups on the aromatic ring.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Evidence for Dissociation and Ionization in Shock Compressed Nitrogen to 800 GPa

Triple bonding in the nitrogen molecule (N 2 ) is among the strongest chemical bonds with a dissociation enthalpy of 9.8 eV/molecule. Nitrogen is therefore an excellent test bed for theoretical and numerical methods aimed at understanding how bonding evolves under the influence of the extreme pressures and temperatures of the warm dense matter regime. Here, we report laser-driven shock experiments on fluid molecular nitrogen up to 800 GPa and 4.0 g/cm 3 . Line-imaging velocimetry measurements and impedance matching method with a quartz reference yield shock equation of state data of initially precompressed nitrogen. Comparison with numerical simulations using path integral Monte Carlo and density functional theory molecular dynamics reveals clear signatures of chemical dissociation and the onset of L-shell ionization. Furthermore, combining data along multiple shock Hugoniot curves starting from densities between 0.76 and 1.29 g/cm 3 , our study documents how pressure and density affect these changes in chemical bonding and provides benchmarks for future theoretical developments in this regime, with applications for planetary interior modeling, high energy density science, and inertial confinement fusion research.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2020

Effects of CELF Pretreatment Severity on Lignin Structure and the Lignin-Based Polyurethane Properties

Conversion of technical lignin into performance biopolymers such as polyurethane offers environmental and economic advantages when combined with production of biofuels from biomass sugars, presenting significant interest toward studying the role of pretreatment on lignin structure and functionality. Co-solvent enhanced lignocellulosic fractionation (CELF) pretreatment, employing acidic aqueous tetrahydrofuran (THF) mixtures, was developed to effectively break down the lignin-carbohydrate matrix and promote extraction of lignin from lignocellulosic biomass with desirable purity and yield. In this study, we report the effects of CELF pretreatment reaction severity on the molecular structure of CELF-extracted lignin and its impact toward the mechanical properties of the resulting lignin-based polyurethanes. Reaction temperature was found to play the most significant role, compared to reaction time and acidity, in manipulating structural features such as molecular weight, functionality and intra-polymer structure. At the severe reaction conditions at 180°C, the order of reactivity for primary lignin interlinkages characterized by semiquantitative HSQC NMR analysis were found to be β-ether > phenylcoumaran (β−5′) > resinol (β−β′) facilitating a high degree of depolymerization and yielding a high frequency of free phenolics and reduced aliphatic hydroxyl groups. All side-chain interlinkages were depleted converting guaiacyl subunits into condensed forms, while still retaining uncondensed syringyl subunits. Under the mild 150°C temperature reaction, CELF lignin had higher molecular weight and retained more β-ether interlinkages. The results from CELF lignin-based polyurethane synthesis indicated that the tensile properties depended on the miscibility of CELF lignin with other components and low molecular weight cuts improved the dispersion of lignin in the polyurethane network. Pre-mixing of CELF with poly(ethylene glycol) (PEG) reduced the brittleness and improved the ductility of the CELF lignin-PEG polyurethanes.

59 BASIC BIOLOGICAL SCIENCES↗