Impact of Atomic Defects on Ceria Surfaces on Chemical Mechanical Polishing of Silica Glass Surfaces.

This study investigates the impact of atomic defects, such as oxygen vacancies and Ce3+ ions, on cerium oxide (ceria) surfaces during chemical mechanical polishing (CMP) for silica glass finishing. Using density functional theory (DFT) and reactive molecular dynamics simulations, the interaction of orthosilicic molecules and silica glass with dry and wet ceria surfaces is explored. Defects alter the surface reactivity, leading to the dissociation of orthosilicic acid on oxygen vacancies, forming a strong Si-O-Ce bond. Hydroxylated surfaces exhibit easier oxygen vacancy formation and thermodynamically favored substitution of hydroxyl groups with orthosilicic acid. A new ReaxFF library for silica/ceria interfaces with defects is validated using DFT outcomes. Reactive MD simulations demonstrate that ceria surfaces with 30% Ce3+ ions on (111) planes exhibit higher polishing efficiency, attributed to increased Si-O-Ce bond formation. The simultaneous presence of oxygen vacancies and various acidic and basic sites on ceria surfaces enhances the polishing efficiency, involving acid-base reactions with silica. Defective surfaces show superior efficiency by removing silicate chains, contrasting with nondefective surfaces removing isolated orthosilicate units. This study provides insights into optimizing CMP processes for high-precision glass industry surface finishing.

Atomistic mechanism of structural and volume relaxation below glass transition temperature in a soda-lime silicate glass revealed by Raman spectroscopy and its DFT calculations.

To elucidate the atomistic origin of volume relaxation in soda-lime silicate glass annealed below the glass transition temperature (Tg), the experimental and calculated Raman spectra were compared. By decomposing the calculated Raman spectra into specific groups of atoms, the Raman peaks at 800, 950, 1050, 1100, and 1150 cm-1 were attributed to oxygen and silicon in Si-O-Si, non-bridging oxygen in the Q2 unit, bridging oxygen in low-angle Si-O-Si, non-bridging oxygen in the Q4 unit, and bridging oxygen in high-angle Si-O-Si, respectively. Based on these attributions, we found that by decreasing the fictive temperature by annealing below Tg - 70 K, a homogenization reaction Q2 + Q4 → 2Q3 and an increase in average Si-O-Si angle occurred simultaneously. By molecular dynamics simulation, we clarified how the experimentally demonstrated increase in average Si-O-Si angle contributes to volume shrinkage; increasing Si-O-Si angles can expand the space inside the rings, and Na can be inserted into the ring center.

Machine learning molecular dynamics reveals the structural origin of the first sharp diffraction peak in high-density silica glasses.

The first sharp diffraction peak (FSDP) in the total structure factor has long been regarded as a characteristic feature of medium-range order (MRO) in amorphous materials with a polyhedron network, and its underlying structural origin is a subject of ongoing debate. In this study, we utilized machine learning molecular dynamics (MLMD) simulations to explore the origin of FSDP in two typical high-density silica glasses: silica glass under pressure and permanently densified glass. Our MLMD simulations accurately reproduce the structural properties of high-density silica glasses observed in experiments, including changes in the FSDP intensity depending on the compression temperature. By analyzing the simulated silica glass structures, we uncover the structural origin responsible for the changes in the MRO at high density in terms of the periodicity between the ring centers and the shape of the rings. The reduction or enhancement of MRO in the high-density silica glasses can be attributed to how the rings deform under compression.

Three-dimensional ordering of water molecules reflecting hydroxyl groups on sapphire (001) and α-quartz (100) surfaces.

Water molecules on oxide surfaces influence the chemical reactivity and molecular adsorption behavior of oxides. Herein, three-dimensional atomic force microscopy (3D-AFM) and molecular dynamics simulations are used to visualize the surface hydroxyl (OH) groups and their hydration structures on sapphire (001) and α-quartz (100) surfaces at the atomic-scale. The obtained results revealed that the spatial density distributions and hydrogen-bonding strengths of surface OH groups affect their local hydration structures. In particular, the force curves obtained by 3D-AFM suggest that the hydration forces of water molecules intensify at sites where water molecules strongly interact with the surface OH groups. The insights obtained in this study deepen our understanding of the affinities of Al2O3 and SiO2 for water molecules and contribute to the use of 3D-AFM in the investigation of atomic-scale hydration structures on various surfaces, thereby benefiting a wide range of research fields dealing with solid-liquid interfaces.

New Atomistic Insights on the Chemical Mechanical Polishing of Silica Glass with Ceria Nanoparticles.

Reactive molecular dynamics simulations have been used to simulate the chemical mechanical polishing (CMP) process of silica glass surfaces with the ceria (111) and (100) surfaces, which are predominantly found in ceria nanoparticles. Since it is known that an alteration layer is formed at the glass surface as a consequence of the chemical interactions with the slurry solutions used for polishing, we have created several glass surface models with different degrees of hydroxylation and porosity for investigating their morphology and chemistry after the interaction with acidic, neutral, and basic water solutions and the ceria surfaces. Both the chemical and mechanical effects under different pressure and temperature conditions have been studied and clarified. According to the simulation results, we have found that the silica slab with a higher degree of hydroxylation (thicker alteration layer) is more reactive, suggesting that proper chemical treatment is fundamental to augment the polishing efficiency. The reactivity between the silica and ceria (111) surfaces is higher at neutral pH since more OH groups present at the two surfaces increased the Si-O-Ce bonds formed at the interface. Usually, an outermost tetrahedral silicate unit connected to the rest of the silicate network through a single bond was removed during the polishing simulations. We observed that higher pressure and temperature accelerated the removal of more SiO4 units. However, excessively high pressure was found to be detrimental since the heterogeneous detachment of SiO4 units led to rougher surfaces and breakage of the Si-O-Si bond, even in the bulk of the glass. Despite the lower concentration of Ce ions at the surface resulting in the lower amount of Si-O-Ce formed, the (100) ceria surface was intrinsically more reactive than (111). The different atomic-scale mechanisms of silica removal at the two ceria surfaces were described and discussed.

Evidence of Multiple Crystallization Pathways in Lithium Disilicate: A Metadynamics Investigation.

Metadynamics simulations driven by using two X-ray diffraction peaks identified three alternative crystallization pathways of the lithium disilicate crystal from the melt. The most favorable one passes through the formation of disordered layered structures undergoing internal ordering in a second step. The second pathway involves the formation of phase-separated structures composed of nuclei of ß-cristobalite crystals surrounded by lithium-rich phases in which metasilicate chains are formed. The conversion of these structures to the stable lithium disilicate crystal involves an intermediate structure whose silicate layers are connected by silicate rings with the energy barrier of 2.5 kJ/mol per formula unit (f.u.). The third pathway is highly unlikely because of the huge energy barrier involved (20 kJ/mol per f.u.). This path also involves the passage through a phase-separated structure of an indefinite silica region surrounded mainly by amorphous lithium oxide.

Composition Effect on Interfacial Reactions of Sodium Aluminosilicate Glasses in Aqueous Solution.

Understanding the underlying reaction mechanisms responsible for aluminosilicate glass dissolution in aqueous environments is crucial for designing glasses for technological applications ranging from architecture windows and touch screens to nuclear waste disposal. This study investigated the glass composition effect on the interfacial reactions of sodium aluminosilicate (NAS) glasses using molecular dynamics (MD) simulations with recently developed reactive potentials. Glass-water interfacial models of six NAS glasses with varying Al2O3/Na2O ratios were investigated for up to 4 nanoseconds (ns) to elucidate the interfacial reaction mechanisms at ambient temperature. The results showed that the coordination defects, such as undercoordinated Si and Al, as well as non-bridging oxygens (NBOs) accumulated at the glass surfaces, play a crucial role in the initial hydration reaction process of the glasses. They promote the formation of silanol (Si-OH) and aluminol (Al-OH) species together with the Na+⇔ H+ ion-exchange reactions. The z-density profiles of H2O and H+ ions affirmed the water/H+ propagation into the glass up to 2 nanometers after 4 ns reactions. The penetration depth depends on the composition and shows a nonlinear dependence, suggesting that the subsequent water penetration, particularly into the bulk glass, is supported by the availability of random channels. Aluminol formations, including Al-OH or Al-OH2 near the surface, were found to form mainly through the hydrolysis of Al-O-Al bonds and hydration of Al+-NBO- units. While water molecules are involved in initial interfacial reactions, water penetration into the bulk glass region is primarily achieved by proton transfer. Compared to highly mobile proton transfer involving silanol groups, proton transfer associated with [AlO4]- species is much more limited, particularly in the bulk glass region. These new insights into the role of aluminum in interfacial reactions of the NAS glasses can help to understand the initial dissolution mechanisms and in designing more durable glasses.

Composition Dependence of the Atomic Structures and Properties of Sodium Aluminosilicate Glasses: Molecular Dynamics Simulations with Reactive and Nonreactive Potentials.

Understanding the composition-structure-property relations of glass materials is essential for their technological applications. In this study, the structures and properties of a series of sodium aluminosilicate glasses with varying Al2O3/Na2O ratios ((35 - x)Na2O-xAl2O3-65SiO2, x = 0, 5, 10, 15, 17.5, 20) covering peralkaline to peraluminous compositions, have been studied by using molecular dynamics simulations with two types of interatomic potentials: a fixed partial charge pairwise potential (Teter) and a reactive diffusive charge reactive potential (DCRP). The short and medium structural features such as bond lengths, coordination numbers, Qn distributions, and ring size distributions were obtained and compared with experimental data. It was found that silicon remained fourfold-coordinated throughout the compositional range, while a noticeable amount of fivefold-coordinated aluminum together with oxygen triclusters (TBO) are present in compositions with higher Al2O3 contents (RAl/Na > 1). In addition, the simulation results from both potentials show a certain level of violation of the Al avoidance rule by exhibiting a non-negligible amount of [AlOx]-[AlOx] polyhedral connections. Neutron and X-ray diffraction structure factors of the simulated glasses were calculated and compared with available experimental data. The mechanical properties, including Bulk, Shear, and Young's modulus, were calculated and found to increase with increasing RAl/Na, in good agreement with the experiments. Correlations of the properties with glass structures as a function of glass compositions and the advantages as well as potential issues of the two sets of potentials in modeling sodium aluminosilicate glasses are discussed in the context of features of glass structures and the prospect of future simulations of glass-water reactions.

Prognostic value of pre-transplantation total metabolic tumor volume on 18fluoro-2-deoxy-D-glucose positron emission tomography-computed tomography in relapsed and refractory aggressive lymphoma.

Relapsed and refractory aggressive lymphoma have a poor prognosis. High-dose chemotherapy followed by autologous hematopoietic stem cell transplantation (auto-HSCT) is effective in chemosensitive patients. Allogeneic hematopoietic stem cell transplantation (allo-HSCT) is among the few options for non-chemosensitive patients. 18Fluoro-2-deoxy-D-glucose positron emission tomography-computed tomography (18FDG-PET/CT) is the standard tool for evaluating response to chemotherapy and residual tumor volume. However, accurate assessment of residual tumor volume is not currently being achieved in clinical practice, and its value in prognostic and therapeutic stratification remains unclear. To answer this question, we investigated the efficacy of quantitative indicators, including total metabolic tumor volume (TMTV), in predicting prognosis after auto-HSCT and allo-HSCT. We retrospectively analyzed 39 patients who received auto-HSCT and 28 who received allo-HSCT. In the auto-HSCT group, patients with a higher TMTV had a poor prognosis due to greater risk of relapse. In the allo-HSCT group, patients with a higher TMTV had a lower progression-free survival rate and a significantly higher relapse rate. Neither Deauville score nor other clinical parameters were associated with prognosis in either group. Therefore, pre-transplant TMTV on PET is effective for prognostic prediction and therapeutic decision-making for relapsed or refractory aggressive lymphoma.

Adsorption characteristics of peptides on ω-functionalized self-assembled monolayers: a molecular dynamics study.

Molecular dynamics simulations were employed to investigate the adsorption behavior of a variety of amino-acid side-chain analogs (SCAs) and a ß-hairpin (HP7) peptide on a series of liquid-like self-assembled monolayers (SAMs) with terminal functional groups of -OH, -OCH3, -CH3, and -CF3. The relationships between the adsorption free energy of the SCAs and the interfacial properties of water on the SAMs were examined to determine the acute predictors of protein adsorption on the SAM surfaces. The structural changes of HP7 on the SAM surfaces were also investigated to understand the relationship between the surface nature and protein denaturation. It was found that the adsorption free energy of the SCAs was linearly related to the surface hydrophobicity, which was computed as the free energy of cavity formation near the SAM-water interfaces. In addition, the hydrophobic -CH3 and -CF3 SAMs produced substantial conformational changes in HP7 because of the strong hydrophobic attractions to the nonpolar side chains. The hydrophilic surface terminated by -OH also promoted structural changes in HP7 resulting from the formation of hydrogen bonds between the hydrophilic tail and HP7. Consequently, the moderate amphiphilic surface terminated by -OCH3 avoided the denaturation of HP7 most efficiently, thus improving the biocompatibility of the surface. In conclusion, these results provide a deep understanding of protein adsorption for a wide range of polymeric surfaces, and they can potentially aid the design of appropriate biocompatible coatings for medical applications.

Biasing crystallization in fused silica: An assessment of optimal metadynamics parameters.

Metadynamics (MetaD) is a useful technique to study rare events such as crystallization. It has been only recently applied to study nucleation and crystallization in glass-forming liquids such as silicates, but the optimal set of parameters to drive crystallization and obtain converged free energy surfaces is still unexplored. In this work, we systematically investigated the effects of the simulation conditions to efficiently study the thermodynamics and mechanism of crystallization in highly viscous systems. As a prototype system, we used fused silica, which easily crystallizes to ß-cristobalite through MetaD simulations, owing to its simple microstructure. We investigated the influence of the height, width, and bias factor used to define the biasing Gaussian potential, as well as the effects of the temperature and system size on the results. Among these parameters, the bias factor and temperature seem to be most effective in sampling the free energy landscape of melt to crystal transition and reaching convergence more quickly. We also demonstrate that the temperature rescaling from T > Tm is a reliable approach to recover free energy surfaces below Tm, provided that the temperature gap is below 600 K and the configurational space has been properly sampled. Finally, albeit a complete crystallization is hard to achieve with large simulation boxes, these can be reliably and effectively exploited to study the first stages of nucleation.

H2O2adsorption and dissociation on various CeO2(111) surface models: a first-principles study.

Periodic density functional theory (DFT) calculations using the hybrid PBE0 functional and atom-centered Gaussian functions as basis sets were carried out to investigate the absorption and the first steps involved in the decomposition of hydrogen peroxide (H2O2) on three different models of the ceria (111) surface. One of the models is a clean surface, and the others are defective and partially hydroxylated ceria surfaces. On the clean surface, we found that the minimum energy path of hydrogen peroxide decomposition involves a three-step process, i.e., adsorption, deprotonation, and formation of the peroxide anion, stabilized through its interaction with the surface at a Ce (IV) site, with activation barriers of less than about 0.5 eV. The subsequent formation of superoxide anions and molecular oxygen species is attributed to electron transfer from the reactants to the Ce (IV) ions underneath. On the defective surface, H2O2dissociation is an energetically downhill reaction thermodynamically driven by the healing of the O vacancies, after the reduction and decomposition of H2O2into oxygen and water. On the hydroxylated surface, H2O2is first adsorbed by forming a favorable H-bond and then undergoes heterolytic dissociation, forming two hydroxyl groups at two vicinal Ce sites.

A new universal force-field for the Li2S-P2S5 system.

Lithium thiophosphate electrolyte is a promising material for application in all-solid-state batteries. Ab initio molecular dynamics (AIMD) simulations have been used to investigate the ion conduction mechanisms in single-crystalline and glassy compounds. However, the complexity of real materials (e.g., materials with grain boundaries and multiphase glass-ceramics) causes AIMD simulations to have high computational cost. To overcome this computational limitation, we developed a new interatomic potential for classical molecular dynamics (CMD) simulations of Li solid-state electrolytes. The training datasets were generated from representative sulfide electrolytes (ß-Li3PS4, γ-Li3PS4, Li4P2S6, Li7P3S11, and Li7PS6 crystals and 70Li2S-30P2S5 glass). Using the functional forms of the Class II and Stillinger-Weber potentials, all parameters were optimized by minimizing the differences in forces on atoms, stresses, and potential energies between the CMD and AIMD results. Subsequent validation showed that the optimized parameters can reproduce the dynamics of Li+ as well as the structures of the crystalline and glassy materials. The ionic conductivity of Li7P3S11 crystal was approximately five times that of the isostoichiometric 70Li2S-30P2S5 glass, indicating that CMD simulations using the developed force-field accurately reproduced the effective conduction path in Li7P3S11 from AIMD. The developed force-field parameters make it possible to simulate complex materials including amorphous-crystalline interfaces and multiphase glass-ceramics in the CMD framework.

Mathematical model and numerical analysis method for dynamic fracture in a residual stress field.

Residual stress field is a self-equilibrium state of stress in the bulk solid material with the inhomogeneous field of the inelastic deformations. The high level of tensile residual stress often leads to dynamic fracture resulting in the instantaneous and catastrophic destruction of the materials because the cracks are fed with the strain energy initially stored in the bulk materials due to the residual stress. The dissipation of the strain energy with crack growth results in the release and the redistribution of the residual stress. In this paper, we propose an effective mathematical model and a numerical analysis method for dynamic fracture in residual stress field. We formulate the dynamic behavior of solid continuum with residual stress field in the context of particle discretization scheme finite element method. This formulation enables the appropriate evaluation of (i) release and redistribution of residual stress due to dynamic propagation of the cracks and (ii) the effect of the elastic wave on crack propagation, which are the most substantial problems on dynamic fracture in residual stress field. We perform the experiments and the simulations of dynamic fracture process in chemically tempered glass sheets with residual stress field to validate the proposed numerical analysis method. The simulation results show remarkable agreement with the experiments of the catastrophic failure of the glass sheets with residual stress field in all aspects of crack behavior. These results indicate that the proposed model and method can rigorously evaluate the release and the autonomous redistribution of the residual stress in the dynamic fracture process.

Structural origin of thermal shrinkage in soda-lime silicate glass below the glass transition temperature: A theoretical investigation by microsecond timescale molecular dynamics simulations.

Microscopic dynamical features in the relaxation of glass structures are one of the most important unsolved problems in condensed matter physics. Although the structural relaxation processes in the vicinity of glass transition temperature are phenomenologically expressed by the Kohlrausch-Williams-Watts function and the relaxation time can be successfully interpreted by Adam-Gibbs theory and/or Narayanaswamy's model, the atomic rearrangement, which is the origin of the volume change, and its driving force have not been elucidated. Using the microsecond time-scale molecular dynamics simulations, this study provides insights to quantitatively determine the origin of the thermal shrinkage below Tg in a soda-lime silicate glass. We found that during annealing below Tg, Na ions penetrate into the six-membered silicate rings, which remedies the acute O-O-O angles of the energetically unstable rings. The ring structure change makes the space to possess the cation inside the rings, but the ring volume is eventually reduced, which results in thermal shrinkage of the soda-lime silica glass. In conclusion, the dynamical structural relaxation due to the cation displacement evokes the overall volume relaxation at low temperature in the glassy material.

Simulation of Catastrophic Failure in a Residual Stress Field.

Residual stress has been empirically utilized for industrial applications to control material strength and shape of fragments. The interaction between the dynamically growing cracks and the residual stress field is sufficiently complicated to prevent us from building effective models. To rigorously evaluate the release and redistribution of residual stress in the dynamic fracture process, we develop a mathematical model and a numerical analysis method for the dynamic fracture in a residual stress field. Our methodology is simple and rigorous and applicable regardless of materials and scales.

Development and Application of a ReaxFF Reactive Force Field for Cerium Oxide/Water Interfaces.

Ceria (CeO2) is a well-known catalytic oxide with many environmental, energy production, and industrial applications, most of them involving water as a reactant, byproduct, solvent, or simple spectator. In this work, we parameterized a Ce/O/H ReaxFF for the study of ceria and ceria/water interfaces. The parameters were fitted to an ab initio training set obtained at the DFT/PBE0 level, including the structures, cohesive energies, and elastic properties of the crystalline phases Ce, CeO2, and Ce2O3; the O-defective structures and energies of vacancy formation on CeO2 bulk and CeO2 (111) surface, as well as the absorption and reaction energies of H2 and H2O molecules on CeO2 (111). The new potential reproduced reasonably well all the fitted properties as well as the relative stabilities of the different ceria surfaces, the oxygen vacancies formation, and the energies and structures of associative and dissociative water molecules on them. Molecular dynamics simulations of the liquid water on the CeO2 (111) and CeO2 (100) surfaces were carried out to study the coverage and the mechanism of water dissociation. After equilibration, on average, 35% of surface sites of CeO2 (111) are hydroxylated whereas 15% of them are saturated with molecular water associatively adsorbed. As for the CeO2 (100) surface, we observed that water preferentially dissociates covering 90% of the available surface sites in excellent agreement with recent experimental findings.

Self-assembly of the cationic surfactant n-hexadecyl-trimethylammonium chloride in methyltrimethoxysilane aqueous solution: classical and reactive molecular dynamics simulations.

A flexible aerogel polymerized from methyltrimethoxysilane (MTMS) shows great promise as a high-performance insulator owing to its substantially low thermal conductivity and mechanical flexibility, attributed to its porous microstructure and organic-inorganic hybridization, respectively, which promote its industrial applications. Conventionally, the cationic surfactant n-hexadecyltrimethylammonium chloride (CTAC) is utilized to experimentally control the nanoscale microstructure and, consequently, the flexibility of the MTMS aerogel; however, the mechanism through which CTAC prevents MTMS aggregation in the solution is not yet fully understood. This study unravels the role of CTAC in preventing MTMS aggregation in aqueous solution using both classical and reactive molecular dynamics simulations. We found that CTAC molecules can form self-aggregates even when the polymerization of MTMS progresses and then the MTMS-derived oligomer turns to be hydrophobic in aqueous solution. In summary, the self-assemblies of CTAC disperse among the MTMS associations and effectively prevent MTMS clustering, and this is considered as the key mechanism underlying the formation of a flexible microstructure of the hybrid aerogel.

Effects of Side-Chain Spacing and Length on Hydration States of Poly(2-methoxyethyl acrylate) Analogues: A Molecular Dynamics Study.

Hydration states of polymers are known to directly influence the adsorption of biomolecules. Particularly, intermediate water (IW) has been found able to prevent protein adsorption. Experimental studies have examined the IW content and nonthrombogenicity of poly(2-methoxyethyl acrylate) analogues with different side-chain spacings and lengths, which are HPx (x is the number of backbone carbons in a monomer) and PMCyA (y is the number of carbons in-between ester and ether oxygens of the side-chain) series, respectively. HPx was reported to possess more IW content but lower nonthrombogenicity compared to PMCyA with analogous composition. To understand the reason for the conflict, molecular dynamics simulations were conducted to elucidate the difference in the properties between the HPx and PMCyA. Simulation results showed that the presence of more methylene groups in the side chain more effectively prohibits water penetration in the polymer than those in the polymer backbone, causing a lower IW content in the PMCyA. At a high water content, the methoxy oxygen in the shorter side chain of the HPx cannot effectively bind water compared to that in the PMCyA side chain. HPx side chains may have more room to contact with molecules other than water (e.g., proteins), causing experimentally less nonthrombogenicity of HPx than that of PMCyA. In summary, theoretical simulations successfully explained the difference in the effects of side-chain spacing and length in atomistic scale.

Outcomes of methotrexate-associated lymphoproliferative disorders in rheumatoid arthritis patients treated with disease-modifying anti-rheumatic drugs.

Recently, the use of targeted synthetic or biological disease-modifying anti-rheumatic drugs (ts/bDMARDs) in addition to conventional synthetic (cs)DMARDs including methotrexate (MTX) for rheumatoid arthritis (RA) has increased. However, whether ts/bDMARDs are associated with the development and clinicopathological features of MTX-associated lymphoproliferative disorder (MTX-LPD) in patients with RA remains unknown. Therefore, we evaluated the clinical outcomes of 121 patients with MTX-LPD. Results showed that prior use of ts/bDMARDs was not associated with the different histopathological subtypes of MTX-LPD. Patients with polymorphic-type LPD had a better event-free survival than those with diffuse large B-cell lymphoma (DLBCL), classical Hodgkin lymphoma and peripheral T-cell lymphoma. The pathological subtype of lymphoma could predict the clinical outcome of MTX-LPD. In patients with DLBCL, the use of tumour necrosis factor-alpha (TNF-α) inhibitors prior to MTX-LPD onset was associated with a higher non-relapse mortality. Further, patients with RA previously treated with Janus kinase (JAK) inhibitors more commonly required chemotherapy than those treated with csDMARDs alone, indicating disease aggressiveness. Hence, special caution should be observed when managing patients with MTX-LPD previously treated with JAK or TNF-α inhibitors for RA.