2017.11.11 02:19 *Eugene1026* **iPhone X Wallpapers**

Find gorgeous wallpapers for your shiny new iPhone with OLED display!

2011.03.03 19:14 *boratborat* **Absolutely, Wildly, Unapologetically, SHAMELESS**

2011.09.22 20:12 **House of Lies**

2023.03.27 18:58 *Practical-Mix1818* **Ocd is bothering me and need answers**

I have a cowlick on my hairline that makes it grow a different direction! I pull up my hairline alot just worrying about it. My doctor and my friends all say I'm nuts but I need opinions. Am I balding? Eveyr9ne says I'm nuts but my cowlick side just bothers me. Thanks guys.

submitted by Practical-Mix1818 to Hairloss [link] [comments]
2023.03.27 18:57 *Legendary_Xrayz* **WB Mega Venusaur add me 9904 4533 7926 and 5842 0809 7343**

submitted by Legendary_Xrayz to PokemonGoRaids [link] [comments] |

2023.03.27 18:56 *StepwiseUndrape574* **GTA 5 Online Modded vehicles Shop ⚾ FURYMODZ.com**

submitted by StepwiseUndrape574 to gta5moddedoutfits_ [link] [comments] |

2023.03.27 18:56 *serenity8989* **Four months later, hitting rock bottom again**

I can barely breathe. Function. Survive.

Everyone around me is in these happy relationships and I can barely hear them talk without bursting into tears.

I'm 34. Childless. Going to lose my dog. Losing my home. Losing my husband.

And he blames me for it all. Says I treated him like shit. Gaslit him. His examples...don't feel like they should be divorce worthy but I guess they are to him.

I didn't want to buy a house until we were healthy. He thought that was fucked up. So we bought a house and he resents me for saying I didn't want it.

I once told my sister I thought he took too many selfies. He read my texts and thought that was fucked up cuz I was talking shit about him.

I asked him to stop drinking. I was trying to control him. He could be himself drunk - he told the truth drunk. He said he drank because he was so unhappy in his life with me.

I went out with my friends and went on solo travel trips he didn't want to go on. I was just passionate about travel and we had talked about this? But it was used against me. I wouldn't go if he didn't want me to. I wanted to go with him.

He said I disrespected him. I didnt value him. I took him for granted. How, how, how. I would've fixed it all.

He says we tried counselling and it didn't work. We went to four sessions two years ago! How could you not want to try again.

My mind is a whirl. My therapist obviously only hears my side. Says he's a narcissist and shows BPD tendencies.

All I hear is that my husband, my person doesn't love me anymore. That he left me.

I was doing good four months later. Felt healthy-ish. Thought I'd sign up for a dating site, see what was out there. Haven't been dating in 12 years. Deleted it a week later (wow online dating can be horrific). My ex finds out and although he says he doesn't care, says he won't share custody of the dog if I'm dating. And then launches into a vitriol of the same stuff he said when he was divorcing me.

I treated him like shit. That is his narrative. And it's probably true to him.

I just wanted to be a good wife. I just put up a boundary around the drinking and the spending - that was it. Maybe I did raise my voice sometimes. Maybe I did say mean things during fights. But I don't remember it that way? I want to take accountability? I want to apologize and work on it ? just give me time. Just - fuck. I feel lost.

I truly did and do love him, no matter how much he seethes in anger towards me. And everytime he tells me how much i was the worst wife ever, it's like another knife cut.

submitted by serenity8989 to Divorce [link] [comments]
Everyone around me is in these happy relationships and I can barely hear them talk without bursting into tears.

I'm 34. Childless. Going to lose my dog. Losing my home. Losing my husband.

And he blames me for it all. Says I treated him like shit. Gaslit him. His examples...don't feel like they should be divorce worthy but I guess they are to him.

I didn't want to buy a house until we were healthy. He thought that was fucked up. So we bought a house and he resents me for saying I didn't want it.

I once told my sister I thought he took too many selfies. He read my texts and thought that was fucked up cuz I was talking shit about him.

I asked him to stop drinking. I was trying to control him. He could be himself drunk - he told the truth drunk. He said he drank because he was so unhappy in his life with me.

I went out with my friends and went on solo travel trips he didn't want to go on. I was just passionate about travel and we had talked about this? But it was used against me. I wouldn't go if he didn't want me to. I wanted to go with him.

He said I disrespected him. I didnt value him. I took him for granted. How, how, how. I would've fixed it all.

He says we tried counselling and it didn't work. We went to four sessions two years ago! How could you not want to try again.

My mind is a whirl. My therapist obviously only hears my side. Says he's a narcissist and shows BPD tendencies.

All I hear is that my husband, my person doesn't love me anymore. That he left me.

I was doing good four months later. Felt healthy-ish. Thought I'd sign up for a dating site, see what was out there. Haven't been dating in 12 years. Deleted it a week later (wow online dating can be horrific). My ex finds out and although he says he doesn't care, says he won't share custody of the dog if I'm dating. And then launches into a vitriol of the same stuff he said when he was divorcing me.

I treated him like shit. That is his narrative. And it's probably true to him.

I just wanted to be a good wife. I just put up a boundary around the drinking and the spending - that was it. Maybe I did raise my voice sometimes. Maybe I did say mean things during fights. But I don't remember it that way? I want to take accountability? I want to apologize and work on it ? just give me time. Just - fuck. I feel lost.

I truly did and do love him, no matter how much he seethes in anger towards me. And everytime he tells me how much i was the worst wife ever, it's like another knife cut.

2023.03.27 18:56 *Legendary_Xrayz* **WB Mega Venusaur add me 9904 4533 7926 and 5842 0809 7343 adding first 15**

submitted by Legendary_Xrayz to PokemonGoRaids [link] [comments] |

2023.03.27 18:56 *Many-Rice7351* **Hahaha Hannah hasn’t even dated 5 women 🤣**

submitted by Many-Rice7351 to kaylaannxoxosnark [link] [comments] |

2023.03.27 18:56 *Legendary_Xrayz* **WB Mega Venusaur add me 9904 4533 7926 and 5842 0809 7343**

submitted by Legendary_Xrayz to PokemonGoRaids [link] [comments] |

2023.03.27 18:55 *VladyComan* **KITGG vs JJSPORT vs FC24SHOPT: in hand review**

I m not making reviews at all, this is my honest, transparent review about these 3 sites which i orderd almost on the same day. I ordered the same 3 jersey with same specs on every on of these 3 sites. Kitgg: 78,7$ - jjsport: 61$ - fc24shop: 73,8$ ->final prices on 3 shirts, shipping included Kitgg: ordered on 08 march, on 13 march the order was on the way, received at my house on 27 march; Jjsport: ordered on 07 march, on 17 march the order was on the way, received at my house on 27 march; Fc24shop: ordered on 07 march, on 12 march the order was on the way, received at my house on 24 march I ll write them down 1 by 1 and on the final of the post i ll share with you guyz my conclusion.

KITTGG: stitching almost perfect, no flaws on the shirts at all, very very good materials on every shirts, patches all rigtht, like the perfect order

JJSPORT: 1. the “Spotify” from the barca shirt is too high on the front side, too close to club crest and nike logo. “La liga” patch, name, number and down sponsor from back of the shirt, perfect 2. The name on back of the baryern shirt is missing a letter, the letter N from COMAN is missing, making the shirt unwerable from my pov. The man city jersey was good, no big flaws at all

Fc24SHOP: this order was my second order from this site, in the first order was 3 players version jersey with no flaws, but this on was awful. 1. Man city shirt has the wrong premier league patch, it must be that “Champions 21/22” patch but instead is “Champions 20/21” patch. 2. Barca shirt was awful: spotify sponsor was bad printed, almost not glued at all, ready to peel off if i put it in the washing machine. The second problem with this barca shirt was the bad placement of the “la liga” patch on the sleeve. It was printed too high on the sleeve and when i m wearing it is like the patch is on shoulder, not on arm. Bayern shirt was okey, nothing to complain about.

Conclusion/opinion: i don t know if was unlucky on 2/3 orders but my money in the future goes to kittgg, the quality was very good, no flaws, stiching very ok, no big problems at all. I talked with fc24shop to solve the problem and they said on the next order, i will recieve the barca shirt for free. With jjsport i wait for the answer by the moment i write this post.

submitted by VladyComan to Soccer00 [link] [comments]
KITTGG: stitching almost perfect, no flaws on the shirts at all, very very good materials on every shirts, patches all rigtht, like the perfect order

JJSPORT: 1. the “Spotify” from the barca shirt is too high on the front side, too close to club crest and nike logo. “La liga” patch, name, number and down sponsor from back of the shirt, perfect 2. The name on back of the baryern shirt is missing a letter, the letter N from COMAN is missing, making the shirt unwerable from my pov. The man city jersey was good, no big flaws at all

Fc24SHOP: this order was my second order from this site, in the first order was 3 players version jersey with no flaws, but this on was awful. 1. Man city shirt has the wrong premier league patch, it must be that “Champions 21/22” patch but instead is “Champions 20/21” patch. 2. Barca shirt was awful: spotify sponsor was bad printed, almost not glued at all, ready to peel off if i put it in the washing machine. The second problem with this barca shirt was the bad placement of the “la liga” patch on the sleeve. It was printed too high on the sleeve and when i m wearing it is like the patch is on shoulder, not on arm. Bayern shirt was okey, nothing to complain about.

Conclusion/opinion: i don t know if was unlucky on 2/3 orders but my money in the future goes to kittgg, the quality was very good, no flaws, stiching very ok, no big problems at all. I talked with fc24shop to solve the problem and they said on the next order, i will recieve the barca shirt for free. With jjsport i wait for the answer by the moment i write this post.

2023.03.27 18:55 *Personal_Ad_7717* **45% Off Calculator, Deli Standard Function Desktop Calculators with 12 Digit Large LCD Display and Sensitive Button, Solar Battery Dual Power Office Calculator, Pink xO**

submitted by Personal_Ad_7717 to AmazonPromoCodesMO [link] [comments] |

2023.03.27 18:55 *Wild_native01* **Is there anyway to remove all the unnecessary clutter from the screen I can hardly see the map**

submitted by Wild_native01 to doordash [link] [comments] |

2023.03.27 18:55 *protawa069* **I am safe**

submitted by protawa069 to beastboyshub [link] [comments] |

2023.03.27 18:54 *ZeroFPS_hk* **[RANT] I decided to do wormhole today.**

Had a guild weekly quest plus it's been a really long time since I touched wormhole so I thought I'd give it a try. The problem is the wormhole stages are immune to 2 random elements and they don't tell you which before you enter so I had to use Lin as my main (sole) dps while bringing along Lyra, since her skill in damage resonance gives you shields which is unaffected by healing debuff nor whether or not the monsters have phys immunity. That was grindy and painful but eventually I made it to the second last stage, brought a healing weapon to heal myself back to full health, and am ready for the final stage.

I got dumpstered (stupid horse). Good news is I found out the boss isn't immune to my element (ice), so I can just switch back to my normal team and rechallenge. At the defeat screen there where two buttons, I didn't want to immediately retry (I have to change my team first) so I hit the other button.

It reset my entire fucking progress for the entire fucking chapter fml like omfg how hard is it to add a confirmation popup warning me I have to redo fucking everything from stage 1 christ jesus I casually pressed 1 button and all the grind was for nothing fucking hell with so much shit on the game lately and most of my friends and guildmates quitting the game during the past few months (which is my main motivation for the fucking game I want to interact and talk shit with other players in a multiplayer game is that too much to ask) this seemingly insignificant blunder felt like the hotta developers personally flew from china all the way to my house just to personally slap me in the face ahdbaskdvajvq3!3q92u2jrqkjbafw08rqwr ok rant over but I'm still so fucking mad. Holy.

submitted by ZeroFPS_hk to TowerofFantasy [link] [comments]
I got dumpstered (stupid horse). Good news is I found out the boss isn't immune to my element (ice), so I can just switch back to my normal team and rechallenge. At the defeat screen there where two buttons, I didn't want to immediately retry (I have to change my team first) so I hit the other button.

It reset my entire fucking progress for the entire fucking chapter fml like omfg how hard is it to add a confirmation popup warning me I have to redo fucking everything from stage 1 christ jesus I casually pressed 1 button and all the grind was for nothing fucking hell with so much shit on the game lately and most of my friends and guildmates quitting the game during the past few months (which is my main motivation for the fucking game I want to interact and talk shit with other players in a multiplayer game is that too much to ask) this seemingly insignificant blunder felt like the hotta developers personally flew from china all the way to my house just to personally slap me in the face ahdbaskdvajvq3!3q92u2jrqkjbafw08rqwr ok rant over but I'm still so fucking mad. Holy.

2023.03.27 18:54 *justhere0001* **[US][SELLING]Initial D vol 1-33, future diary vol 1-10,rurouni kenshin 1-28,Battle Royale 1-15, ultimo 1-11, steins gate Hardcover and steins gate0 vol 1-3,blood on the tracks 1-10, animal land 1-6,8,10-12, no gun life 1-13, noragami 1-23, hell's paradise vol 1-13, and Yu-gi-oh Art Box**

Hi Everyone on here today to sell some sets. Prices are negotiable let me know what your interested in and shoot me an offer!

https://imgur.com/a/xMiEp0A

Initial D vol 1-33 ( some staining on a few vols please see pics and if you would like more pics please let me know) $880 plus shipping

Yu-gi-oh Art Box $120 plus shipping

future diary vol 1-10 $325 plus shipping

rurouni kenshin 1-28 $180 plus shipping

Battle Royale 1-15 $300 plus shipping

ultimo 1-11 $250 plus shipping

steins gate Hardcover and steins gate0 vol 1-3 $160 plus shipping (all sealed in plastic still)

blood on the tracks 1-10 $50 plus shipping

animal land 1-6,8,10-12 $110 plus shipping

no gun life 1-13 $100 plus shipping

noragami 1-23 $180 plus shipping

hell's paradise vol 1-13 $100 plus shipping

submitted by justhere0001 to mangaswap [link] [comments]
https://imgur.com/a/xMiEp0A

Initial D vol 1-33 ( some staining on a few vols please see pics and if you would like more pics please let me know) $880 plus shipping

Yu-gi-oh Art Box $120 plus shipping

future diary vol 1-10 $325 plus shipping

rurouni kenshin 1-28 $180 plus shipping

Battle Royale 1-15 $300 plus shipping

ultimo 1-11 $250 plus shipping

steins gate Hardcover and steins gate0 vol 1-3 $160 plus shipping (all sealed in plastic still)

blood on the tracks 1-10 $50 plus shipping

animal land 1-6,8,10-12 $110 plus shipping

no gun life 1-13 $100 plus shipping

noragami 1-23 $180 plus shipping

hell's paradise vol 1-13 $100 plus shipping

2023.03.27 18:54 *InvestigatorAlive216* **[discussion] Double the internet and Internet Archive is still slow as dirt 😂**

submitted by InvestigatorAlive216 to macbookpro [link] [comments]

2023.03.27 18:53 *Ermete89* **I found him in a free pull. What about him?**

submitted by Ermete89 to EmpiresAndPuzzles [link] [comments] |

2023.03.27 18:53 *Spooky_Looop* **Can someone please read my chart? I’m so lost and I can’t understand myself. Hoping this helps :)**

submitted by Spooky_Looop to astrologyreadings [link] [comments] |

2023.03.27 18:52 *SamuraiJordan89* **Custom armiger build thoughts?**

submitted by SamuraiJordan89 to ImperialKnights [link] [comments] |

2023.03.27 18:52 *gamermilkyway* **local club in latvia gonna get kanye on stage**

submitted by gamermilkyway to GoodAssSub [link] [comments] |

2023.03.27 18:50 *URBANnebula64* **Python Coding: finding the energy levels of harmonic potential embedded in infinite square well**

I have tried for hours to fix this code but it doesn't work. I am trying ro generate energy level and plots of wavefunction but they are both wrong for some reason. The curves are just flat lines that spike up at the very end. This could previously worked when using a potential that just retuned 0. But when trying to use a harmonic potential it fails completely. Im not sure if my E_guess is wrong of what else might be wrong:

%matplotlib inline

# importing modules

import numpy as np

import matplotlib.pyplot as plt

m = 9.109383702 * 10**-31 # kg, electron mass

hbar = 1.054571817 * 10**-34 # J.s constant

e = 1.602176634 * 10**-19 ## A⋅s, electron charge

d = 5 * 10**-9 # side length of cubic quantum dot, meters

c = 2.998e8 # speed of light in m/s

xstart = -d/2 # start position,meters

xend = d/2 # end position

N = 2000 # number of points for Runge-Kutta

h = (xend - xstart)/N # size of Runge-Kutta steps

tpoints = np.arange(xstart, xend, h)

ϕ = 1

def V(x, potential_function):

'''

This function returns the potential energy of a particle in a 1D potential well.

Parameters:

x : float

The position of the particle in meters.

potential_function : callable

A function that takes the position x as an argument and returns the potential energy at that position.

Returns:

float

The potential energy of the particle at the given position x.

'''

return potential_function(x)

V0 = 700*e

def H(x): # define harmonic potential well

return V0*(x**2)/((d/2)**2)

def f(r, x, E, potential_function):

psi = r[0]

phi = r[1]

d_psi = phi

d_phi = ((2*m)/hbar**2)*(V(x, potential_function)-E)*psi

return np.array([d_psi, d_phi])

def RungeKutta2d(r, x, f, E, potential_function):

xpoints = []

ypoints = []

for t in tpoints:

xpoints.append(r[0])

ypoints.append(r[1])

k1 = h * f(r, t, E, potential_function)

k2 = h * f(r + 0.5 * k1, t + 0.5 * h, E,potential_function)

k3 = h * f(r + 0.5 * k2, t + 0.5 * h, E, potential_function)

k4 = h * f(r + k3, t + h, E, potential_function)

r = r + (k1 + 2 * k2 + 2 * k3 + k4) / 6

xpoints.append(r[0])

ypoints.append(r[1])

return np.array([xpoints, ypoints])

fig, ax = plt.subplots(1, 2, figsize=(12, 6), gridspec_kw={'wspace':0.3})

# Find energy states for n = 2, 3, 4, 5

ax[0].axvline(x=-d/2,c='#5f5f5f',ls='-',lw=2.5)

ax[0].axvline(x=d/2,c='#5f5f5f',ls='-',lw=2.5)

ax[1].axvline(x=-d/2,c='#5f5f5f',ls='-',lw=2.5)

ax[1].axvline(x=d/2,c='#5f5f5f',ls='-',lw=2.5)

for n in range(1,4):

# Define initial energy guesses

E_guess = (1/2 * hbar * np.sqrt((4 * V0)/(m * d**2)) * (n - 1/2))

E1 = E_guess - 3 * 10**-21

E2 = E_guess

# Solve for the first guess

# Here we set the initial conditions in a separate array r

r = np.array([0, ϕ])

psi1 = RungeKutta2d(r, tpoints, f, E1, H)[0, N]

psi2 = RungeKutta2d(r, tpoints, f, E2, H)[0, N]

### now for the secant method to converge on the right answer:

tolerance = e/100000 # set the tolerance for convergence

while abs(E2-E1) > tolerance:

E3 = E2 - psi2*(E2-E1)/(psi2-psi1)

# update initial energies for the next iteration

E1 = E2

E2 = E3

# and recalculate wavefunctions

psi1 = RungeKutta2d(r, tpoints, f, E1, H)[0, N]

psi2 = RungeKutta2d(r, tpoints, f, E2, H)[0, N]

soln = RungeKutta2d(r, tpoints, f, E3,H)

psi = soln[0][:len(tpoints)]

# Normalize wavefunction and plot result

if n % 2 == 0:

I = h*(0.5*psi[0]**2 + np.sum(psi[1:-1]**2) + 0.5*psi[-1]**2)

psi_norm = psi/np.sqrt(I)

ax[0].plot(tpoints, psi_norm, label=f'n = numerical {n}')

ax[0].set_xlabel('Position (m)')

ax[0].set_ylabel('Wavefunction')

ax[0].set_title("Even wavefunctions")

ax[0].legend()

else:

I = h*(0.5*psi[0]**2 + np.sum(psi[1:-1]**2) + 0.5*psi[-1]**2)

psi_norm = psi/np.sqrt(I)

ax[1].plot(tpoints, psi_norm, label=f'n = numerical {n}')

ax[1].set_xlabel('Position (m)')

ax[1].set_ylabel('Wavefunction')

ax[1].set_title("Odd wavefunctions")

ax[1].legend()

print(f"Energy level n = {n}: E = {E3}")

submitted by URBANnebula64 to PhysicsHelp [link] [comments]
%matplotlib inline

# importing modules

import numpy as np

import matplotlib.pyplot as plt

m = 9.109383702 * 10**-31 # kg, electron mass

hbar = 1.054571817 * 10**-34 # J.s constant

e = 1.602176634 * 10**-19 ## A⋅s, electron charge

d = 5 * 10**-9 # side length of cubic quantum dot, meters

c = 2.998e8 # speed of light in m/s

xstart = -d/2 # start position,meters

xend = d/2 # end position

N = 2000 # number of points for Runge-Kutta

h = (xend - xstart)/N # size of Runge-Kutta steps

tpoints = np.arange(xstart, xend, h)

ϕ = 1

def V(x, potential_function):

'''

This function returns the potential energy of a particle in a 1D potential well.

Parameters:

x : float

The position of the particle in meters.

potential_function : callable

A function that takes the position x as an argument and returns the potential energy at that position.

Returns:

float

The potential energy of the particle at the given position x.

'''

return potential_function(x)

V0 = 700*e

def H(x): # define harmonic potential well

return V0*(x**2)/((d/2)**2)

def f(r, x, E, potential_function):

psi = r[0]

phi = r[1]

d_psi = phi

d_phi = ((2*m)/hbar**2)*(V(x, potential_function)-E)*psi

return np.array([d_psi, d_phi])

def RungeKutta2d(r, x, f, E, potential_function):

xpoints = []

ypoints = []

for t in tpoints:

xpoints.append(r[0])

ypoints.append(r[1])

k1 = h * f(r, t, E, potential_function)

k2 = h * f(r + 0.5 * k1, t + 0.5 * h, E,potential_function)

k3 = h * f(r + 0.5 * k2, t + 0.5 * h, E, potential_function)

k4 = h * f(r + k3, t + h, E, potential_function)

r = r + (k1 + 2 * k2 + 2 * k3 + k4) / 6

xpoints.append(r[0])

ypoints.append(r[1])

return np.array([xpoints, ypoints])

fig, ax = plt.subplots(1, 2, figsize=(12, 6), gridspec_kw={'wspace':0.3})

# Find energy states for n = 2, 3, 4, 5

ax[0].axvline(x=-d/2,c='#5f5f5f',ls='-',lw=2.5)

ax[0].axvline(x=d/2,c='#5f5f5f',ls='-',lw=2.5)

ax[1].axvline(x=-d/2,c='#5f5f5f',ls='-',lw=2.5)

ax[1].axvline(x=d/2,c='#5f5f5f',ls='-',lw=2.5)

for n in range(1,4):

# Define initial energy guesses

E_guess = (1/2 * hbar * np.sqrt((4 * V0)/(m * d**2)) * (n - 1/2))

E1 = E_guess - 3 * 10**-21

E2 = E_guess

# Solve for the first guess

# Here we set the initial conditions in a separate array r

r = np.array([0, ϕ])

psi1 = RungeKutta2d(r, tpoints, f, E1, H)[0, N]

psi2 = RungeKutta2d(r, tpoints, f, E2, H)[0, N]

### now for the secant method to converge on the right answer:

tolerance = e/100000 # set the tolerance for convergence

while abs(E2-E1) > tolerance:

E3 = E2 - psi2*(E2-E1)/(psi2-psi1)

# update initial energies for the next iteration

E1 = E2

E2 = E3

# and recalculate wavefunctions

psi1 = RungeKutta2d(r, tpoints, f, E1, H)[0, N]

psi2 = RungeKutta2d(r, tpoints, f, E2, H)[0, N]

soln = RungeKutta2d(r, tpoints, f, E3,H)

psi = soln[0][:len(tpoints)]

# Normalize wavefunction and plot result

if n % 2 == 0:

I = h*(0.5*psi[0]**2 + np.sum(psi[1:-1]**2) + 0.5*psi[-1]**2)

psi_norm = psi/np.sqrt(I)

ax[0].plot(tpoints, psi_norm, label=f'n = numerical {n}')

ax[0].set_xlabel('Position (m)')

ax[0].set_ylabel('Wavefunction')

ax[0].set_title("Even wavefunctions")

ax[0].legend()

else:

I = h*(0.5*psi[0]**2 + np.sum(psi[1:-1]**2) + 0.5*psi[-1]**2)

psi_norm = psi/np.sqrt(I)

ax[1].plot(tpoints, psi_norm, label=f'n = numerical {n}')

ax[1].set_xlabel('Position (m)')

ax[1].set_ylabel('Wavefunction')

ax[1].set_title("Odd wavefunctions")

ax[1].legend()

print(f"Energy level n = {n}: E = {E3}")

2023.03.27 18:49 *Naive-Assumption-661* **Sova really said "hold my bow"**

submitted by Naive-Assumption-661 to Eggwick [link] [comments] |

2023.03.27 18:48 *URBANnebula64* **Python: finding energy level of harmonic potential embedded in an infinite square well**

I have tried for hours to fix this code but it doesn't work. I am trying ro generate energy level and plots of wavefunction but they are both wrong for some reason. The curves are just flat lines that spike up at the very end. This could previously worked when using a potential that just retuned 0. But when trying to use a harmonic potential it fails completely. Im not sure if my E_guess is wrong of what else might be wrong:

%matplotlib inline

# importing modules

import numpy as np

import matplotlib.pyplot as plt

m = 9.109383702 * 10**-31 # kg, electron mass

hbar = 1.054571817 * 10**-34 # J.s constant

e = 1.602176634 * 10**-19 ## A⋅s, electron charge

d = 5 * 10**-9 # side length of cubic quantum dot, meters

c = 2.998e8 # speed of light in m/s

xstart = -d/2 # start position,meters

xend = d/2 # end position

N = 2000 # number of points for Runge-Kutta

h = (xend - xstart)/N # size of Runge-Kutta steps

tpoints = np.arange(xstart, xend, h)

ϕ = 1

def V(x, potential_function):

'''

This function returns the potential energy of a particle in a 1D potential well.

Parameters:

x : float

The position of the particle in meters.

potential_function : callable

A function that takes the position x as an argument and returns the potential energy at that position.

Returns:

float

The potential energy of the particle at the given position x.

'''

return potential_function(x)

V0 = 700*e

def H(x): # define harmonic potential well

return V0*(x**2)/((d/2)**2)

def f(r, x, E, potential_function):

psi = r[0]

phi = r[1]

d_psi = phi

d_phi = ((2*m)/hbar**2)*(V(x, potential_function)-E)*psi

return np.array([d_psi, d_phi])

def RungeKutta2d(r, x, f, E, potential_function):

xpoints = []

ypoints = []

for t in tpoints:

xpoints.append(r[0])

ypoints.append(r[1])

k1 = h * f(r, t, E, potential_function)

k2 = h * f(r + 0.5 * k1, t + 0.5 * h, E,potential_function)

k3 = h * f(r + 0.5 * k2, t + 0.5 * h, E, potential_function)

k4 = h * f(r + k3, t + h, E, potential_function)

r = r + (k1 + 2 * k2 + 2 * k3 + k4) / 6

xpoints.append(r[0])

ypoints.append(r[1])

return np.array([xpoints, ypoints])

fig, ax = plt.subplots(1, 2, figsize=(12, 6), gridspec_kw={'wspace':0.3})

# Find energy states for n = 2, 3, 4, 5

ax[0].axvline(x=-d/2,c='#5f5f5f',ls='-',lw=2.5)

ax[0].axvline(x=d/2,c='#5f5f5f',ls='-',lw=2.5)

ax[1].axvline(x=-d/2,c='#5f5f5f',ls='-',lw=2.5)

ax[1].axvline(x=d/2,c='#5f5f5f',ls='-',lw=2.5)

for n in range(1,4):

# Define initial energy guesses

E_guess = (1/2 * hbar * np.sqrt((4 * V0)/(m * d**2)) * (n - 1/2))

E1 = E_guess - 3 * 10**-21

E2 = E_guess

# Solve for the first guess

# Here we set the initial conditions in a separate array r

r = np.array([0, ϕ])

psi1 = RungeKutta2d(r, tpoints, f, E1, H)[0, N]

psi2 = RungeKutta2d(r, tpoints, f, E2, H)[0, N]

### now for the secant method to converge on the right answer:

tolerance = e/100000 # set the tolerance for convergence

while abs(E2-E1) > tolerance:

E3 = E2 - psi2*(E2-E1)/(psi2-psi1)

# update initial energies for the next iteration

E1 = E2

E2 = E3

# and recalculate wavefunctions

psi1 = RungeKutta2d(r, tpoints, f, E1, H)[0, N]

psi2 = RungeKutta2d(r, tpoints, f, E2, H)[0, N]

soln = RungeKutta2d(r, tpoints, f, E3,H)

psi = soln[0][:len(tpoints)]

# Normalize wavefunction and plot result

if n % 2 == 0:

I = h*(0.5*psi[0]**2 + np.sum(psi[1:-1]**2) + 0.5*psi[-1]**2)

psi_norm = psi/np.sqrt(I)

ax[0].plot(tpoints, psi_norm, label=f'n = numerical {n}')

ax[0].set_xlabel('Position (m)')

ax[0].set_ylabel('Wavefunction')

ax[0].set_title("Even wavefunctions")

ax[0].legend()

else:

I = h*(0.5*psi[0]**2 + np.sum(psi[1:-1]**2) + 0.5*psi[-1]**2)

psi_norm = psi/np.sqrt(I)

ax[1].plot(tpoints, psi_norm, label=f'n = numerical {n}')

ax[1].set_xlabel('Position (m)')

ax[1].set_ylabel('Wavefunction')

ax[1].set_title("Odd wavefunctions")

ax[1].legend()

print(f"Energy level n = {n}: E = {E3}")

submitted by URBANnebula64 to CodingHelp [link] [comments]
%matplotlib inline

# importing modules

import numpy as np

import matplotlib.pyplot as plt

m = 9.109383702 * 10**-31 # kg, electron mass

hbar = 1.054571817 * 10**-34 # J.s constant

e = 1.602176634 * 10**-19 ## A⋅s, electron charge

d = 5 * 10**-9 # side length of cubic quantum dot, meters

c = 2.998e8 # speed of light in m/s

xstart = -d/2 # start position,meters

xend = d/2 # end position

N = 2000 # number of points for Runge-Kutta

h = (xend - xstart)/N # size of Runge-Kutta steps

tpoints = np.arange(xstart, xend, h)

ϕ = 1

def V(x, potential_function):

'''

This function returns the potential energy of a particle in a 1D potential well.

Parameters:

x : float

The position of the particle in meters.

potential_function : callable

A function that takes the position x as an argument and returns the potential energy at that position.

Returns:

float

The potential energy of the particle at the given position x.

'''

return potential_function(x)

V0 = 700*e

def H(x): # define harmonic potential well

return V0*(x**2)/((d/2)**2)

def f(r, x, E, potential_function):

psi = r[0]

phi = r[1]

d_psi = phi

d_phi = ((2*m)/hbar**2)*(V(x, potential_function)-E)*psi

return np.array([d_psi, d_phi])

def RungeKutta2d(r, x, f, E, potential_function):

xpoints = []

ypoints = []

for t in tpoints:

xpoints.append(r[0])

ypoints.append(r[1])

k1 = h * f(r, t, E, potential_function)

k2 = h * f(r + 0.5 * k1, t + 0.5 * h, E,potential_function)

k3 = h * f(r + 0.5 * k2, t + 0.5 * h, E, potential_function)

k4 = h * f(r + k3, t + h, E, potential_function)

r = r + (k1 + 2 * k2 + 2 * k3 + k4) / 6

xpoints.append(r[0])

ypoints.append(r[1])

return np.array([xpoints, ypoints])

fig, ax = plt.subplots(1, 2, figsize=(12, 6), gridspec_kw={'wspace':0.3})

# Find energy states for n = 2, 3, 4, 5

ax[0].axvline(x=-d/2,c='#5f5f5f',ls='-',lw=2.5)

ax[0].axvline(x=d/2,c='#5f5f5f',ls='-',lw=2.5)

ax[1].axvline(x=-d/2,c='#5f5f5f',ls='-',lw=2.5)

ax[1].axvline(x=d/2,c='#5f5f5f',ls='-',lw=2.5)

for n in range(1,4):

# Define initial energy guesses

E_guess = (1/2 * hbar * np.sqrt((4 * V0)/(m * d**2)) * (n - 1/2))

E1 = E_guess - 3 * 10**-21

E2 = E_guess

# Solve for the first guess

# Here we set the initial conditions in a separate array r

r = np.array([0, ϕ])

psi1 = RungeKutta2d(r, tpoints, f, E1, H)[0, N]

psi2 = RungeKutta2d(r, tpoints, f, E2, H)[0, N]

### now for the secant method to converge on the right answer:

tolerance = e/100000 # set the tolerance for convergence

while abs(E2-E1) > tolerance:

E3 = E2 - psi2*(E2-E1)/(psi2-psi1)

# update initial energies for the next iteration

E1 = E2

E2 = E3

# and recalculate wavefunctions

psi1 = RungeKutta2d(r, tpoints, f, E1, H)[0, N]

psi2 = RungeKutta2d(r, tpoints, f, E2, H)[0, N]

soln = RungeKutta2d(r, tpoints, f, E3,H)

psi = soln[0][:len(tpoints)]

# Normalize wavefunction and plot result

if n % 2 == 0:

I = h*(0.5*psi[0]**2 + np.sum(psi[1:-1]**2) + 0.5*psi[-1]**2)

psi_norm = psi/np.sqrt(I)

ax[0].plot(tpoints, psi_norm, label=f'n = numerical {n}')

ax[0].set_xlabel('Position (m)')

ax[0].set_ylabel('Wavefunction')

ax[0].set_title("Even wavefunctions")

ax[0].legend()

else:

I = h*(0.5*psi[0]**2 + np.sum(psi[1:-1]**2) + 0.5*psi[-1]**2)

psi_norm = psi/np.sqrt(I)

ax[1].plot(tpoints, psi_norm, label=f'n = numerical {n}')

ax[1].set_xlabel('Position (m)')

ax[1].set_ylabel('Wavefunction')

ax[1].set_title("Odd wavefunctions")

ax[1].legend()

print(f"Energy level n = {n}: E = {E3}")

2023.03.27 18:47 *jeannie15* **Nice filling salald for 548**

submitted by jeannie15 to 1200isplenty [link] [comments] |