Is there a code implementation in R of the Louvian method for community detection for a user defined modularity function?

I am working on a problem that needs to do community detection with a specific modularity function in mind. I need to have code implementation which maximizes my specific modularity function across communities. If any other methods exist, please let me know.