Modulator photoelectric combined simulation method and device, computer equipment and storage medium

Abstract

The application relates to a modulator photoelectric combined simulation method and device, computer equipment and a storage medium, the method comprising the following steps: acquiring an input radio frequency signal of the modulator; inputting the radio frequency signal into a radio frequency time domain coupling mode equation to acquire a forward propagation radio frequency wave amplitude and a backward propagation radio frequency wave amplitude in the modulator; acquiring a radio frequency voltage distribution in the modulator according to the forward propagation radio frequency wave amplitude and the backward propagation radio frequency wave amplitude; acquiring a first light wave amplitude and a second light wave amplitude in the modulator according to the radio frequency voltage distribution and a light wave time domain coupling mode equation; acquiring an output optical signal of the modulator according to the first light wave amplitude and the second light wave amplitude; and obtaining a simulated eye diagram according to the optical signal. The application can simultaneously realize electrical simulation and optical simulation in the modulator.

Description

Technical Field

[0001] This invention relates to the field of simulation, and in particular to a modulator optoelectronic co-simulation method, apparatus, computer equipment, and storage medium. Background Technology

[0002] In the fields of optical communication and optical information processing, modulators are key components for encoding signals onto optical waves. Modulators can be manufactured using a variety of materials, and simulation is crucial for their design and manufacturing. It allows for the preliminary evaluation of modulator performance based on eye diagrams before manufacturing, providing guidance for design improvements.

[0003] Modulators themselves involve technical fields such as optics and circuits, so the simulation of modulators requires both electrical and optical simulations, hence the term optoelectronic joint simulation.

[0004] In the existing technology, the simulation of modulators mainly relies on two methods: the segmented circuit method and the transmission line circuit method, but these methods have significant drawbacks.

[0005] In the segmented circuit method, optical devices are transformed into equivalent circuit devices, and the optical waveguide and electrodes are segmented into an array of equivalent circuit devices to model the propagation of radio frequency waves and optical waves. Simulation is then performed in software such as Verilog-A or SPICE to achieve optoelectronic co-simulation. However, due to the significant physical differences between current-voltage (IV) signals and wave signals, this method struggles to handle the characteristics of wave signals, such as multi-dimensional and bidirectional propagation. In the transmission line circuit method, electrodes are treated using an equivalent circuit model of the transmission line. For electrode dispersion, additional equivalent circuit devices need to be introduced into the transmission line, and convolution calculations are performed in each computational step. Each convolution requires accessing all stored data, making it extremely time-consuming. For silicon photonic modulators, the nonlinear IV equations of the PN junction are also introduced, further increasing computational overhead and complexity. Therefore, this method is inefficient in handling the nonlinear IV relationships and dispersion of electrodes.

[0006] In addition, both of the above methods are material-dependent because they require the creation of specific circuit models for each material, thus lacking versatility.

[0007] It is evident that existing optoelectronic simulation methods are inefficient and highly dependent on materials. Summary of the Invention

[0008] To solve the above-mentioned technical problems, or at least partially solve them, the present invention provides a modulator optoelectronic co-simulation device, computer equipment, and storage medium.

[0009] In a first aspect, the present invention provides a modulator optoelectronic co-simulation method, wherein the modulator receives a radio frequency signal and outputs an optical signal, the method comprising:

[0010] Acquire the radio frequency signal input to the modulator;

[0011] The radio frequency signal is input into the radio frequency time-domain coupled mode equation to obtain the forward propagation radio frequency wave amplitude and the reverse propagation radio frequency wave amplitude in the modulator;

[0012] The radio frequency voltage distribution within the modulator is obtained based on the forward propagation radio frequency wave amplitude and the reverse propagation radio frequency wave amplitude.

[0013] Based on the radio frequency voltage distribution and the optical wave time-domain coupling mode equation, the first optical wave amplitude and the second optical wave amplitude in the modulator are obtained;

[0014] Based on the first optical wave amplitude and the second optical wave amplitude, the output optical signal of the modulator is obtained, and a simulated eye diagram is obtained based on the optical signal.

[0015] Optionally, before inputting the radio frequency signal into the radio frequency time-domain coupled mode equation to obtain the forward propagation radio frequency wave amplitude and the reverse propagation radio frequency wave amplitude within the modulator, the method further includes:

[0016] Obtain the initial radio frequency equation.

[0017] Determine the nonlinear coefficients and boundary conditions of the initial radio frequency equation.

[0018] The initial radio frequency equations that determine the nonlinear coefficients are used as the radio frequency time-domain coupling mode equations;

[0019] Before obtaining the first and second optical wave amplitudes within the modulator based on the radio frequency voltage distribution and the optical wave time-domain coupling mode equation, the method further includes:

[0020] Obtain the initial light wave equation.

[0021] Determine the nonlinear coefficients and boundary conditions in the initial light wave equation.

[0022] The initial light wave equation with determined nonlinear coefficients is used as the light wave time-domain coupled mode equation.

[0023] Optionally, the initial radio frequency equation is:

[0024]

[0025]

[0026] V rf,F(z,t)=A rf,F (z,t)exp(iβ rf z-iω rf t)+cc

[0027] t rf,R (z,t)=A rf,R (z,t)exp(-iβ rf z-iω rf t)+cc

[0028] Where, α rf A is the radio frequency wave attenuation constant. rf,F (z,t) represents the amplitude of the forward-propagating radio frequency wave, v g,rf For radio frequency group velocity, A rf,R (z,t) represents the amplitude of the reverse-propagating radio frequency wave, V dc V is the electrostatic field voltage. rf,F (z,t) represents the voltage of the forward-propagating radio frequency wave, V rf,R (z,t) represents the voltage of the reverse propagating radio frequency wave, β rf Let z be the phase constant of the radio frequency wave, z be the spatial coordinate, t be the time, and ω be the phase constant of the radio frequency wave. rf Here, ωc represents the angular frequency of the radio frequency wave, and cc indicates that the preceding expression is taken as its complex conjugate. For a positive RF wave, the second-order nonlinearity, For a positive RF wave, the third-order nonlinearity, For a positive RF wave, the fourth-order nonlinearity, The second-order nonlinearity of the reflected RF wave, The third-order nonlinearity of the reflected RF wave, For the fourth-order nonlinearity of the reflected RF wave, Λ dr,12 Λ r,4 Λ dr,22 Λ dr,22 Λ dr,14 and Λ dr,32 is a nonlinear coefficient.

[0029] Optionally, the initial light wave equation is:

[0030]

[0031]

[0032] Where, α opt Let v be the light wave attenuation constant. g,opt Let A be the group velocity of light. opt (z,t) represents the amplitude of the light wave. This is a second-order nonlinearity of light waves. For the third-order nonlinearity of light waves, V dcV is the electrostatic field voltage. rf,F For forward propagation of radio frequency wave voltage, V rf,R Let ε0 be the voltage of the reverse propagating radio frequency wave, ε0 be the vacuum permittivity, and h be... The power exponent, j is The power exponent, k is The power exponent, For electrostatic field distribution, For radio frequency mode field distribution, For the optical wave mode field distribution, Λ dp,12 Λ rp,11 Λ dp,22 Λ rdp,111 and Λ rp,12 These are nonlinear coefficients.

[0033] Optionally, the nonlinear coefficients and boundary conditions of the initial radio frequency equation are determined, and the determination of the nonlinear coefficients and boundary conditions in the initial optical wave equation is carried out in the following manner:

[0034]

[0035] Γ dp,12 =Λ rp,12

[0036] Λ dp,22 =Λ rp,22 =Λ drp,112

[0037] Optionally, the boundary conditions are configured as follows:

[0038] A rf,R (L,t)=Γ L A rf,F (L,t)exp(2iβ rf L)

[0039]

[0040]

[0041] Λ drp,hjk These are nonlinear coefficients. For electrostatic field distribution, For radio frequency mode field distribution, Let h be the mode field distribution of the light wave. The power exponent, j is The power exponent, k is The exponent of the power, where h, j, and k are positive integers and h + j + k - 1 ≥ 2, and d is used to refer to r is used to refer to p is used to refer to χ is the nonlinear polarizability of the material, Δneff,rf k is the effective refractive index of radio frequency waves. 0,rf Let A be the vacuum wave vector of the radio frequency wave. rf,R (L,t) represents A on the coordinate (L,t). rf,R A rf,F (L,t) represents A on the coordinate (L,t). rf,F ,β rf Let ω be the phase constant of the radio frequency wave, L be the length of the traveling wave region of the modulator, and ω be the phase constant of the radio frequency wave. rf V is the angular frequency of the radio frequency wave. in (t) represents the radio frequency voltage signal input to the modulator, Z c R is the equivalent radio frequency waveguide impedance of the electrode. g R is the resistance value of the internal resistance of the equivalent radio frequency waveguide of the electrode. L Λ is the resistance value of the load resistance of the equivalent radio frequency waveguide of the electrode. dp,12 Λ rp,12 Λ dp,22 Λ rdp,112 Λ rp,p2 Λ drp,hjk Λ dr,12 Λ dr,22 and Λ dr,32 ε is a nonlinear coefficient, and ε0 is the vacuum permittivity.

[0042] Optionally, the step of inputting the radio frequency signal into the radio frequency time-domain coupled mode equation to obtain the forward propagation radio frequency wave amplitude and the reverse propagation radio frequency wave amplitude within the modulator is performed in the following manner:

[0043]

[0044] Δz rf =Δt / β 1,rf

[0045] Δz opt =Δt / β 1,opt

[0046]

[0047] Δt represents the time interval after discretization of time t, and α rf V is the radio frequency wave attenuation constant. dc Let β be the electrostatic field voltage. 1,rf and β 1,opt v is the first derivative of the propagation constant. g,rf v is the group velocity of the RF. g,opt Let Λ be the group velocity of the light wave. dr,12 Λ r,4 Λ dr,22 Λ dr,32 and Λdr,14 These are nonlinear coefficients.

[0048] Optionally, the first and second optical wave amplitudes within the modulator are obtained based on the radio frequency voltage distribution and the optical wave time-domain coupling mode equation, in the following manner:

[0049] A opt，1 (z+Δz opt ,t+Δt)=exp(Δφ1(z,t))A opt，1 (z, t)

[0050] A opt，2 (z+Δz opt ,t+Δt)=exp(Δφ2(z,t))A opt，2 (z, t)

[0051] Δφ1(z,t)=-α opt Δz opt +2iΛ dp，12 (V rf，F (z,t)+V rf，R (z,t)+V dc )Δz opt +3iΛ dp，22 (V rf，F (z,t)+V rf，R (z,t)+V dc ) 2 Δz opt

[0052] Δφ2(z,t)=-α opt Δz opt +2iΛ dp，12 (-V rf，F (z,t)-V rf，R (z,t)+V dc )Δz opt

[0053] Among them, A opt，1 (z, t) represents the amplitude of the first light wave, A opt，2 (z, t) represents the amplitude of the second light wave, Δz opt To represent the interval of discretization of spatial coordinate z, α opt V is the light wave attenuation constant. rf，F (z, t) represents the voltage of the forward-propagating radio frequency wave, V rf，R (z, t) represents the voltage of the reverse propagating radio frequency wave, V dc For electrostatic field voltage, Λ dp，12 Λ dp，22 and Λ dp，12 These are nonlinear coefficients.

[0054] Secondly, a modulator-optical co-simulation device is provided, wherein the modulator receives radio frequency signals and outputs optical signals, and the device includes:

[0055] A signal acquisition unit is used to acquire the radio frequency signal input to the modulator;

[0056] The radio frequency simulation unit is used to input the radio frequency signal into the radio frequency time-domain coupled mode equation to obtain the forward propagation radio frequency wave amplitude and the reverse propagation radio frequency wave amplitude in the modulator;

[0057] The radio frequency simulation unit is also used to obtain the radio frequency voltage distribution in the modulator based on the forward propagation radio frequency wave amplitude and the reverse propagation radio frequency wave amplitude;

[0058] The optical wave simulation unit is used to obtain the first optical wave amplitude and the second optical wave amplitude in the modulator based on the radio frequency voltage distribution and the optical wave time-domain coupling mode equation.

[0059] An eye diagram generation unit is used to obtain the optical signal output by the modulator based on the first optical wave amplitude and the second optical wave amplitude, and to obtain a simulated eye diagram based on the optical signal.

[0060] Thirdly, a computer device is provided, including a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor, when executing the computer program, implements the method as described in any of the preceding claims.

[0061] Fourthly, a computer-readable storage medium is provided having a computer program stored thereon, which, when executed by a processor, implements the method as described in any of the preceding claims.

[0062] This invention provides a modulator optoelectronic co-simulation device, computer equipment, and storage medium. The modulator receives a radio frequency (RF) signal and outputs an optical signal. The method includes: acquiring the RF signal input to the modulator; inputting the RF signal into an RF time-domain coupled-mode equation to acquire the forward propagation RF wave amplitude and the reverse propagation RF wave amplitude within the modulator; acquiring the RF voltage distribution within the modulator based on the forward propagation RF wave amplitude and the reverse propagation RF wave amplitude; acquiring the first optical wave amplitude and the second optical wave amplitude within the modulator based on the RF voltage distribution and the optical wave time-domain coupled-mode equation; acquiring the output optical signal of the modulator based on the first optical wave amplitude and the second optical wave amplitude; and obtaining a simulated eye diagram based on the optical signal. The method of this invention simultaneously achieves electrical and optical simulations in the modulator design process, thereby enabling efficient and accurate design and simulation of the modulator, promoting the development of optical communication and optical computing systems. Furthermore, the method of this invention does not rely on specific materials or require the establishment of specific circuit models, thus having wider applicability and higher efficiency. Attached Figure Description

[0063] The accompanying drawings, which are incorporated in and form part of this specification, illustrate embodiments consistent with the invention and, together with the description, serve to explain the principles of the invention.

[0064] To more clearly illustrate the technical solutions in the embodiments of the present invention or the prior art, the drawings used in the description of the embodiments or the prior art will be briefly introduced below. Obviously, those skilled in the art can obtain other drawings based on these drawings without creative effort.

[0065] Figure 1 The diagram shows the application environment of the modulator optoelectronic co-simulation method according to an embodiment of the present invention.

[0066] Figure 2 The diagram shown is a schematic flowchart of the modulator optoelectronic co-simulation method according to an embodiment of the present invention.

[0067] Figure 3 The diagram shown is a schematic representation of a radio frequency circuit model according to an embodiment of the present invention.

[0068] Figure 4 The figure shows the calculation A according to an embodiment of the present invention. rf,F and A rf,R Numerical calculation process of the time-domain coupled mode equations of radio frequency waves;

[0069] Figure 5 The figure shows the calculation A according to an embodiment of the present invention. opt Numerical calculation process of the optical wave temporal coupled mode equations;

[0070] Figure 6 The image shows a comparison between the eye diagram obtained by the simulation method of this invention and the eye diagram measured experimentally.

[0071] Figure 7 The diagram shown is a structural block diagram of the modulator optoelectronic co-simulation device according to an embodiment of the present invention.

[0072] Figure 8 This is an internal structural diagram of a computer device in an embodiment of the present invention. Detailed Implementation

[0073] To make the objectives, technical solutions, and advantages of the embodiments of the present invention clearer, the technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of the present invention, not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention.

[0074] Figure 1 This is a diagram illustrating the application environment of the modulator-optoelectronic co-simulation method in one embodiment. (Refer to...) Figure 1 This modulator optoelectronic co-simulation method is applied to a modulator optoelectronic co-simulation system. The method includes a terminal 110 and / or a server 120. The terminal 110 and server 120 are connected via a network. The terminal 110 can be a desktop terminal or a mobile terminal; the mobile terminal can be at least one of a mobile phone, tablet computer, or laptop computer. The server 120 can be a standalone server or a server cluster consisting of multiple servers.

[0075] like Figure 2 As shown, in one embodiment, a modulator optoelectronic co-simulation method is provided. This embodiment mainly applies this method to the above-mentioned... Figure 1 Taking server 120 as an example, the modulator receives radio frequency signals and outputs optical signals. The method of this embodiment is used to simulate the output optical signal of the modulator.

[0076] Reference Figure 2 The modulator optoelectronic co-simulation method includes:

[0077] Step 210: Obtain the radio frequency signal input to the modulator;

[0078] Step 220: Input the radio frequency signal into the radio frequency time-domain coupled mode equation to obtain the forward propagation radio frequency wave amplitude and the reverse propagation radio frequency wave amplitude in the modulator;

[0079] Step 230: Obtain the radio frequency voltage distribution within the modulator based on the forward propagation radio frequency wave amplitude and the reverse propagation radio frequency wave amplitude;

[0080] Step 240: Based on the radio frequency voltage distribution and the optical wave time-domain coupling mode equation, obtain the first optical wave amplitude and the second optical wave amplitude in the modulator;

[0081] Step 250: Based on the first optical wave amplitude and the second optical wave amplitude, obtain the optical signal output by the modulator, and obtain a simulated eye diagram based on the optical signal.

[0082] In this embodiment of the invention, time-coupled-mode equations (TCMEs) include radio frequency time-coupled-mode equations and optical wave time-coupled-mode equations.

[0083] Modulators typically consist of two parts: radio frequency electrodes and optical waveguides. In this embodiment of the invention, the optical waveguide is designed as a Mach-Zehnder interferometer (MZI). Due to the phase difference between the two arms of the MZI, the amplitude of the optical signal can be modulated according to the interference of light. To generate the phase difference, the refractive index of the optical waveguide in both arms needs to be changed. Specifically, in silicon optical modulators, a voltage is applied between the G and S electrodes. The voltage on the PN junction changes the distribution of electrons and holes, and the refractive index of the optical waveguide changes according to the plasmon dispersion effect. These PN junctions are arranged in a laterally push-pull configuration, thereby halving the capacitance contributed by the PN junction. The PN junction is formed by doping N-type and P-type regions in the optical waveguide. The N++ and P++ doped regions achieve electrical connection between the electrodes and the PN junction. Similarly, in lithium niobate modulators, a voltage is applied between the G and S electrodes, and the refractive index of the optical waveguide changes according to the Pockels effect of lithium niobate material. In this embodiment of the invention, the time-domain coupled-mode equation (TCME) is derived from Maxwell's equations to describe the nonlinear effects of radio frequency waves and optical waves. Secondly, the nonlinear coefficients in the TCME are obtained by solving and fitting the eigenmodes of radio frequency waves and optical waves. Finally, the difference form of the TCME on a time-space two-dimensional grid is introduced to realize the numerical calculation of the TCME.

[0084] In the method of this embodiment of the invention, the input radio frequency (RF) signal of the modulator is acquired; the RF signal is input into the RF time-domain coupled mode equation to acquire the forward propagation RF wave amplitude and the reverse propagation RF wave amplitude within the modulator; based on the forward propagation RF wave amplitude and the reverse propagation RF wave amplitude, the RF voltage distribution within the modulator is acquired; based on the RF voltage distribution and the optical wave time-domain coupled mode equation, the first optical wave amplitude and the second optical wave amplitude within the modulator are acquired; based on the first optical wave amplitude and the second optical wave amplitude, the output optical signal of the modulator is acquired; and a simulated eye diagram is obtained based on the optical signal. The method of this embodiment of the invention simultaneously realizes electrical and optical simulations in the modulator design process, thereby enabling efficient and accurate design and simulation of the modulator, promoting the development of optical communication and optical computing systems. Furthermore, the method of this embodiment of the invention does not rely on specific materials or require the establishment of specific circuit models, thus having wider applicability and higher efficiency.

[0085] In this embodiment of the invention, before step 220, before inputting the radio frequency signal into the radio frequency time-domain coupled mode equation to obtain the forward propagation radio frequency wave amplitude and the reverse propagation radio frequency wave amplitude within the modulator, the method further includes:

[0086] Obtain the initial radio frequency equation.

[0087] Determine the nonlinear coefficients and boundary conditions of the initial radio frequency equation.

[0088] The initial radio frequency equations that determine the nonlinear coefficients are used as the radio frequency time-domain coupling mode equations;

[0089] The method of this embodiment of the invention, before step 240, before obtaining the first optical wave amplitude and the second optical wave amplitude in the modulator based on the radio frequency voltage distribution and the optical wave time-domain coupling mode equation, further includes:

[0090] Obtain the initial light wave equation.

[0091] Determine the nonlinear coefficients and boundary conditions in the initial light wave equation.

[0092] The initial light wave equation with determined nonlinear coefficients is used as the light wave time-domain coupled mode equation.

[0093] In this embodiment of the invention, the initial radio frequency equation is:

[0094]

[0095] V rf,F (z,t)=A rf,F (z,t)exp(iβ rf z-iωrf t)+cc

[0096] V rf,R (z,t)=A rf,R (z,t)exp(-iβ rf z-iω rf t)+cc

[0097] Where, α rf A is the radio frequency wave attenuation constant. rf,F (z,t) represents the amplitude of the forward-propagating radio frequency wave, v g，rf For radio frequency group velocity, A rf，R (z, t) represents the amplitude of the reverse-propagating radio frequency wave, V dc V is the electrostatic field voltage. rf，F (z, t) represents the voltage of the forward-propagating radio frequency wave, V rf，R (z, t) represents the voltage of the reverse propagating radio frequency wave, β rf Let z be the phase constant of the radio frequency wave, z be the spatial coordinate, t be the time, and ω be the phase constant of the radio frequency wave. rf Here, ωc represents the angular frequency of the radio frequency wave, and cc indicates that the preceding expression is taken as its complex conjugate. For a positive RF wave, the second-order nonlinearity, For a positive RF wave, the third-order nonlinearity, For a positive RF wave, the fourth-order nonlinearity, The second-order nonlinearity of the reflected RF wave, The third-order nonlinearity of the reflected RF wave, For the fourth-order nonlinearity of the reflected RF wave, Λ dr，12 Λ r，4 Λ dr，22 Λ dr，22 Λ dr，14 and Λ dr，32 These are nonlinear coefficients.

[0098] The electromagnetic field in the modulator can be approximated as consisting only of the static voltage V. dc The resulting static electric field distribution E dc Light wave amplitude A opt Model field distribution E opt The amplitude A of the forward-propagating radio frequency wave rf，F and the model field distribution E rf The amplitude A of the reverse-propagating radio frequency wave rf，R The mode field distribution of the reverse-propagating radio frequency wave is the same as that of the forward-propagating radio frequency wave. In this embodiment of the invention, the time-domain coupled mode equation describing the nonlinear effects within the modulator can be derived from Maxwell's equations. The time-domain coupled mode equation with unconfirmed nonlinear parameters can be the initial radio frequency equation.

[0099] In this embodiment of the invention, the initial light wave equation is:

[0100]

[0101] Where, α opt Let v be the light wave attenuation constant. g，opt Let A be the group velocity of light. opt (z, t) represents the amplitude of the light wave. This is a second-order nonlinearity of light waves. For the third-order nonlinearity of light waves, V dc V is the electrostatic field voltage. rf,F For forward propagation of radio frequency wave voltage, V rf,R Let ε0 be the voltage of the reverse propagating radio frequency wave, ε0 be the vacuum permittivity, and h be... The power exponent, j is The power exponent, k is The power exponent, For electrostatic field distribution, For radio frequency mode field distribution, For the optical wave mode field distribution, Λ dp,12 Λ rp,12 Λ dp,22 Λ rdp,112 and Λ rp,22 These are nonlinear coefficients.

[0102] In this embodiment of the invention, the determination of the nonlinear coefficients and boundary conditions of the initial radio frequency equation, and the determination of the nonlinear coefficients and boundary conditions in the initial optical wave equation, are carried out in the following manner:

[0103]

[0104] Λ dp,12 =Λ rp,12

[0105] Λ dp,22 =Λ rp,22 =Λ drp,112

[0106] The boundary conditions are defined as follows:

[0107] A rf,R (L,t)=Γ L A rf,F (L,t)exp(2iβ rf L)

[0108]

[0109] Λ drp,hjk These are nonlinear coefficients. For electrostatic field distribution, For radio frequency mode field distribution, Let h be the mode field distribution of the light wave. The power exponent, j is The power exponent, k is The exponent of the power, where h, j, and k are positive integers and h + j + k - 1 ≥ 2, and d is used to refer to r is used to refer to p is used to refer to χ is the nonlinear polarizability of the material, Δn eff,rf k is the effective refractive index of radio frequency waves. 0,rf Let A be the vacuum wave vector of the radio frequency wave. rf,R (L,t) represents A on the coordinate (L,t). rf,R A rf,F (L,t) represents A on the coordinate (L,t). rf,F ,β rf Let ω be the phase constant of the radio frequency wave, L be the length of the traveling wave region of the modulator, and ω be the phase constant of the radio frequency wave. rf V is the angular frequency of the radio frequency wave. in (t) represents the radio frequency voltage signal input to the modulator, Z c R is the equivalent radio frequency waveguide impedance of the electrode. g R is the resistance value of the internal resistance of the equivalent radio frequency waveguide of the electrode. L Λ is the resistance value of the load resistance of the equivalent radio frequency waveguide of the electrode. dp,12 Λ rp,12 Λ dp,22 Λ rdp,112 Λ rp,22 Λ drp,hjk Λ dr,12 Λ dr,22 and Λ dr,32 ε is a nonlinear coefficient, and ε0 is the vacuum permittivity.

[0110] The boundary conditions described in the above embodiments of the present invention are only one approach. Other approaches exist in other embodiments of the present invention, which will not be elaborated here.

[0111] In this embodiment of the invention, the nonlinear coefficient can be based on the change Δn of the effective refractive index of the radio frequency wave and the optical wave. eff,rf or Δn eff,opt With static voltage V dc The fitting is obtained, and the specific calculation method is as described above. Compared with calculating through the definition of nonlinear coefficients, this method is simpler and easier to combine with experimental measurements. Furthermore, some materials lack experimentally measured material parameters χ, such as PN junctions. The method of this embodiment can circumvent this problem, thus making it more universal.

[0112] Figure 3 The diagram shown is a schematic representation of a radio frequency circuit model according to an embodiment of the present invention. Figure 3 As shown, the boundary conditions can be derived based on this model. The electrode is considered to have an impedance Z. c A radio frequency waveguide connected to a load resistor R L and an internal resistor R g The signal source. Boundary conditions describe the reflection and transmission of radio frequency waves at the interfaces at both ends of the radio frequency waveguide.

[0113] In this embodiment of the invention, step 220, which involves inputting the radio frequency signal into the radio frequency time-domain coupled mode equation to obtain the forward propagation radio frequency wave amplitude and the reverse propagation radio frequency wave amplitude within the modulator, is performed in the following manner:

[0114]

[0115] Δz rf =Δt / β 1,rf

[0116] Δz opt =Δt / β 1,opt

[0117]

[0118] Δt represents the time interval after discretization of time t, and α rf V is the radio frequency wave attenuation constant. dc Let β be the electrostatic field voltage. 1,rf and β 1,opt v is the first derivative of the propagation constant. g,rf v is the group velocity of the RF. g,opt Let Λ be the group velocity of the light wave. dr,12 Λ r,4 Λ dr,22 Λ dr,32 and Λ dr,14 These are nonlinear coefficients.

[0119] In this embodiment of the invention, step 240, obtaining the first and second optical wave amplitudes within the modulator based on the radio frequency voltage distribution and the optical wave time-domain coupling mode equation, is performed in the following manner:

[0120] A opt，1 (z+Δz opt ,t+Δt)=exp(Δφ1(z,t))A opt，1 (z, t)

[0121] A opt,2 (z+Δz opt ,t+Δt)=exp(Δφ2(z,t))A opt,2 (z, t)

[0122] Δφ1(z,t)=-α opt Δz opt +2iΛ dp，12 (V rf，F (z,t)+V rf，R (z,t)+V dc )Δz opt +3iΛ dp,22 (V rf，F (z,t)+V rf，R (z,t)+V dc ) z Δz opt

[0123] Δφ2(z,t)=-α opt Δz opt +2iΛ dp，12 (-V rf，F (z,t)-V rf，R (z,t)+V dc )Δz opt

[0124] Among them, A opt，1 (z, t) represents the amplitude of the first light wave, A opt,2 (z, t) represents the amplitude of the second light wave, Δz opt To represent the interval of discretization of spatial coordinate z, α opt V is the light wave attenuation constant. rf，F (z, t) represents the voltage of the forward-propagating radio frequency wave, V rf，R (z, t) represents the voltage of the reverse propagating radio frequency wave, V dc For electrostatic field voltage, Λ dp，12 Λ dp,22 and Λ dp，12 These are nonlinear coefficients.

[0125] In this embodiment of the invention, the specific solution based on the TCME relies on a difference form within a two-dimensional spatiotemporal grid. The TCMEs for radio frequency waves and light waves each have their own two-dimensional spatiotemporal grids, where the grid size Δt is the same for both radio frequency waves and light waves. The grid sizes Δz for radio frequency waves and light waves are different. rf =Δt / β 1，rf and Δz opt =Δt / β 1，opt .

[0126] Figure 4 Demonstrates calculation A rf,F and A rf，R The numerical calculation process for radio frequency wave TCME is described. This calculation is based on a spatial interval of ΔZ. rf A two-dimensional spatiotemporal grid with time intervals Δt, where each circle represents a storage A.rf,F and A rf，R The nodes represent values, and all circles are initialized to zero. Black or gray arrows indicate calculations of A. rf，F Or A rf，R The order is as follows. Boundary conditions are indicated by green arrows. Orange arrows mark the input voltage signal V. in The algorithm calculates A rf,F Initially, based on the boundary conditions on the left side of the RF circuit, by V in The excitation. Due to impedance mismatch, according to the boundary conditions on the right side of the RF circuit, A rf,R Then from A rf，F Excitation. In the simulation, Δt = 0.5 ps was used to ensure numerical convergence. Finally, V on the grid was obtained from the forward and backward propagating RF wave voltages. rf，F and V rf，R The distribution of .

[0127] Figure 5 Demonstrates calculation A opt The numerical calculation process for optical waves using TCME. Input optical wave A opt,in It is divided into two paths, namely the two arms of the MZI (Mach-Zehnder interferometer, abbreviated as MZI), namely A opt,1 and A opt,2 The black arrow indicates the calculation of A. opt,1 and A opt,2 The order is shown, and the dashed green arrows mark the nonlinear terms in the light wave TCME. Figure 5 V in rf V represents rf,F V rf,R and V dc Due to Δz opt ≠Δz rf We use interpolation to resample V along z = 0 to z = L at t. rf Distribution. Considering that optical waves oscillate faster than radio frequency waves, we do not use a finer grid size Δz. opt Instead, A opt,1 and A opt,2 The differential form is transformed into the exponential differential form. The final light wave output by the modulator is obtained through coherent superposition of A... opt,1 and A opt,2 Get A opt,out .

[0128] In this embodiment of the invention, the method for obtaining the eye diagram can be any method in the prior art, and will not be described in detail here.

[0129] Figure 6The image shows a comparison between the eye diagram obtained by the simulation method of this invention and the experimentally measured eye diagram. It can be seen that the eye diagram obtained by the simulation method of this invention is as follows: Figure 6 (a) Eye diagram as measured experimentally Figure 6 (b) Very close.

[0130] Figure 6 In the eye diagram, the red and blue colors are due to device or system settings and do not affect the comparison structure of the eye diagram.

[0131] The method of this invention can simultaneously perform electrical and optical simulations, thereby enabling efficient and accurate design and simulation of modulators, promoting the development of optical communication and optical computing systems. Furthermore, in this method, the RF electrodes are treated as RF waveguides rather than equivalent circuit models; the silicon PN junction is considered an equivalent fourth-order nonlinear medium; and the nonlinear IV relationship of the PN junction is replaced by self-phase modulation (SPM) and cross-phase modulation (XPM) of the RF wave. In this method, the dispersion of the RF electrodes is automatically resolved by using the group velocity of the RF wave in the RF waveguide, avoiding the convolution calculations required in previous methods.

[0132] like Figure 7 As shown, the present invention also provides a modulator optoelectronic co-simulation device, the device comprising:

[0133] The signal acquisition unit 710 is used to acquire the radio frequency signal input to the modulator;

[0134] The radio frequency simulation unit 720 is used to input the radio frequency signal into the radio frequency time-domain coupled mode equation to obtain the forward propagation radio frequency wave amplitude and the reverse propagation radio frequency wave amplitude in the modulator;

[0135] The radio frequency simulation unit 720 is also used to obtain the radio frequency voltage distribution in the modulator based on the forward propagation radio frequency wave amplitude and the reverse propagation radio frequency wave amplitude;

[0136] The optical wave simulation unit 730 is used to obtain the first optical wave amplitude and the second optical wave amplitude in the modulator based on the radio frequency voltage distribution and the optical wave time-domain coupling mode equation.

[0137] The eye diagram generation unit 740 is used to obtain the optical signal output by the modulator based on the first optical wave amplitude and the second optical wave amplitude, and to obtain a simulated eye diagram based on the optical signal.

[0138] In this embodiment of the invention, the radio frequency simulation unit 720 is further configured to:

[0139] Obtain the initial radio frequency equation.

[0140] Determine the nonlinear coefficients and boundary conditions of the initial radio frequency equation.

[0141] The initial radio frequency equations for determining the nonlinear coefficients are used as the radio frequency time-domain coupled mode equations.

[0142] In this embodiment of the invention, the initial radio frequency equation is:

[0143]

[0144] V rf,F (z,t)=A rf,F (z,t)exp(iβ rf z-iω rf t)+cc

[0145] V rf,R (z,t)=A rf,R (z,t)exp(-iβ rf z-iω rf t)+cc

[0146] Where, α rf A is the radio frequency wave attenuation constant. rf，F (z, t) represents the amplitude of the forward-propagating radio frequency wave, v g，rf For radio frequency group velocity, A rf，R (z, t) represents the amplitude of the reverse-propagating radio frequency wave, V dc V is the electrostatic field voltage. rf，F (z, t) represents the voltage of the forward-propagating radio frequency wave, V rf，R (z, t) represents the voltage of the reverse propagating radio frequency wave, β rf Let z be the phase constant of the radio frequency wave, z be the spatial coordinate, t be the time, and ω be the phase constant of the radio frequency wave. rf Here, ωc represents the angular frequency of the radio frequency wave, and cc indicates that the preceding expression is taken as its complex conjugate. For a positive RF wave, the second-order nonlinearity, For a positive RF wave, the third-order nonlinearity, For a positive RF wave, the fourth-order nonlinearity, The second-order nonlinearity of the reflected RF wave, The third-order nonlinearity of the reflected RF wave, For the fourth-order nonlinearity of the reflected RF wave, Λ dr，12 Λ r，4 Λ dr，22 Λ dr,22 Λ dr,14 and Λ dr,32 These are nonlinear coefficients.

[0147] In this embodiment of the invention, the optical wave simulation unit 730 is further used for:

[0148] Obtain the initial light wave equation;

[0149] Determine the nonlinear coefficients and boundary conditions in the initial light wave equation;

[0150] The initial light wave equation with determined nonlinear coefficients is used as the light wave time-domain coupled mode equation.

[0151] In this embodiment of the invention, the initial light wave equation is:

[0152]

[0153] Where, α opt Let v be the light wave attenuation constant. g，opt Let A be the group velocity of light. opt (z, t) represents the amplitude of the light wave. This is a second-order nonlinearity of light waves. For the third-order nonlinearity of light waves, V dc V is the electrostatic field voltage. rf，F For forward propagation of radio frequency wave voltage, V rf，R Let ε0 be the voltage of the reverse propagating radio frequency wave, ε0 be the vacuum permittivity, and h be... The power exponent, j is The power exponent, k is The power exponent, For electrostatic field distribution, For radio frequency mode field distribution, For the optical wave mode field distribution, Λ dp,12 Λ rp,12 Λ dp,22 Λ rdp,112 and Λ rp,22 These are nonlinear coefficients.

[0154] In this embodiment of the invention, the determination of the nonlinear coefficients and boundary conditions of the initial radio frequency equation, and the determination of the nonlinear coefficients and boundary conditions in the initial optical wave equation, are carried out in the following manner:

[0155]

[0156] Λ dp,12 =Λ rp,12

[0157] Λ dp,22 =Λ rp,22 =Λ drp,112

[0158] The boundary conditions are defined as follows:

[0159] A rf,R (L,t)=Γ L A rf,F(L,t)exp(2iβ rf L)

[0160]

[0161] Λ drp,hjk These are nonlinear coefficients. For electrostatic field distribution, For radio frequency mode field distribution, Let h be the mode field distribution of the light wave. The power exponent, j is The power exponent, k is The exponent of the power, where h, j, and k are positive integers and h + j + k - 1 ≥ 2, and d is used to refer to r is used to refer to p is used to refer to χ is the nonlinear polarizability of the material, Δn eff，rf k is the effective refractive index of radio frequency waves. 0，rf Let A be the vacuum wave vector of the radio frequency wave. rf，R (L, t) represents A on the coordinate (L, t). rf，R A rf，F (L, t) represents A on the coordinate (L, t). rf，F ,β rf Let ω be the phase constant of the radio frequency wave, L be the length of the traveling wave region of the modulator, and ω be the phase constant of the radio frequency wave. rf V is the angular frequency of the radio frequency wave. in (t) represents the radio frequency voltage signal input to the modulator, Z c R is the equivalent radio frequency waveguide impedance of the electrode. g R is the resistance value of the internal resistance of the equivalent radio frequency waveguide of the electrode. L Λ is the resistance value of the load resistance of the equivalent radio frequency waveguide of the electrode. dp，12 Λ rp，12 Λ dp，22 Λ rdp，112 Λ rp，22 Λ drp，hjk Λ dr，12 Λ dr，22 and Λ dr.32 ε is a nonlinear coefficient, and ε0 is the vacuum permittivity.

[0162] In this embodiment of the invention, the radio frequency simulation unit 720 is further configured to obtain the forward propagation radio frequency wave amplitude and the reverse propagation radio frequency wave amplitude in the following manner:

[0163]

[0164] Δz rf =Δt / β 1，rf

[0165] Δz opt =Δt / β 1，opt

[0166]

[0167] Δt represents the time interval after discretization of time t, and α rf V is the radio frequency wave attenuation constant. dc Let β be the electrostatic field voltage. 1，rf and β 1，opt v is the first derivative of the propagation constant. g，rf v is the group velocity of the RF. g，opt Let Λ be the group velocity of the light wave. dr，12 Λ r，4 Λ dr,22 A dr，32 and Λ dr，14 These are nonlinear coefficients.

[0168] In this embodiment of the invention, the optical wave simulation unit 730 is further configured to obtain the first optical wave amplitude and the second optical wave amplitude in the following manner:

[0169] A opt，1 (z+Δz opt ,t+Δt)=exp(Δφ1(z,t))A opt，1 (z, t)

[0170] A opt,2 (z+Δz opt ,t+Δt)=exp(Δφ2(z,t))A opt,2 (z, t)

[0171] Δφ1(z,t)=-α opt Δz opt +2iΛ dp，12 (V rf，F (z,t)+V rf，R (z,t)+V dc )Δz opt +3iΛ dp,22 (V rf，F (z,t)+V rf，R (z,t)+V dc ) 2 Δz opt

[0172] Δφ2(z,t)=-α opt Δz opt +2iΛ dp，12 (-V rf，F (z,t)-V rf，R (z,t)+V dc )Δz opt

[0173] Among them, A opt，1 (z, t) represents the amplitude of the first light wave, A opt,2 (z, t) represents the amplitude of the second light wave, Δz opt To represent the interval of discretization of spatial coordinate z, α opt V is the light wave attenuation constant. rf，F (z, t) represents the voltage of the forward-propagating radio frequency wave, V rf，R (z, t) represents the voltage of the reverse propagating radio frequency wave, V dc For electrostatic field voltage, Λ dp,12 Λ dp,22 and Λ dp,12 These are nonlinear coefficients.

[0174] The method of this invention simultaneously performs electrical and optical simulations during the modulator design process, thereby enabling efficient and accurate design and simulation of the modulator and promoting the development of optical communication and optical computing systems. Furthermore, the method of this invention does not rely on specific materials or require the establishment of specific circuit models, thus possessing wider applicability and higher efficiency.

[0175] This invention also provides a computer device, including a memory, a processor, and a computer program stored in the memory and executable on the processor. When the processor executes the computer program, it implements the following method: acquiring an input radio frequency (RF) signal from the modulator; inputting the RF signal into an RF time-domain coupled-mode equation to acquire the forward propagation RF wave amplitude and the reverse propagation RF wave amplitude within the modulator; acquiring the RF voltage distribution within the modulator based on the forward propagation RF wave amplitude and the reverse propagation RF wave amplitude; acquiring a first optical wave amplitude and a second optical wave amplitude within the modulator based on the RF voltage distribution and the optical wave time-domain coupled-mode equation; acquiring the output optical signal from the modulator based on the first optical wave amplitude and the second optical wave amplitude; and obtaining a simulated eye diagram based on the optical signal.

[0176] This invention also provides a computer-readable storage medium storing a computer program thereon. When executed by a processor, the computer program implements the following method: acquiring the input radio frequency (RF) signal of the modulator; inputting the RF signal into the RF time-domain coupled-mode equation to acquire the forward propagation RF wave amplitude and the reverse propagation RF wave amplitude within the modulator; acquiring the RF voltage distribution within the modulator based on the forward propagation RF wave amplitude and the reverse propagation RF wave amplitude; acquiring the first optical wave amplitude and the second optical wave amplitude within the modulator based on the RF voltage distribution and the optical wave time-domain coupled-mode equation; acquiring the output optical signal of the modulator based on the first optical wave amplitude and the second optical wave amplitude; and obtaining a simulated eye diagram based on the optical signal.

[0177] The aforementioned modulator optoelectronic co-simulation method achieves the beneficial effect of solving the technical problems raised in the background art.

[0178] Figure 2 This is a flowchart illustrating the modulator-optoelectronic co-simulation method in one embodiment. It should be understood that, although... Figure 2 The steps in the flowchart are shown sequentially as indicated by the arrows, but these steps are not necessarily executed in the order indicated by the arrows. Unless otherwise specified herein, there is no strict order in which these steps are executed, and they can be performed in other orders. Figure 2 At least some of the steps in the process may include multiple sub-steps or multiple stages. These sub-steps or stages are not necessarily completed at the same time, but can be executed at different times. The execution order of these sub-steps or stages is not necessarily sequential, but can be executed in turn or alternately with other steps or at least some of the sub-steps or stages of other steps.

[0179] Figure 8 An internal structural diagram of a computer device in one embodiment is shown. Specifically, this computer device may be... Figure 1 Server 120 in the middle. For example... Figure 8As shown, the computer device includes a processor, memory, network interface, input device, and display screen connected via a system bus. The memory includes non-volatile storage media and internal memory. The non-volatile storage media stores an operating system and may also store a computer program. When executed by the processor, this computer program enables the processor to implement a modulator-optical co-simulation method. The internal memory may also store a computer program, which, when executed by the processor, enables the processor to execute the modulator-optical co-simulation method. The display screen can be a liquid crystal display (LCD) or an electronic ink display. The input device can be a touch layer covering the display screen, buttons, a trackball, or a touchpad mounted on the computer device's casing, or an external keyboard, touchpad, or mouse.

[0180] Those skilled in the art will understand that Figure 8 The structure shown is merely a block diagram of a portion of the structure related to the present invention and does not constitute a limitation on the computer device to which the present invention is applied. A specific computer device may include more or fewer components than those shown in the figure, or combine certain components, or have different component arrangements.

[0181] Those skilled in the art will understand that all or part of the processes in the methods of the above embodiments can be implemented by a computer program instructing related hardware. The program can be stored in a non-volatile computer-readable storage medium, and when executed, it can include the processes of the embodiments of the above methods. Any references to memory, storage, databases, or other media used in the embodiments provided by this invention can include non-volatile and / or volatile memory. Non-volatile memory can include read-only memory (ROM), programmable ROM (PROM), electrically programmable ROM (EPROM), electrically erasable programmable ROM (EEPROM), or flash memory. Volatile memory can include random access memory (RAM) or external cache memory. By way of illustration and not limitation, RAM is available in various forms, such as static RAM (SRAM), dynamic RAM (DRAM), synchronous DRAM (SDRAM), dual data rate SDRAM (DDRSDRAM), enhanced SDRAM (ESDRAM), synchronous link DRAM (SLDRAM), RAMbus direct RAM (RDRAM), direct memory bus dynamic RAM (DRDRAM), and RAMbus dynamic RAM (RDRAM), etc.

[0182] It should be noted that, in this document, relational terms such as "first" and "second" are used merely to distinguish one entity or operation from another, and do not necessarily require or imply any such actual relationship or order between these entities or operations. Furthermore, the terms "comprising," "including," or any other variations thereof are intended to cover non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements includes not only those elements but also other elements not expressly listed, or elements inherent to such a process, method, article, or apparatus. Without further limitations, an element defined by the phrase "comprising one..." does not exclude the presence of other identical elements in the process, method, article, or apparatus that includes said element.

[0183] The above description is merely a specific embodiment of the present invention, enabling those skilled in the art to understand or implement the invention. Various modifications to these embodiments will be readily apparent to those skilled in the art, and the general principles defined herein may be implemented in other embodiments without departing from the spirit or scope of the invention. Therefore, the present invention is not to be limited to the embodiments shown herein, but is to be accorded the widest scope consistent with the principles and novel features claimed herein.

Claims

1. A modulator optoelectronic co-simulation method, characterized in that, The modulator receives a radio frequency signal and outputs an optical signal; the method includes: Acquire the radio frequency signal input to the modulator; The radio frequency signal is input into the radio frequency time-domain coupled mode equation to obtain the forward propagation radio frequency wave amplitude and the reverse propagation radio frequency wave amplitude in the modulator; The radio frequency voltage distribution within the modulator is obtained based on the forward propagation radio frequency wave amplitude and the reverse propagation radio frequency wave amplitude. Based on the radio frequency voltage distribution and the optical wave time-domain coupling mode equation, the first optical wave amplitude and the second optical wave amplitude in the modulator are obtained; Based on the first optical wave amplitude and the second optical wave amplitude, the output optical signal of the modulator is obtained, and a simulated eye diagram is obtained based on the optical signal; Before inputting the radio frequency signal into the radio frequency time-domain coupled mode equation to obtain the forward propagation radio frequency wave amplitude and the reverse propagation radio frequency wave amplitude within the modulator, the method further includes: Obtain the initial radio frequency equation. Determine the nonlinear coefficients and boundary conditions of the initial radio frequency equation. The initial radio frequency equations that determine the nonlinear coefficients are used as the radio frequency time-domain coupling mode equations; Before obtaining the first and second optical wave amplitudes within the modulator based on the radio frequency voltage distribution and the optical wave time-domain coupling mode equation, the method further includes: Obtain the initial light wave equation. Determine the nonlinear coefficients and boundary conditions in the initial light wave equation. The initial light wave equation with determined nonlinear coefficients is used as the light wave time-domain coupling mode equation; The initial radio frequency equation is: Where, α rf A is the radio frequency wave attenuation constant. rf,F (z,t) represents the amplitude of the forward-propagating radio frequency wave, v g,rf For radio frequency group velocity, A rf,R (z,t) represents the amplitude of the reverse-propagating radio frequency wave, V dc V is the electrostatic field voltage. rf,F (z,t) represents the voltage of the forward-propagating radio frequency wave, V rf,R (z,t) represents the voltage of the reverse propagating radio frequency wave, β rf Let z be the phase constant of the radio frequency wave, z be the spatial coordinate, t be the time, and ω be the phase constant of the radio frequency wave. rf Here, ωc represents the angular frequency of the radio frequency wave, and cc indicates that the preceding expression is taken as its complex conjugate. For a positive RF wave, the second-order nonlinearity, For a positive RF wave, the third-order nonlinearity, For a positive RF wave, the fourth-order nonlinearity, The second-order nonlinearity of the reflected RF wave, The third-order nonlinearity of the reflected RF wave, For the fourth-order nonlinearity of the reflected RF wave, Λ dr,12 Λ r,4 Λ dr,22 Λ dr,22 Λ dr,14 and Λ dr,32 These are nonlinear coefficients; The initial light wave equation is: Where, α opt Let v be the light wave attenuation constant. g,opt Let A be the group velocity of light. opt (z,t) represents the amplitude of the light wave. This is a second-order nonlinearity of light waves. For the third-order nonlinearity of light waves, V dc V is the electrostatic field voltage. rf,F For forward propagation of radio frequency wave voltage, V rf,R Let ε0 be the voltage of the reverse propagating radio frequency wave, ε0 be the vacuum permittivity, and h be... The power exponent, j is The power exponent, k is The power exponent, For electrostatic field distribution, For radio frequency mode field distribution, For the optical wave mode field distribution, Λ dp,12 Λ rp,12 Λ dp,22 Λ rdp,112 and Λ rp,22 These are nonlinear coefficients.

2. The method according to claim 1, characterized in that, The determination of the nonlinear coefficients and boundary conditions of the initial radio frequency equation, and the determination of the nonlinear coefficients and boundary conditions in the initial optical wave equation, are carried out in the following manner: L dp,12 =L rp,12 L dp,22 =L tp,22 =L drp,112 The boundary conditions are defined as follows: From rf,R (L,t)=Γ L From rf,F (L,t)exp(2iβ rf L) Λ drp,hjk These are nonlinear coefficients. For electrostatic field distribution, For radio frequency mode field distribution, Let h be the mode field distribution of the light wave. The power exponent, j is The power exponent, k is The exponent of the power, where h, j, and k are positive integers and h + j + k - 1 ≥ 2, and d is used to refer to r is used to refer to p is used to refer to χ is the nonlinear polarizability of the material, Δn eff,rf k is the effective refractive index of radio frequency waves. 0,rf Let A be the vacuum wave vector of the radio frequency wave. rf,R (L,t) represents A on the coordinate (L,t). rf,R A rf,F (L,t) represents A on the coordinate (L,t). rf,F ,β rf Let ω be the phase constant of the radio frequency wave, L be the length of the traveling wave region of the modulator, and ω be the phase constant of the radio frequency wave. rf V is the angular frequency of the radio frequency wave. in (t) represents the radio frequency voltage signal input to the modulator, Z c R is the equivalent radio frequency waveguide impedance of the electrode. g R is the resistance value of the internal resistance of the equivalent radio frequency waveguide of the electrode. L Λ is the resistance value of the load resistance of the equivalent radio frequency waveguide of the electrode. dp,12 Λ rp,12 Λ dp,22 Λ rdp,112 Λ rp,22 Λ drp,hjk Λ dr,12 Λ dr,22 and Λ dr,32 ε is a nonlinear coefficient, and ε0 is the vacuum permittivity.

3. The method according to claim 1, characterized in that, The step of inputting the radio frequency signal into the radio frequency time-domain coupled mode equation to obtain the forward propagation radio frequency wave amplitude and the reverse propagation radio frequency wave amplitude within the modulator is performed in the following manner: Δz rf =Δt / β 1,rf Δz opt =Δt / β 1,opt Δt represents the time interval after discretization with respect to time t, and α rf V is the radio frequency wave attenuation constant. dc Let β be the electrostatic field voltage. 1,rf and β 1,opt v is the first derivative of the propagation constant. g,rf v is the group velocity of the RF. g,opt Let Λ be the group velocity of the light wave. dr,12 Λ r,4 Λ dr,22 Λ dr,32 and Λ dr,14 These are nonlinear coefficients.

4. The method according to claim 1, characterized in that, The first and second optical wave amplitudes within the modulator are obtained based on the radio frequency voltage distribution and the optical wave time-domain coupling mode equation, in the following manner: A opt,1 (z+Δz opt ,t+Δt)=exp(Δφ1(z,t))A opt,1 (z,t) A opt,2 (z+Δz opt ,t+Δt)=exp(Δφ2(z,t))A opt,2 (z,t) Δφ1(z,t)=-α opt Δz opt +2iΛ dp,12 (V rf,F (z,t)+V rf,R (z,t)+V dc )Δz opt +3iΛ dp,22 (V rf,F (z,t)+V rf,R (z,t)+V dc ) 2 Δz opt Δφ2(z,t)=-α opt Δz opt +2iΛ dp,12 (-V rf,F (z,t)-V rf,R (z,t)+V dc )Δz opt Among them, A opt,1 (z,t) represents the amplitude of the first light wave, A opt,2 (z,t) represents the amplitude of the second light wave, Δz opt To represent the interval of discretization of spatial coordinate z, α opt V is the light wave attenuation constant. rf,F (z,t) represents the voltage of the forward-propagating radio frequency wave, V rf,R (z,t) represents the voltage of the reverse propagating radio frequency wave, V dc For electrostatic field voltage, Λ dp,12 Λ dp,22 and Λ dp,12 These are nonlinear coefficients.

5. A modulator optoelectronic co-simulation device, characterized in that, The modulator receives radio frequency signals and outputs optical signals; the device includes: A signal acquisition unit is used to acquire the radio frequency signal input to the modulator; The radio frequency simulation unit is used to input the radio frequency signal into the radio frequency time-domain coupled mode equation to obtain the forward propagation radio frequency wave amplitude and the reverse propagation radio frequency wave amplitude in the modulator; The radio frequency simulation unit is also used to obtain the radio frequency voltage distribution in the modulator based on the forward propagation radio frequency wave amplitude and the reverse propagation radio frequency wave amplitude; The optical wave simulation unit is used to obtain the first optical wave amplitude and the second optical wave amplitude in the modulator based on the radio frequency voltage distribution and the optical wave time-domain coupling mode equation. An eye diagram generation unit is used to obtain the optical signal output by the modulator based on the first optical wave amplitude and the second optical wave amplitude, and to obtain a simulated eye diagram based on the optical signal. The radio frequency simulation unit is also used for: Obtain the initial radio frequency equation. Determine the nonlinear coefficients and boundary conditions of the initial radio frequency equation. The initial radio frequency equations that determine the nonlinear coefficients are used as the radio frequency time-domain coupling mode equations; The optical wave simulation unit is also used for: Obtain the initial light wave equation. Determine the nonlinear coefficients and boundary conditions in the initial light wave equation. The initial light wave equation with determined nonlinear coefficients is used as the light wave time-domain coupling mode equation; The initial radio frequency equation is: V rf,F (z,t)=A rf,F (z,t)exp(iβ rf z-iω rf t)+c.c. V rf,R (z,t)=A rf,R (z,t)exp(-iβ rf z-iω rf t)+c.c. Where, α rf A is the radio frequency wave attenuation constant. rf,F (z,t) represents the amplitude of the forward-propagating radio frequency wave, v g,rf For radio frequency group velocity, A rf,R (z,t) represents the amplitude of the reverse-propagating radio frequency wave, V dc V is the electrostatic field voltage. rf,F (z,t) represents the voltage of the forward-propagating radio frequency wave, V rf,R (z,t) represents the voltage of the reverse propagating radio frequency wave, β rf Let z be the phase constant of the radio frequency wave, z be the spatial coordinate, t be the time, and ω be the phase constant of the radio frequency wave. rf Here, ωc represents the angular frequency of the radio frequency wave, and cc indicates that the preceding expression is taken as its complex conjugate. For a positive RF wave, the second-order nonlinearity, For a positive RF wave, the third-order nonlinearity, For a positive RF wave, the fourth-order nonlinearity, The second-order nonlinearity of the reflected RF wave, The third-order nonlinearity of the reflected RF wave, For the fourth-order nonlinearity of the reflected RF wave, Λ dr,12 Λ r,4 Λ dr,22 Λ dr,22 Λ dr,14 and Λ dr,32 These are nonlinear coefficients; The initial light wave equation is: Where, α opt Let v be the light wave attenuation constant. g,opt Let A be the group velocity of light. opt (z,t) represents the amplitude of the light wave. This is a second-order nonlinearity of light waves. For the third-order nonlinearity of light waves, V dc V is the electrostatic field voltage. rf,F For forward propagation of radio frequency wave voltage, V rf,R Let ε0 be the voltage of the reverse propagating radio frequency wave, ε0 be the vacuum permittivity, and h be... The power exponent, j is The power exponent, k is The power exponent, For electrostatic field distribution, For radio frequency mode field distribution, For the optical wave mode field distribution, Λ dp,12 Λ rp,12 Λ dp,22 Λ rdp,112 and Λ rp,22 These are nonlinear coefficients.

6. A computer device, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that, When the processor executes the computer program, it implements the method of any one of claims 1 to 4.

7. A computer-readable storage medium having a computer program stored thereon, characterized in that, When the computer program is executed by a processor, it implements the method of any one of claims 1 to 4.

Images