Nearest-neighbor Interactions Research Articles

Three-dimensional checkerboard spin structure on a breathing pyrochlore lattice

The standard approach to realize a spin-liquid state is through magnetically frustrated states, relying on ingredients such as the lattice geometry, dimensionality, and magnetic interaction type of the spins. While Heisenberg spins on a pyrochlore lattice with only antiferromagnetic nearest-neighbor interactions are theoretically proven disordered, spins in real systems generally include longer-range interactions. The spatial correlations at longer distances typically stabilize a long-range order rather than enhancing a spin-liquid state. Both states can, however, be destroyed by short-range static correlations introduced by chemical disorder. Here, using disorder-free specimens with a clear long-range antiferromagnetic order, we refine the spin structure of the Heisenberg spinel ZnFe2O4 through neutron magnetic diffraction. The unique wave vector (1,0,12) leads to a spin structure that can be viewed as alternatively stacked ferromagnetic and antiferromagnetic tetrahedra in a three-dimensional checkerboard form. Stable coexistence of these opposing types of clusters is enabled by the bipartite breathing pyrochlore crystal structure, leading to a second-order phase transition at 10 K. The diffraction intensity of ZnFe2O4 is an exact complement to the inelastic scattering intensity of several chromate spinel systems which are regarded as model classical spin liquids. Our results challenge this attribution, and suggest instead of the six-spin ring mode, spin excitations in chromate spinels are closely related to the (1,0,12) type of spin order and the four-spin ferromagnetic cluster locally at one tetrahedron. Published by the American Physical Society 2024

The Effect of Polymer-Solvent Interaction on the Swelling of Polymer Matrix Tablets: A Magnetic Resonance Microscopy Study Complemented by Bond Fluctuation Model Simulations.

Polymer matrix tablets are an important drug-delivery system widely used for oral drug administration. Understanding the tablet hydration process, both experimentally and theoretically, is, thus, very important for the development of drug delivery systems that exhibit high drug loading capacity and controlled release potential. In this study, we used magnetic resonance microscopy (MRM) to nondestructively and dynamically analyze the water hydration process of xanthan-based tablets. The swelling process was characterized by well-resolved fronts of erosion, swelling, and penetration. The experimental results were complemented by numerical simulations of the polymer matrix hydration process. In the simulations, the polymer tablet matrix was modeled as an assembly of interacting chains with embedded drug particles, while its hydration process was mediated by interaction with solvent particles. The swelling dynamics were modeled within a Monte Carlo-based bond fluctuation model (BFM) that elegantly accounted for steric and nearest-neighbor interactions. This study provides an efficient experimental-theoretical approach for the study of polymer matrix swelling processes.

Persistence of weak ferromagnetism in antiferromagnetic systems on the body-centered octahedral lattice.

The possibility of existence of weak ferromagnetism in antiferromagnetic systems on the body-centered octahedral lattice is investigated in the framework of the corresponding spin-1/2J_{1}-J_{2} model in the recursive-lattice approach. The exact solution of the model is found and its phase diagram is determined. The magnetic and thermodynamic properties of all phases are studied, the nature of all phase transitions is established, and an equationthat determines the positions of all second-order phase transitions of the model is found. The magnetic and entropy properties of all ground states of the model are also determined. It is shown that the weak ferromagnetism predicted earlier in the pure spin-1/2 antiferromagnetic system on the octahedral lattice remains present even in the case of the model on the body-centered octahedral lattice with antiferromagnetic as well as ferromagnetic interactions between the central site of each elementary octahedron and each of its vertices. However, the presence of the interacting central site suppresses the weak ferromagnetism at very low temperatures, where the standard antiferromagnetic phase emerges. At the same time, the temperature region with the antiferromagnetic phase increases with simultaneous decreasing of the region, where the weak ferromagnetism can be observed, when the strength of the interaction between the central site of each elementary octahedron and each of its vertices increases. Moreover, the phenomenon of weak ferromagnetism disappears completely when this interaction becomes sufficiently stronger than the nearest-neighbor antiferromagnetic interaction between spin variables placed in vertices of each elementary octahedron of the lattice. Moreover, the possibility of the existence of the classical spin-liquid behavior in such magnetic systems is also discussed.

Digital quantum simulation of non-perturbative dynamics of open systems with orthogonal polynomials

Classical non-perturbative simulations of open quantum systems' dynamics face several scalability problems, namely, exponential scaling of the computational effort as a function of either the time length of the simulation or the size of the open system. In this work, we propose the use of the Time Evolving Density operator with Orthogonal Polynomials Algorithm (TEDOPA) on a quantum computer, which we term as Quantum TEDOPA (Q-TEDOPA), to simulate non-perturbative dynamics of open quantum systems linearly coupled to a bosonic environment (continuous phonon bath). By performing a change of basis of the Hamiltonian, the TEDOPA yields a chain of harmonic oscillators with only local nearest-neighbour interactions, making this algorithm suitable for implementation on quantum devices with limited qubit connectivity such as superconducting quantum processors. We analyse in detail the implementation of the TEDOPA on a quantum device and show that exponential scalings of computational resources can potentially be avoided for time-evolution simulations of the systems considered in this work. We applied the proposed method to the simulation of the exciton transport between two light-harvesting molecules in the regime of moderate coupling strength to a non-Markovian harmonic oscillator environment on an IBMQ device. Applications of the Q-TEDOPA span problems which can not be solved by perturbation techniques belonging to different areas, such as the dynamics of quantum biological systems and strongly correlated condensed matter systems.

Transition to period-3 synchronized state in coupled gauss maps.

We study coupled Gauss maps in one dimension with nearest-neighbor interactions. We observe transitions from spatiotemporal chaos to period-3 states in a coarse-grained sense and synchronized period-3 states. Synchronized fixed points are frequently observed in one dimension. However, synchronized periodic states are rare. The obvious reason is that it is very easy to create defects in one dimension. We characterize all transitions using the following order parameter. Let x∗ be the fixed point of the map. The values above (below) x∗ are classified as +1 (-1) spins. We expect all sites to return to the same band after three time steps for a coarse-grained periodic or three-period state. We define the flip rate F(t) as the fraction of sites i such that si(3t-3)≠si(t). It is zero in the coarse-grained periodic state. This state may or may not be synchronized. We observe three different transitions. (a) If different sites reach different bands, the transition is in the directed-percolation universality class. (b) If all sites reach the same band, we find an Ising-type transition. (c) A synchronized period-3 state where a new exponent is observed. We also study the finite-size scaling at critical points. The exponents obtained indicate that the synchronized period-3 transition is in a new universality class.

Modeling the effects of near-neighbor cooperative interactions on contractile dynamics in murine and porcine myocardium

Modeling the effects of near-neighbor cooperative interactions on contractile dynamics in murine and porcine myocardium

Inducing Ferrimagnetic Exchange in 1D-FeSe2 Chains Using Heteroleptic Amine Complexes: [Fe(en)(tren)][FeSe2]2.

[Fe(en)(tren)][FeSe2]2 (en = ethylenediamine, C2H8N2, tren = tris(2-aminoethyl)amine, C6H18N4) has been synthesized by a mixed-ligand solvothermal method. Its crystal structure contains heteroleptic [Fe(en)(tren)]2+ complexes with distorted octahedral coordination, incorporated between 1D-FeSe2 chains composed of edge-sharing FeSe4 tetrahedra. The twisted octahedral coordination environment of the Fe-amine complex leads to partial dimerization of Fe-Fe distances in the FeSe2 chains so that the FeSe4 polyhedra deviate strongly from the regular tetrahedral geometry. 57Fe Mössbauer spectroscopy reveals oxidation states of +3 for the Fechain atoms and +2 for the Fecomplex atoms. The close proximity of Fe atoms in the chains promotes ferromagnetic nearest neighbor interactions, as indicated by a positive Weiss constant, θ = +53.8(6) K, derived from the Curie-Weiss fitting. Magnetometry and heat capacity reveal two consecutive magnetic transitions below 10 K. DFT calculations suggest that the ordering observed at 4 K is due to antiferromagnetic intrachain interactions in the 1D-FeSe2 chains. The combination of two different ligands creates an asymmetric coordination environment that induces changes in the structure of the Fe-Se fragments. This synthetic strategy opens new ways to explore the effects of ligand field strength on the structure of both Fe-amine complexes and surrounding Fe-Se chains.

Investigation of Nearest Neighbor Interactions in Nickel Substituted Cobalt Using 59Co IFNMR

Co100-x Nix (x = 10, 20, 30, 40) alloys were synthesized via the chemical method by reducing nitrates of cobalt and nickel using hydrazine as a reducing agent. The proportion of solvents, ethanol, and water was maintained constant for all the samples in the series. The structure and morphology of the samples were examined based on X-Ray diffraction (XRD) data and Scanning Electron Microscope (SEM) images. From XRD it was evident that the system had a mixture of hexagonal close packed (hcp) and face centered cubic (fcc) phases. The amount of hcp phase in the composition gradually reduced as the concentration of nickel increased. The effect of substitution of nickel on cobalt was studied using a home built 59Co internal field nuclear magnetic resonance (IFNMR) spectrometer. Further, the saturation magnetization and hence magnetic field anisotropy values, obtained from M vs H loops, are found to vary with nickel concentration. The NMR results show a broadening of the NMR signals in comparison with the pure phase of cobalt. The NMR spectra and the number of peaks were assigned to possible cobalt-nickel neighbors, and the percentage of such combinations is interpreted in terms of the frequency and intensity of the peaks. This percentage of different nickel-cobalt combinations was compared with the random distribution obtained using the binomial distribution function. This comparison shows a clear trend compared to that analysed by the IFNMR experiment.

Electron doping as a handle to increase the Curie temperature in ferrimagnetic Mn3Si2X6 (X = Se, Te).

By analysing the results of ab initio simulations performed for Mn3Si2X6 (X = Se, Te), we first discuss the analogies and the differences in electronic and magnetic properties arising from the anion substitution, in terms of size, electronegativity, band widths of p electrons and spin-orbit coupling strengths. For example, through mean-field theory and simulations based on density functional theory, we demonstrate that magnetic frustration, known to be present in Mn3Si2Te6, also exists in Mn3Si2Se6 and leading to a ferrimagnetic ground state. Building on these results, we propose a strategy, electronic doping, to reduce the frustration and thus to increase the Curie temperature (TC). To this end, we first study the effect of electronic doping on the electronic structure and magnetic properties and discuss the differences in the two compounds, along with their causes. Secondly, we perform Monte-Carlo simulations, considering from the first to the fifth nearest-neighbor magnetic interactions and single-ion anisotropy, and show that electron doping efficiently raises the TC.

Two-dimensional ferromagnetism induced by Jahn-Teller distortion and orbital order

As the dimension of the system decreases, the quantum confinement effect and electronic correlation interaction inside the material will be correspondingly enhanced, often resulting in some novel physical properties. Recently, the freestanding perovskite oxide films down to the monolayer limit have been successfully prepared and can be transferred to any desired substrate, providing a great opportunity to explore the functional properties of two-dimensional perovskites. In perovskite materials, Jahn-Teller distortion and orbital order often cause a variety of correlated electronic behaviors. However, unlike van der Waals materials that retain their structural and chemical bonding characteristics when reduced to the monolayer limit, perovskite materials may undergo structural reconstruction when reduced to two dimensions. Therefore, whether Jahn-Teller distortion and related effects exist in the perovskite monolayer limit, and whether two-dimensional perovskite can exhibit some new properties different from its bulk phase, have become urgent issues to be solved. In this work, perovskite fluoride KCuF<sub>3</sub> and its monolayer have been comparatively studied by means of first-principles calculation, symmetry analysis, and Monte Carlo simulation methods, to reveal the change in lattice dynamics, structural, electronic, and magnetic properties caused by dimensionality reduction in perovskites. The results show that the cooperative Jahn-Teller distortion and the in-plane staggered orbital order occurring in the KCuF<sub>3</sub> bulk can be retained to the monolayer limit. However, unlike the bulk phase, the Jahn-Teller distortion mode appears as a soft mode of the prototype phase in the monolayer, and the insulating property of the monolayer does not rely on the emergence of the Jahn-Teller distortion, but is related to the enhancement of the electronic correlation effect. The staggered orbital order causes the nearest-neighbor exchange interaction to be ferromagnetic, resulting in the monolayer being a two-dimensional ferromagnetic insulator, different from the antiferromagnetic phase in the bulk. Monte Carlo simulations predict that the Curie temperature of the monolayer is about 5 K, which is much lower than the Néel temperature of the bulk phase, indicating that the disappearance of interlayer coupling leads to a significant reduction in the magnetic phase transition temperature. This work provides guidance and reference for the study of two-dimensional perovskite materials and the design of perovskite-based two-dimensional ferromagnets.

Theoretical study of long-range ordered vacancy distribution in two-dimensional boron structures

<sec>In a two-dimensional boron structure, the ordered high-concentration vacancy distribution can enhance structural stability and significantly modulates material properties. Based on recent experimental progress, herein we particularly focus on the two-dimensional boron structures with a striped distribution of hexagonal vacancies, in order to explore the formation of long-period boron structures.</sec><sec>Utilizing the structures of alloy generation and recognition (SAGAR) program developed by our group, we eliminate duplicate structures according to the structural symmetry to reduce computational cost. An effective model system is proposed to investigate the effect of vacancy distribution on the stability of the system, where the interactions between vacancies are utilized for estimating the total energy. By selecting structures with appropriate concentrations and combining first-principles calculations, the parameters in the model are fitted for different vacancy neighbor interactions, which can be further used to predict stable structures at various vacancy concentrations. The feasibility of model analysis is emphasized for structural screening, showing the good agreement between the parameterized model and the first-principles calculations.</sec><sec>Interestingly, under the same vacancy concentration, stable boron structures with different cell sizes exhibit distinct vacancy distributions, indicating a trend of long-period distribution for ground state structures. To address this phenomenon, when the stable candidate structures from the 1/6 series are dominant in number within the computable range and the changes in neighbor statistics can be clearly seen, we select the structures from this concentration series for detailed calculations.</sec><sec>The calculation results indicate that the convergence of the average energy is primarily influenced by the interaction between the fourth nearest neighbor and the sixth nearest neighbor. When considering only these two neighbors, the system energy changes with the increase of cell size as follows: the average energy of structures with a cell size being an even multiple of the minimum cell size keeps unchanged, while the average energy of structure with a cell size being an odd multiple of the minimum cell size gradually decreases, eventually converging to a stable value. When including the interactions between the ninth nearest neighbor and the tenth nearest neighbor, the average energy of structures with a cell size being an even times the minimum cell size also decreases gradually. The average energy decreases with oscillations, with the magnitude gradually diminishing and eventually stabilizing. This discovery reveals that the enhanced stability of long-period structures is attributed to the competitive interactions among different neighboring vacancies.</sec>

Pressure-induced structural and magnetic ordering transitions in the J1-J2 square lattice antiferromagnets AMoOPO4Cl (A = K, Rb).

By means of ab initio density functional theory calculations taking into account electronic correlation and van der Waals force, we conducted comprehensive studies of the electronic and magnetic properties, as well as structural and magnetic ordering evolution under pressure of the square lattice antiferromagnets AMoOPO4Cl (A = K, Rb) containing Mo5+ ions with , theoretically predicted as the potential candidates for achieving quantum phases, existing in the boundary regimes for square lattice magnets. Our results indicate that the columnar antiferromagnetic ordering, experimentally determined, is the magnetic ground state of the ambient P4/nmm phase, stabilized by the predominant antiferromagnetic next nearest neighbor interaction J2 in the diagonal directions of the square lattice, regardless of the effective Hubbard amendment values. More importantly, the P4/n phase, involving the mutual twisting of the MoO5Cl and PO4 polyhedra, satisfactorily reproduces the experimentally observed structural transition and the subsequent magnetic ordering transition from columnar antiferromagnetic ordering to Néel antiferromagnetic one, identified to be the appropriate high pressure structure. Furthermore, the mechanism underlined responsible for the magnetic ordering transition at high pressure has been disclosed in terms of density of states and spin density isosurface analysis across the transition. The loss of mirror plane symmetry in the P4/n phase activates the P 3s orbitals to participate in the magnetic interaction, giving rise to a competitive ferromagnetic superexchange interaction, in addition to antiferromagnetic direct one, and consequently initiating the magnetic ordering transition. The insights revealed here not only deepen our understanding of the electronic properties and structural and magnetic ordering transitions under high pressure of square lattice antiferromagnets AMoOPO4Cl (A = K, Rb), but also push the boundaries of knowledge by recognizing the role of nonmagnetic ions P 3s in magnetic exchange coupling.

Persistent current in a mesoscopic Holstein-Hubbard ring with Dresselhaus interaction

The effect of electron-phonon coupling, onsite repulsive Coulomb interaction and temperature on the persistent current in a quantum ring is studied in the presence of Dresselhaus spin-orbit interaction. The quantum ring threaded by the Aharonov-Bohm flux is modelled by the one-dimensional Holstein-Hubbard-Dresselhaus Hamiltonian. The electron-phonon interaction and Dresselhaus spin-orbit interaction are decoupled by employing the Lang-Firsov coherent transformation and a unitary transformation respectively. Thereafter, a self-consistent diagonalization technique is performed numerically at the Hartree-Fock level to obtain the effective electronic energy and current. It is shown that the intrinsic Dresselhaus spin-orbit interaction enhances the persistent charge and spin currents significantly. On the other hand, the persistent current is reduced by the onsite and nearest-neighbour electron-phonon interaction and Coulomb interaction. Also, the behaviour of the currents is modified by temperature. The spin-splitting of persistent spin current is enhanced considerably by Dresselhaus spin-orbit interaction and this splitting is tuneable in different regimes of magnetic flux, temperature, chemical potential and the interactions present in the system.

GROUND STATES FOR THE SOS MODEL WITH COMPETING BINARY INTERACTIONS ON A CAYLEY TREE OF ORDER THREE

We consider a SOS (solid-on-solid) model with nearest-neighbor interaction prolonged next-nearest-neighbor interaction and one level next-nearest-neighbor interaction where the spin takes values in the set on a Cayley tree of order three. In the paper, we study translation-invariant and periodic ground states of the model SOS.

Signature of elementary excitations in Se-substituted layered triangular-lattice antiferromagnet CuCrS2

Geometrical spin frustration in a layered triangular-lattice system prohibits the long-range magnetic order and could give rise to unusual spin-disordered states of matter in which spins are strongly correlated and highly entangled with low-energy excitations. Here we investigate the low-temperature magnetic and specific heat properties of the layered triangular-lattice system CuCr(S1-xSex)2 within the entire Se-for-S substitution range. In magnetization measurements, the x = 0 phase exhibits a sharp cusp-like antiferromagnetic transition at 38 K, while for all the Se-substituted samples a broad maximum along with a spin-glass-like freezing at low temperatures is seen. The Curie-Weiss temperature systematically increases from −122 K for x = 0 to + 15 K for x = 1. These observations suggest that the magnetic ground state is determined by the competition between antiferromagnetic direct exchange interactions and ferromagnetic super-exchange interactions. For x = 0, the magnetic contribution to heat capacity exhibits a sharp jump at the magnetic transition temperature indicative of 3D long-range magnetic order, while for the x > 0 samples it increases smoothly across the magnetic transition range following a power-law (∼T2.1) behavior. All these results are in line with an existence of unusual gapless spin-liquid like excitations originating from the spin frustrations and competing nearest-neighbor interactions in CuCr(S,Se)2.

Unveiling exotic magnetic phase diagram of a non-Heisenberg quasicrystal approximant

A magnetic phase diagram of the non-Heisenberg Tsai-type 1/1 Au-Ga-Tb approximant crystal (AC) has been established across a wide electron-per-atom (e/a) range via magnetization and powder neutron diffraction measurements. The diagram revealed exotic ferromagnetic (FM) and antiferromagnetic (AFM) orders that originate from the unique local spin icosahedron common to icosahedral quasicrystals (iQCs) and ACs; The noncoplanar whirling AFM order is stabilized as the ground state at the e/a of 1.72 or less whereas a noncoplanar whirling FM order was found at the larger e/a of 1.80, with magnetic moments tangential to the Tb icosahedron in both cases. Moreover, the FM/AFM phase selection rule was unveiled in terms of the nearest neighbour (J1) and next nearest neighbour (J2) interactions by numerical calculations on a non-Heisenberg single icosahedron. The present findings will pave the way for understanding the intriguing magnetic orders of not only non-Heisenberg FM/AFM ACs but also non-Heisenberg FM/AFM iQCs, the latter of which are yet to be discovered.

Ising Ladder with Four-Spin Plaquette Interaction in a Transverse Magnetic Field.

The spin-1/2 quantum transverse Ising model, defined on a ladder structure, with nearest-neighbor and four-spin interaction on a plaquette, was studied by using exact diagonalization on finite ladders together with finite-size-scaling procedures. The quantum phase transition between the ferromagnetic and paramagnetic phases has then been obtained by extrapolating the data to the thermodynamic limit. The critical transverse field decreases as the antiferromagnetic four-spin interaction increases and reaches a multicritical point. However, the exact diagonalization approach was not able to capture the essence of the dimer phase beyond the multicritical transition.

Spin-liquid behavior in the classical antiferromagnetic [formula omitted] spin system on the square-kagome lattice: Exact analysis within recursive-lattice approach

Magnetic and thermodynamic properties of the antiferromagnetic J1−J2−J3 spin system, i.e., the spin system with three different nearest-neighbor antiferromagnetic interactions, on the square-kagome lattice is investigated in the framework of the corresponding exactly solvable spin-1/2 Ising model on the recursive square-kagome lattice. The phase diagram of the model is found and it is shown that, depending on the model parameters, the model exhibits the existence of three long-range order (antiferromagnetic) phases, which are separated from the paramagnetic phase by the second order phase transitions. It is also shown that, in the zero temperature limit, the paramagnetic phase splits into three different ground states realized on the lines that separate antiferromagnetic ground states of the model. Due to the frustration, all these “paramagnetic” ground states are highly macroscopically degenerated with strict hierarchy of their residual entropies. Moreover, the high macroscopic degeneracy of the paramagnetic phase in the vicinity of these ground states together with the presence of strong correlations (given by the proximity of the second order phase transitions) allow one to consider the paramagnetic ground states as well as their immediate paramagnetic surroundings as regions with the classical spin-liquid behavior. It is also shown that the presence of the second order phase transitions in the vicinity of the paramagnetic ground states suppresses expected anomalous (Schottky-type) behavior of the specific heat related to the rapid entropy changes caused by the frustration.

Effects of phosphorylation on Drp1 activation by its receptors, actin, and cardiolipin.

Drp1 is a dynamin family GTPase required for mitochondrial and peroxisomal division. Oligomerization increases Drp1 GTPase activity through interactions between neighboring GTPase domains. In cells, Drp1 is regulated by several factors including Drp1 receptors, actin filaments, cardiolipin, and phosphorylation at two sites: S579 and S600. Commonly, phosphorylation of S579 is considered activating, while S600 phosphorylation is considered inhibiting. However, direct effects of phosphorylation on Drp1 GTPase activity have not been investigated in detail. Here, we compare effects of S579 and S600 phosphorylation on purified Drp1, using phosphomimetic mutants and in vitro phosphorylation. Both phosphomimetic mutants are shifted toward smaller oligomers. Both phosphomimetic mutations maintain basal GTPase activity, but eliminate GTPase stimulation by actin and decrease GTPase stimulation by cardiolipin, Mff, and MiD49. Phosphorylation of S579 by Erk2 produces similar effects. When mixed with wildtype Drp1, both S579D and S600D phosphomimetic mutants reduce the actin-stimulated GTPase activity of Drp1-WT. Conversely, a Drp1 mutant (K38A) lacking GTPase activity stimulates Drp1-WT GTPase activity under both basal and actin-stimulated conditions. These results suggest that the effect of S579 phosphorylation is not to activate Drp1 directly. In addition, our results suggest that nearest neighbor interactions within the Drp1 oligomer affect catalytic activity.

Pair density wave and superconductivity in a kinetically frustrated doped Emery model on a square lattice

The quest to understand the nature of superconductivity in the cuprates has spotlighted the pair density wave (PDW)–a superconducting state characterized by a spatially modulated order parameter. Despite significant advances in understanding PDW properties, conclusively demonstrating its presence in systems pertinent to cuprate superconductors remains elusive. In this study, we present a systematic density-matrix renormalization group study to investigate the Emery model (or the three-band Hubbard model) on two-leg square cylinders with negative electron hopping term tpp between adjacent oxygen sites. Kinetic frustration - introduced by changing the sign of oxygen-oxygen hopping - leads to a much reduced Cu-Cu antiferromagnetic exchange along with an enlarged charge transfer energy that changes the local properties of the model. At light doping levels, our findings reveal a ground state remarkably consistent with a PDW, exhibiting mutually commensurate superconducting (SC), charge, and spin density wave correlations. Intriguingly, the dominant SC pairing is observed between neighboring oxygen sites, diverging from the expected Cu sites in the positive tpp case. When the system incorporates moderate near-neighbor interactions, particularly an attractive Vpp between adjacent oxygen sites, the SC correlations become quasi-long-ranged, accompanied by a pronounced divergence in the PDW susceptibility. When the attractive Vpp increases further, the system gives way to an unconventional d-wave superconductivity.